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ABSTRACT 


Title of Dissertation; SHIFT MEASUREMENTS OF THE STARK-BROADENED 

IONIZED HELIUM LINES AT 1640 AND 1215 X 
James Russell Van Zandt, Doctor of Philosophy, 1976 
Dissertation directed by: Hans R. Griem, Professor of Physics 

Time-resolved measurements were made of the shifts of the 
ionized helium lines at 1640 X (n = 3 -»■ 2) and 1215 X (n = 4 ^ 2) , 
and of the Stark profile of the X 1215 X line. An electromagnetic 
shock tube was used as a light source. The plasma conditions corres- 
ponded to electron temperatures of 'v3.5 eV and electron densities of 
17 -3 

0.8 to 1.8 10 cm . The measured shifts fell between two previous 

estimates of plasma polarization shifts. The measured Stark width of 
the X 1215 X line was up to 30% greater than the theoretical width. 
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CHAPTER I 


INTRODUCTION 

Spectroscopy has long been recognized as an important diagnostic 

tool for both astrophysical plasmas, where it is often the only method 

available, and in the laboratory where, unlike many methods, it does not 

disturb the plasma under study. 

The considerable theoretical and experimental efforts in this 

field have resulted in good understanding of the pressure broadening 

and shifts of spectral lines due to the Stark effect of nearby charged 

perturbers. Particular attention has been paid to the lines of 

hydrogen and the hydrogenic ions, for which the quasistatic and impact 

theories predict considerable broadening but no shifts. However, in 

1962, Berg et al r.eported^ a blue shift for the He II 4686 line, which 

they attributed to the reduction of the Coulomb potential of_the nucleus 

by the polarization of the plasma near the radiating ion. Later measure- 
2 

ments demonstrated this "shift" had been simulated by some unresolved 
Si III lines on the blue wing of the helium line. Greig et al. then 

3 

reported blue shifts of the He II 304 line. Subsequent photographic 
4 5 

measurements ’ did not verify the shift of the 304 line, but higher 
series members (256, 243, etc.) had blue shifts which could have been 
due to plasma polarization. The most recent measurement^ showed blue 
shifts for the 256 and 243 lines, with a greater shift for the 304 line, 
in agreement with Greig 's result. 

The polarization shift is expected to be important for high-Z 
ion lines and may limit wavelength accuracies in, for example, laser- 
produced plasmas. The theoretical treatments of this effect have 
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been unsatisfactoxy, * and no attempts have been made to measure 
shifts of the "Balmer" (or "second Lyman'^ series lines of ionized 
helium, at 1640, 1215, 1084 . . . S. The primary aim of this experiment 
was to look for such shifts, and investigate their possible dependence 
on plasma conditions. A secondary purpose was to measure the Stark 
broadening of the higher series members, to check the theoretical 
calculations . ^ 

A T~tube was chosen as a source because it produces a fairly 
homogeneous plasma^ near local thermal equilibrium (LTE),^^*^^ at 
a density and temperature suitable for the emission of ionized helium 
lines. The line positions were measured relative to nearby impurity 
lines. Plasma conditions were determined from photoelectric measure- 
ments of the He II 4686 line, and plasma reproducibility was checked 
by monitoring the total intensities of the 4686 line and the continuum 
near 4976 S. 

The first chapter of this dissertation has served as an introduction. 
In Chapter 2 some of the relevant results of plasma spectroscopy are 
presented. A description of the experimental apparatus and method 
appears in Chapter 3. The experimental results, with a discussion of 
them and possible errors, are in Chapter 4. 



CHAPTER II 


THEORETICAL BACKGROUND 
A. Line Intensities 


The relative intensities of emission lines depend on the population 

densities of atoms in the upper state and the probability of radiative 

transition to the corresponding lower state. 

In equilibrium, the density of ions of charge Z is related 

to the electron density and the density of atoms In the next lower 

13 

ionization stage according to the Saha equation 



Since nearly all the atoms are in the ground state, the partition 

function Z (T) can usually be replaced by the statistical weight g 

of the ground state. In this case, 1=0, and we have g = 2S+1 (1 for 

z 

He*^, and He^; 2 for and He^) . The correction AE^ ^ to the 

CO 

z-1 13 

ionization energy due to Coulomb interactions in the plasma is 


AE 


z-1 


Ze 

4irCo^D 


where is the plasma Debye length 


14 


‘d = (‘" I 


2 ^- 1/2 


( 2 - 2 ) 


(2-3) 


where is the density of particles with charge q^. Plots of helium 
Ionization stage concentrations as functions of temperature appear 
in Fig. 2-1. Plots of and other plasma properties appear in Fig. 2-2, 


3 
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Fig. 2 
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where the shaded region is typical oi T'-tubes. Note that the plasma 
3 

approximation ?•> 1 (indicating many particles in a Debye sphere) 

is only marginally satisfied. 


The population densities of the state (n,L,S) of a given 


ionization stage is given by the corresponding Boltzmann factor 


13 


\lS “ ®nLS “P 


(-¥)■ 


(2-4) 


together with the normalization condition. The exponential term is 

nearly always much less than one for excited states, justifying the 

earlier statement that most atoms are in the ground state. Note that 

an isolated atom has an infinite number of bound states, whose energies 

2 

tend to the ionization energy E^. When the atom is embedded in a 

plasma, however, the ionization energy is reduced as described above, 

and only a finite number of bound states remain. 

Treating an atom as an electric dipole radiator, the transition 

13 

probability per unit time for spontaneous emission is 


\u = 


3hc-^[4Tre„] ^ i 


av 


(2-5) 


^ 0 ^ 
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which is tabulated for many spectral lines. The sum is over the 

components of the coordinate vector of the radiating electron, and 

the average is over possible final states. Multiplying by the energy -hu 

of the photon, and using (2-4) to relate the upper state population 

density to the ground state population density N (with statistical 

S 

weight g ) we find the total power per unit volume spontaneously 
S 

13 

radiated in the given line to be 


N A, g 
g X.U u 


exp 


\ kT / ’ 


P„ = 2irhc 
Au 


(2-6) 



Since the line intensities are proportional to the concentration 
of atoms in the appropriate ionization stage, the intensity ratio of 
lines of different ionization states is an extremely sensitive function 
of temperature (see Fig. 2-3), and can be used to measure the temperature. 
Note, however, that this measurement depends strongly on the assumption 
of local thermal equilibrixim, which can require a long time and 
considerable distance to establish between states with very different 
energies. 


B. Continuum Intensities 


Plasmas emit continuum radiation due to radiative recombination 
(inverse photoionization), bremsstrahlung , and the formation of negative 
ions. A pseudo-continuum results when the Stark profiles of nearby 
lines overlap. 

The extremely weak bremsstrahlung radiation due to lon-ion and 
nonrelativistic electron-electron collisions can be neglected. That due 
to electron collisions with ions of charge z is given by 


ei 


Ibre 


3c^/6nm 


z\(m,T ) , 
3 " e 


(2-7) 


where is the free-free Gaunt factor, which is usually of order one. 

The radiation from electron-neutral collisions (approximated by 
elastic, billiard-ball type interactions) is given by 
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3c (2nm) 


( 2 - 8 ) 


When an Ion captures a free electron, the binding energy and the 
electron's kinetic energy are given to a photon. For recombination 
into a given orbital (n,L,S), the photon then has the minimum energy 
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-Ro) =* E , ” E 

Z-l,« 2-1,11 


Viewed another way, this restricts the possible final states for the 
electron for a contribution to the continuum at a given frequency. 
The recombination continuum is then, by detailed balancing. 


R 2^ / 4iij) \ y ®z-l,n . 

~ 7 , ’7372 “P (- 1 !• i~ 7 “ “z-l,n • 

c (2innKT^) \ e / n ®z,l 

18 19 

where a , is the photoionization cross section ’ > S 1 § 1 

y XI ) XI Z y 

are statistical weights, and the sum runs from the lowest allowed 
state to the highest bound state (i.e., with energy less than the reduced 
ionization energy calculated from (2-2)). 

Some electronegative atoms (H, N, 0, C, etc.) can capture a free 
electron and form a negative ion, while emitting a continuum as in 


(2-10) 


recombination. The spectral emission coefficient is, similarly,^ 


e = -y 

(j) 


„ *4 3 - /E 

' 3/2 W z (T ) (■^^J ’ 

c (2innKT^)^^ e'' V e / 


( 2 - 11 ) 


where g is the statistical weight of the negative ion (which usually has 
only one bound state), Zq is the partition function of the neutral atom, 

E IS the binding energy of the new electron (generally less than 2 eV) , 

3 . 

20 

and a is the cross section for the inverse process of photodetachment. 
This process is unimportant in hot plasmas, where the density N of 

cL 

neutral atoms is low. 

The pseudo-continuum of lines is generally important only for 

hydrogenic atoms, which are subject to the linear Stark effect, and 

then only near a series limit. The last clearly distinguishable line 

21 

of a series is then given by the Inglis-Teller limit . 

Since both line and continuum intensities increase with electron 
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density, but scale differently with temperature, the ratio of the line 
intensity to that of the nearby continuum can be used to measure the 
temperature. 


C. Radiation Transfer 


In previous sections we have discussed the spectral emission 

coefficient e of the plasma, expressed as power radiated per unit 
(1) 

solid angle, frequency interval, and volume. The experimentally 

measurable quantity is I , the power radiated per unit solid angle , 

(0 

frequency interval, and surface area of plasma observed. In the 

simplest situation, l.e., neglecting scattering, it obeys the differential 

equation^^ 

4-1 = e - k’l , (2-12) 

dx 0) (0 CO (D 

where k' is the effective absorption constant, equal to the actual 
absorption constant minus the induced emission, includes only the 
spontaneous emission. If the plasma is in LTE, the emission follows 
Kirchoff's law^^ 


e = k’B (T) , (2-13) 

OJ b) 0} 

where B is the Planck function. If we further assume the plasma to 
(1) 

be homogeneous, the solution of (2-12) is 


I (S.) = B (T)[l - exp(-k's.)l . 

0) (I) tiJ 


(2-14) 


The quantity k'£ is the "optical depth", and if k*£ « 1, the equation 
OJ w 

reduces to 


I (iO = B (T)k’£ = £ i , 
(1) (1) m u) 


(2-15) 
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as expected. In the opposite limit, >> 1, the plasma radiates 
as a blackbody. Stellar atmospheres have great optical depth at almost 
all wavelengths, while laboratory plasmas are normally optically thin 
except possibly near the centers of some resonance lines. 

D. Line Broadening 


Spectral line broadening in a plasma is a complex phenomenon, 

and no attempt is made here to discuss all the results of investigations 

j .4 22 . 23 , , 10,13 

in atomic spectroscopy, astrophysics _ . _ . . 


, and plasma spectroscopy. 


Only a physical picture of the various effects is presented. 

Let the (frequency-space) spectral line profile I(uj) be proportional 
to the light intensity between m and oj + dm, subject to the normalization 


condition 
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which satisfies the normalization condition 

00 

lf(t)l^dt = 1 , 

and has the Fourier components 

C(w) = , 

y TT LjQ-tiJ+iY 

leading to the Lorentz or dispersion profile 


(2-20) 


(2-21) 


1(03) = ic(u,)l^ = ^ (2-22) 

(mQ-o))S 

The half-half width y, the frequency separation at which the intensity 
is half the maximum, is given by the sum of the transition rates for 
transitions originating from either the upper or lower state of the line^ 


\u ^,\’u 
u 


(2-23) 


9 -1 

Since atomic excited states have relatively long lifetimes (A < 10 sec ), 

this natural broadening is almost always smaller (AX < 10 than the 

other effects we will discuss. It can of course be derived rigorously 

25 

from the quantum theory of radiation. 

When the energy levels of the radiating atoms are well separated, 

compared to mean thermal energies, electron collisions rarely exchange 

energy with the radiator , but change the polarization or phase of the 

emitted light. Although this approximation does not hold, for example, 

7,13 


for neutral helium 


(where there are nearby perturbing levels with 


the same ^n but different L) , it is well satisfied for hydrogenic atoms. 
Assuming the light to be monochromatic between collisions, we have a 
sinusoidal wave train of duration i , with the Fourier components 



13 


C^(.) 


/I f 

j ® 


i(aj 

= -5 

i((jQ-(i))/2Tr 


producing the intensity 


= 


2/1 \ 

Ln ly ((UQ-d})"? } 


If the probability per unit time of a collision is constant, 
the intervals between collisions have the Poisson distribution 


(2-24) 


(2-25) 


P(i)dr = y^e dr . 


(2-26) 


Weighting the intensities (2-25) by the corresponding probabilities, 
we again arrive at the dispersion profile (2-22) j now with width 

The frequency of light emitted by a moving atom is Doppler-shifted 
according to 


1 - 


U3 — ( 1 )- 


2 

c 




(2-27) 


Assuming the atoms have a Maxwellian distribution of velocities, 

- / m \3/2 / m v^\ _ 

f«(y)d V = (-t—rr-f exp l-qr } d v , (2-28) 


the collection will emit light with the Doppler or Gaussian line profile 

,2 

\ 

(2-29) 


13,26 


■ it ■ 




where the characteristic width is 


2Am^ 


y 2 

If 


(2-30) 


and the Doppler half-half width is /Zn4 Aoj^. 
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In contrast to the fast electron impacts, nearby ions can usually 
be considered stationary, and supply only perturbing electric fields. 

These fields perturb the energy levels (each labeled by n, L, S, and J) 
of the radiating atom and usually split them into several sublevels, 
each a linear combination of states of different magnetic quantum 
number m^ Transitions between such sublevels of different principal 

quantum nimiber give rise to the Stark components of a line. Since the 
operator expressing the interaction of the electric field 

and a given electron, has odd parity (therefore no diagonal elements), 
there is usually no first-order interaction, and second-order perturbation 
theory is used. In the hydrogenic case, however, the terms (labeled 
only by n) are degenerate, and a linear effect is found. This problem 
is most conveniently solved in the parabolic coordinates (C,n »<{>)! 


5 = r+z 
n = r-z 

(2-31) 

tanij) = y/x 

,2 2 , 2 . 2 . 

(r = x + y + z), 

27 

where the unperturbed wavefunction is 

u;/ Ir l^ ~ 2n ..m/2 m/2 , ^ 6^““^ 

T(n n m[Cni}») = e S n U U)" Cn) - 

n E + n 2 + |m| + 1 = 1,2,... (2-32) 

m = 0, +1, ..., ±(n-l) , 

27 

and the energy, correct to second order in field strength, is 
4 2 

(in units of m^e /h ) 
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( - I . 3 

2 2^2 7 . '.V 

n 


(2-33) 


- I6 \ - 3(nj^-n2)^-9m^l9]lF|2 

/j 


For a first approximation, we may assume the ions in the plasma are 

uncorrelated In this case, they produce an electric field F with the 

It), 28 


Holtsmark distribution 




( 3/2. 

exp(-x )sin 


0 


(i’') 


xdx , 


(2-34) 


plotted in Fig. 2-4, where the Holtsmark normal field strength produced 


by perturbers with density and charge is 


10 


0 


■ (fe "p) 


2/3 


(2-35) 


Integrating the energies (2-33) over the distributions (2-34) for each 
level (though the effects on the upper level usually predominate) we 
arrive at the Holtsmark profile, shown for the hydrogen line H. in Fig. 

p 

2-5. Profiles of lines subject to the linear Stark effect are usu.illy 

. . , 10. 13 

expressed in terms of the reduced wavelength separation, defined by 


a - 



(2-36) 


Note that where there is no unshifted Stark component (as in hydrogen 
transitions n=4->-2 or 3^1), the low probability of very small fields 
(since H(0) = H'(0) = 0) gives a line profile with a central dip, 
usually partly filled by other effects. 

Finally, observed profiles are broadened by the instrument response 
function of the observing monochromator. According to physical optics. 




S(!i) (statvolc cm 



Fig. 2-5 Holcsmark profile for the Line 
(broadened only by statistically uncorrelated ions) 
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It is a profile like (2-29) folded with the two rectangular slit 
functions, plus a constant background. With wide slits j Gaussian or 
triangular profrl'es are good approximations 

The profile of a line broadened by two independent effects is the 
convolution of the two profiles. 


I(x) = I^(x) ® l 2 (x) 


CO 

I, (x’)I„(x’-x)dx’ 

— 


(2-37) 


and if we assume all of these effects are independent, we may find our 
theoretical profile by convolving all the profiles; 


I, =I ..(2)1, 

theory natural electron^ Doppler ion instrument 


(2-38) 


This assumption of statistical independence is reasonable for plasmas, 
because, e.g., collisions leading to significant changes of radiator 
velocities (Doppler effect) usually involve ions whose direct contribution 
(Stark effect) is insensitive to ion and radiator velocities. 


E. Validity of LTE 

A plasma is in local thermal equilibrium (LTE) if, locally and 

instantaneously, all quantum state population densities (except for 

photon states) correspond to a system in complete thermal equilibrium 

(GTE) which has the same mass density, energy density, and chemical 
13 

composition. Departures from LTE occur when some transitions have 
unbalanced rates, so that some (generally low-energy) states are 
over- or under-populated when compared to the corresponding GTE system. 
In optically thin plasmas, where the rates of radiative excitation 
(photoexcitation and photoionization) are negligible compared to 
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those of radlat^ve de-excitation (spontaneous emission and radiative 

recombination) , the lower-energy states will be overpopulated unless 

collisionnl processes dominate radiative ones. That is, populations 

will be within ”''10% of LTE if coHibional processes are about an order 

of magnitude more important than radiative ones. Since collision cross 

sections are generally 1 irger, and energy gaps smaller, for excited 

states, LTE is most easily satisfied for them An estimate of the 

electron density required for the hydrogcnic level n to be within ' 1 - 10 % 

13 

of LTE with respect to the ion density is 


> (7 ■ lO^^cm-3) 




(2-39) 


The largest gap between atomic energy levels is generally between 

the ground and the first excited states, so the requirements for LTE 

for the ground state are usually the most restrictive. Near LTE, the 

largest transition rates are those to and from the first excited state, 

. 1 3 

and collisional rates can be expected to dominate if 

3 


N > (9 X lO^^cm ^) 
e 


\\) 


(2-40) 


It often happens that the resonance line is optically thick, so 
that radiative de-excitation of the first excited state is balanced by 
photoexcitation I’ho resonance line profile is generally dominated by Doppler 


broadening (for N sufficiently low that electron collisions cannot 

Q 


maintain LTE), so its optical depth can be estimated by 


13 


k’ d % (2 • 10 ^°cm)f,„X 
reson "x- 12 


12 \icT 


1/2 


z-1 
N ,d 
a,l 


(2-41) 


where the resonance line has wavelength absorption strength fj^ 2 ’ 

and the atoms of interest have atomic weight A and ground state density 


miGINAIi PAGE IS 
lOF POOH QUAUK'i 
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If the optical depth of the resonance radiation is greater than 

a,l 

the requirement (2-40) can be relaxed by about an order of magnitude. 
The validity of LIE for ionization stage populations an stationary 
plasmas usually need not be checked separately, since the excited 
states of a given stage are well connected with the ground state of 
the next ionization stage 

In transient plasmas, populations may depart from LTE if equilibrium 

times are long compared to the times over whicli plasma parameters change. 

The lowest transition rates for given stage usually invol“ve the 

collisional excitation of atoms jn the ground state. Assuming hydro- 

1 3 

genic behavior, the equilibrium time is then estimated by 


z-1^ (1.1 X 10 sec cm )z 

^1 " f„T N 

21 e 


-3s .3 r 


,z-l,a 


a a 




(2-42) 

2^1 cL 

where E 2 ’ is the energy of the first excited state and the term in 

brackets is the fraction of atoms or ions that must be excited into 

the next ionization stage If only partial LTE is required (1 e, 

the state with principal quantum number n is in equilibrium with higher 

13 

states) the equilibrium time is much shorter, and is estimated by 


1 V -3, 3 , „ xl/2 /2z E, 

z-1 0 / (4.5 X 10 sec cm )z / kT \ f 1 


\z E„/ \n kT / 



CHAPTER III 


EXPERIMENTAL METHOD 
A. Apparatus 

A«L T~tube and circuit . The plasma studied in this work was 

29 30 

produced in a T-tube similar to those developed by Kolb ’ and used 

3 1 36 

in several previous experiments at the University of Maryland 
37 38 

and elsewhere ’ In this device, illustrated in Fig. 3-1, an aluminum 

(alloy 2024-T4) electrode was sealed into either end of the top of a 

T-shaped tube of high-temperature glass with inside diameter of 16 mm. 

This tube was filled with the test gas at a pressure near .5 Torr (70 

Pascals). A current flowed across the 16 mm gap between the electrodes, 

ionized the gas and ohmically heated it, then returned via a backstrap 

above the T The backstrap current created a transverse magnetic field in 

the current-carrying plasma, and pressure and Lorentz force accelerated 

it down the leg of the T. This luminous front traveled 12 cm down the 

tube at several cm /p sec and struck an adjustable reflecting plate, 

where some of its directed motion was converted to random thermal 

motion. Longer expansion tubes and higher fill pressures are required for the 

formation of a separated shock, but this device produced the high 

17 -3 

temperatures (3.5 eV) and electron densities (2 '10 cm ) needed to 
excite ionized helium lines. The decaying plasma lasted approximately 
one psec 

The circuit used appears in Tig. 3-2. The relatively modest energy 
needed by the tube was supplied by a .5 pF capacitor charged to 40 kV 
(thus storing 400 J) . When charged, this capacitor was disconnected 
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BACKSTRAP 


VACUUM 

MONOCHROMATOR 



Fig. 3-1 T-tube schematic 


POWER SUPPLY 



tNOUCTANCE 


3-2 Experiment circuit diagram 
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from both the high voltage supply and ground, preventing discharges 
from either electrode to the monochromator. The high-voltage circuit 
was enclosed by a copper shield to reduce electromagnetic interference. 

To start the discharge, the nitrogen In a two-electrode pressure 
switch (initially at 30 PSI above atmospheric) was released until its 
dielectric strength was low enough for electron cascade. Since nitrogen 
was used, no ozone or nitrogen oxides were formed, as in a discharge 

I 

in air. The poor control over discharge timing was no problem, since 
the discharge itself triggered the recording system. 

The measured quarter-cycle time was .675 psec, indicating a total 
circuit inductance of 370 nH. A carbon resistor of about .01 damped 
out the oscillations after two cycles. 

The vacuum system is shown in Fig. 3-3, During the experiment, 
valve V3 was closed, while shut-off valve VI and leak valve V2 were opened, 
so the test gas flowed from the inlet, through liquid nitrogen cold trap 
CT3, into the T-tube. It then leaked into the monochromator through 
entrance slit SI, and was removed by pumps DPI and MPl. V2 was adjusted 
so the leak rates into and out of the T-tuhe balanced, and the pressure, 
measured by thermistor gauge G2, stayed at the desired value. 

Between experimental runs, the T-tube was isolated by closing slit 
valve V9 and shut-off valve VI, and kept clean by the small diffusion 
pump DP2. Cold trap CT2 was cooled by a conventional refrigeration 
system and valves V3 and V4 were solenoid-controlled, so this secondary 
pumping system could operate unattended. Since the small pump was not 
forced to pump through a slit, it proved more effective than the large 
pump at outgassing the T-tube and associated plumbing. 
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G1 ionization gauge Veeco RG-83 
G2,G3 thermistor gauge CVC GT-340A 
G4,G5 Gauge thermocouple 

G6 Gauge cold cathode discharge 
VI Metering valve, 1/4 in. 

V2 Screw valve Veeco, 3/8 in. 

V3 Solenoid valve Veeco, 3/4 in. 

V4 Solenoid valve Veeco, 3/4 in. 

V5 Gate valve 

V6 Solenoid valve 

V7 Solenoid valve 

V8 Entrance slit valve 

V9 Exit slit valve 

VIO Air inlet valve 

DEL Diffusion pump, NRG, 6 in. 

DE2 Diffusion pump, 2 in. 

Cn Liquid nitrogen cold trap 

C12 Freon cold trap 

CT3 Liquid nitrogen cold trap 

MPl Fore pump DuoSeal 1397 

MP2 Fore pump DuoSeal serial 16025-2 

51 Entrance slit 

52 exit slit 

PM photomultiplier tube 

Fig. 3-3 Schematic of vacuum system 
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A. 2 VUy monochromator and detector . The optical arrangement is 

shown in Fig. 3-4. A McPherson 225 one-meter monochromator scanned the 

ultraviolet lines shot-to-shot. Its 50 p entrance slit was flush with 

the wall of the T-tube, about .5 ram from the reflector. Since the plasma 

conditions changed sharply as the reflector was moved, the position was 

chosen which gave the most reproducible plasma. A 1200 iLnes/mm Pt-coated 

grating, with speed about f/13.6, focused the light onto a 30 p exit 

slit, for a measured reciprocal dispersion of 8.3 X/mm (4.2 X/iran in 

second order) and an approximately Gaussian instrument response function 

of width 'U.41 X ('v- ,19 X in second order). The light then fell on a 

39 

p-terphenyl coated disc, causing it to fluoresce. These visible photons 
left the vacuum chamber through a quartz window and were detected by an 
EMI 6522 photomultiplier. For some work, a 2 mm thick Mgp 2 filter was 
placed between the exit slit and the fluorescent screen to remove light 
from second order, since it transmitted 40% of the light at-1215 X but 
essentially none below 1100 X.*^^ The exit slit, screen, and PM tube 
were replaced by a film holder for photographic work. The instrument 
function and wavelength calibration were checked using a low-pressure 
Tanaka lamp. 

A. 3 Visible monochromators and detectors . For diagnosis of the 
plasma conditions, three Jarrell-Ash visible-light monochromators were 
used One 1/ 2-meter focal length monochromator, with instrument width 
.4 X, scanned the He II 4686 line shot-to-shot to determine the 

electron density (from line width)^*^ and temperature (from line; 
continuum ratio)^^ The reproducibility of the plasma was monitored on each 
shot by two 1/ 2-meter instruments, one for the continuum at 4976 X 
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(sensitive to election density), and one for the He II 4686 line 
(sensitive to temperature, and used for later data processing). 

PM tube response was checked using neutral density filters and 
pulses from a light emitting diode, and was found linear for signals 
of up to .2 V (1.1 mA) with a PM supply voltage of 900 V. 

Each PM tube housing was insulated from its monochromator, and 
signals were taken from both the anode (negative pulse) and last 
dynode (positive pulse), carried by shielded, coaxial cables terminated 
by 90 resistors, subtracted to suppress noise, amplified, digitized, 
and stored electronically For details on the waveform recorder, see 
Appendix A. 


B. Data Reduction 


The best-fit values of the four parameters (line intensity I, 

line position background intensity B, and electron density N^) are 

found using the following procedure. Assume we have the n measurements 

y (X.) and the corresponding theoretical intensities T = -—-T 
11 ^ 0 

where T(a) is the theoretical profile after convolution with the instrument 
profile G(a) 



T(a) 


09 

f 


S(a“a* )G(a' )da' 


—CO 


(3-1) 


and the instrument function has been transformed into a-space. The best- 
fit values minimize the sum 


2 

a 


1 

n-4 


I [y^ - (IT^+B)]^ , 


(3-2) 


giving the conditions 
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so I and B are found by solving the linear system 



(3-3) 


The computer program "guesses" an electron density to use for the 

transformxng of the instrument function, convolves the theoretical and 

2 

instrument profiles, then finds a from (3-2) (subject to (3-3)) for 

many values of N and When the best values are found, the new N 

e 0 e 

is used to again transform the instrument function. The entire convolution 
and fit are repeated until successive values of are sufficiently close, 
e.g., within 2% of each other. A general discussion of least-square fitting 
when the functional parameters do not occur linearly (e.g., and N ) 

U 6 

appears as Appendix B. Details on the computer programs appear in 
Appendix C. 



CHAPTER IV 


RESULTS AND DISCUSSION 
A. Results 

Examples of photoelectric measurements of the emission profiles 

of the ionized helium lines at 4686, 1640 and 1215 R are shown in 

Figp,4-1 through 4-3, In each case, the solid line is the best-fit 

9 10 

theoretical curve of Kepple, ’ convolved with the Instrument 
profile (taken to be Gaussian) , Dashed lines are the best-fit 
continuum levels, determined primarily by points far from line center, 
which are not shown. Crosses represent points not used in the best-fit 
procedure. 

The 4686 line was found to be unshifted, as in previous experi- 
2 

ments. Its profile was in good agreement with theory, and the 
plasma electron density and temperature were deduced from -its width 
and line. continuum ratio, respectively. 

The position of the 1640 line was measured relative to the 
A1 II 1670 line, and a fairly constant red shift of .11 X was found. 
These shift measurements can be found in Table 4-1 and Fig. 4-4, 

No conclusions could be drawn about the Stark width of this helium line, 
because the observed profile was dominated by instrument broadening. 

The relative positions of the He II 1215 and Si III 1210 lines 
were measured photoelectrically. The helium line was found to have 
a red shift of approximately .19 8., increasing as the density and 
temperature fell at the end of the discharge. The halfwidth of the 
1215 line was also determined as a part of the best-fit procedure. 


29 



Relative Intensity 



Relative Intensity 



Fig. 4-2 


Measured and best-fit profile for He II ^ 4686 line 
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Fig. 4-4 Estimated and measured shifts of He II X 1640 S line 
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Table 4-1 Plasma Conditions, Shifts of Hell \ 1640 % 
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These data are shown in Fig. A-5 and Table 4-2. 

B. Discussion of Possible Errors 

B.l Impurity Lines . Photographs of spectra near each of the 
helium lines showed many Si, 0, and A1 lines. The Jarrell-Ash 1/2-m 
monochromator could easily resolve the Si III and 0 II lines near 
He II 4686, and photoelectric scans were made using points between 
these impurity lines (see Fig. 4-6) . 

A survey spectrum was taken near the 1640 line using Kodak 
SVTR film in the camera attachment for the McPherson 225 vacuum 
monochromator (see Fig. 4-7). Many Si, 0, and A1 lines were identified, 
in both first and second orders. Fortunately, none of these obscured the 
1640 line. The nearby A1 II 1670 line, chosen as the wavelength 
standard for position measurements of the 1640 line, was partially 
obscured by second order lines of 0 II and 0 III. Photographs using 
an MgF^ filter were then taken, which showed no further problems with 
impurity lines. To eliminate second order lines during photoelectric 
scans, the filter was placed between the exit slit and the scintillating 
disc 

A photographic spectrum near 1215 X showed many 0 II, 0 III, 0 IV 
Si III, and Si IV lines, including the second order 0 IV 608 line on 
the red wing of the helium line (see Fig 4-8) To eliminate these, the 
MgF^ filter was again used Cor both photographic and photoelectric runs. 

The resonance lines of N II at 1084 X prevented any observation 
of the next member of the He series, while the He II 1025 line proved 
too weak for reliable observation. 





Fig. 4-5 Estimated and measured shifts of He II A 1215 X line 
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Shifts and Widths of Hell X 1215 X 


Table 4-2 Plasma Conditions, 
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Fig. 4-6 DensitonJeter scan of spectrum near He II X 4686 S line 


37a 








38 


B«2 Wavelength Standards . All line position measurements were 
made relative to nearby impurity lines, and the accuracy of this pro- 
cedure had to be verified. The Stark shifts of these ion lines are 
expected to be small^ (just as their widths are small), but a 
plasma polarization shift certainly cannot be ruled out a priori. 

To check for such shifts, several line position measurements were made 
on a Grant comparator-microphotometer . The second-order lines were 
found to be shifted with respect to the first-order lines by .10 R, 
but otherwise, shifts were less than the measurement accuracy of 
.05 R. This is consistent with previous measurements, in which no 
shifts were found for the 0 III and N III lines near 300 R. In 
photoelectric (time-resolved) studies, no absolute shifts of the 
reference lines were measured as the plasma cooled, also arguing 
against substantial absolute shifts. Only the statistical errors 
In the measured shxfts are indicated in the tables and figures 

The monochromator wavelength scale was checked by measuring 
photographically the wavelength displacement between settings corres- 
ponding to the centers of the helium and reference lines. The errors 
in both cases were less than the setting error of 0.02 R 

B.3 He II 1215 Asymmetry . The helium 1215 line was expected 

to have a symmetric, double-peak profile (like that of H ), but 

p 

photoelectric scans showed only the peak on the blue side (see Fig 
4-3). This was interpreted as showing reabsorption by hydrogen in a 
cooler boundary layer, since the hydrogen Lyman-ot line lies .50 R 
to the red of the (unshifted) helium line center. To check this 
explanation, two scans were made, using mixtures of helium plus 0 5% 
hydrogen, and helium plus 1.0% deuterium, respectively. The amount of 
absorption increased with the increasing admixture of hydrogen, and, 
in die case of the deuterium, the dip shifted to the blue, as expected 
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The residual concentration of hydrogen was estimated from these runs 
to be approximately 0.2% Since natural, Doppler, and Stark broadening 
are all very small for the hydrogen line (<.l S) , points near the 
dip were merely excluded from the fitting procedure. 

B.4 Departure from ITE . Temperature determination from a 

helium ion line: continuum ratio requires that LTE holds also for the 

ion ground state populations, so that the line Intensity (proportional 

to the population in the excited state) and the continuum intensity 

(due mainly to recombination radiation) both have their equilibrium 

values. The equilibration time for atomic states can be estimated 

from (2-42) to be only a few nanoseconds, for both neutral and ionized 

helium On the other hand, the recombination times (into the ground 

43 44 

states) are estimated to be ’ 2 psec for formation 

of singly ionized helium and 20 psec for neutral helium Singly and 
doubly ionized states are, then, expected to be overpopulated, 
simulating a temperature higher than the true electron temperature. 

For the validity of complete LTE in a stationary plasma with 
temperatures near those in the experiment, Eq. (2-11) gives an 
optical depth of ''-ISO, for the resonance line (He II X304 8). We 
are thus justified in relaxing (2-40) by an order of magnitude, 

18 ■*“3 

and the electron density required for complete LTE is N n,i.4xl0 cm , 

0 

which is not reached in the experiment. On the other hand, the 
requirement (2-39) for partial LTE for the level n=4 (upper state of 
the 4686 8 line) is easily satisfied. 

Since the actual electron density is about an order of magnitude 
lower than that required for complete LTE, and the continuum intensity 
is proportional to the electron density while the line intensity 
IS not, we estimate that the line: continuum ratio may be too high by 
an order of magnitude, compared with the LTE value at the true 
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temperature. This yields a temperature ('''3.5 eV) that is too 
high by about ,5 eV. Similarly, if the neutral excited state 
populatxon density were too low by an order of magnitude, tae intensity 
ratio of an ionized and a neutral line would overestimate the temperature 
by about .5 eV. A measurement of the intensity ratio of the He II 4686 
and the He I 3889 lines was performed, yielding temperatures near 
4.1 eV. Since the two effects (overpopulation of singly ionized 
states due to recombination relaxation during the rapid cooling, 
and overpopulation of excited states of He II due to low collision 
rates) are additive, the true electron temperature is estimated to be 
less than the lower figure by '''20%, i.e., near 3.0 eV. 

OQ 

A previous measurement^ of the absolute intensity of the 

17 -3 

He II 4686 line in a shock-tube plasma at N '\' 10 cm indicated the 

populations of the lower excited states of the ion deviate by perhaps 

a factor of 4 from LTE. However, measurements of temperature in the 

same experiment by Thompson scattering of laser light (which does not 

depend on LTE for atomic states) and the intensity ratios of the 

« 

He II 4686 and He I 5876 lines showed good agreement. 

> 

B.5 Summary of Errors. Possible errors in the determination of 

t 

electron density were judged to be 5% due to statistical fluctuations 
and 10-15% due to theoretical uncertainties.^*^ Errors in temperature 
measurements were estimated to be .1 eV statistical and .2 eV 
theoretical (after applying the 20% correction). These possible 
diagnostic errors were not judged to endanger the principal conclusions 
of the work. The tables and figures indicate only statistical errors. 

Errors in the measurements of the shifts were .05 X or less due 
to statistical fluctuations. Systematic errors due to the shift of the 
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reference lines could not be ruled out, but were shown to be less 
than .05 S and are e<potted to be smaller. 

C. Discussion of Results 


As mentioned in the Introduction (Chapter I), previous shift 
measurements of He ion lines have concentrated on the Lyman-series 
lines (n, =1). In principle, these measurements can be used to 

calculate the energy level perturbations, and the shifts of the 
"Balmer"-series ] ines can be found in turn. Since the agreement 
between the various measurements is so poor, little is learned in 
this way 

The polarization shift is difficult to treat theoretically, 
and only estimates have been made thus far. Conceptually, the 
radiating ion is expected to attract plasma electrons, which partially 
screen the nuclear charge seen by the optical electron. A simple 
classical argument^ gives the wavelength (or wavenumber) shifts of 


the Lyman-series lines to be 

39 2 

4X t.'H 8 Vo" <" , V . 

17 ■ - 3 ’ 4 • 

0 V z 

2 2 

where a^ is the radius of the first Bohr orbit: a^ = Ti /me , and 

V is the interaction energy between the perturbing plasma electron 


and the radiating ion. Since the wave packet of the perturbing 

electron will be comparable in size to the atom, Griem proposes^ 

2 

to use the averaged interaction V=e /r, where r is the characteristic 

2 

distance between the nucleus and the optical electron. r=n Sq/z • 

6 2 
Neiger proposes the modified formula V=(J/2)c /r, which ls tlie 


electrostatic energy of a uniform sphere of charge e and radius r 
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in the field of an equ^l but opposite charge at its center. Burgess 

and Peacock argue that the density of electrons near an ion is 

low enough that their velocities are not in equilibrium with the 

surrounding plasma, being directly related to their electrostatic 

energies. They suggest using the interaction energy at the average 

2 1/3 

perturber-perturber distance, V=e . Note that all these estimates 
predict blue shifts (for the Lyman-serles lines) proportional to N^, 
but decreasing with temperature (since, at high temperature, the 
electron's thermal ’ energy is large compared with the electron-ion 
interaction energy, and it doesn't see the potential well) Denoting 
by the chosen interaction energy when the optical electron has prin- 
cipal quantum number n, and expressing the unperturbed energy levels 
in term of the Rydberg constant R, we find, for the wavenumber shifts 
of the "Balmer"-series lines, 

Av = I n R |(n^-l)exp - (Z^^-Dexp . 

2 

This can be converted to a wavelength shift by multiplying with A^, 

or an energy shift by multiplying by he. Shifts predicted by each 

2 2 2 2 2 1/3 

of these choices for V (ze /n a^, 3/2 ze /n a^, and e ) ^re 

plotted in Figs. 4-4 and 4-5. For both lines, Burgess and Peacock 

predict very small blue shifts, nearly independent of temperature. 

Griem's estimate gives somewhat larger shifts, while the stronger 

interaction proposed by Neiger gives large shifts with strong temperature 

dependence. Using the measured values of the temperature, the data are 

consistent with an interaction energy between those of Griem and Neiger, 

while Burgess and Peacock underestimate the shifts To illustrate the 

effect of the systematic error discussed above in the temperature 
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measurement, the shift predicted by Griem'a formula was recalculated 

using a 20% lower temperature, the results being shown as the dashed 

curve in Figs. 4-4 and 4-5. After this correction, his interaction 

energy gives the best fit to the data. 

The halfwxdth of the 1215 X line was up to 30% greater than that 
9,10 

calculated by Kepple. This is to be compared to a previous theta- 

46 

pinch experiinent, in which the ratio of the widths of the 4686 and 
1215 8 lines agreed with the calculated value. However, this 
experiment was done at a substantially higher temperature, T >10 eV, 
so that the difference may not be significant. 

D. Conclusions and Suggestions 

Shifts have been measured of the first two lines of the "Balraer" 
series of ionized heliimi. They are consistent with a plasma polarization 
shift, where the interaction energy between the radiating ion and the 
plasma electrons is between chose proposed by Griem and Neiger and 
probably closer to the former. 

The Stark width of the 1215 X line of ionized helium has been 

measured, and found to be up to 30% greater than calculated by 
9,10 

Kepple, and increasing as the temperature and density of the 

plasma decreased at the end of the discharge. This is perhaps 

due to an increased interference by the 1215 R line of hydrogen. 

Further studies of the plasma polarization shift might include 

more careful measurements of shifts of the hydrogenic spectra of 

39 

heavier atoms, e.g., C VI 33.8 A. Previous measurements showed no 

shifts, but with a possible error of ,05 R, (In this connection, it 

47 

is interesting to note that measured center wavelengths < e.g.. 
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of helium-like copper (Cu XXVIII) are slightly below theoretically 
predicted values.) An attempt might also be made to observe shifts 
of the higher ”Balmer'*-series members of ionized helium, perhaps in a 
Z-pinch or 0-plnch, with their greater optical depth. 



APPENDIX A 


WAVEFORM RECORDER 

To reduce the error and delay of manual data taking with the usual 
Polaroid oscillographs, a waveform recorder was designed and built for 
this experiment (Fig.A-1). The signal from one of the PM tubes is 
amplified and applied simultaneously to 31 comparators. A voltage 
divider provides reference voltages for the comparators, so for a 
given signal voltage some of the comparators will be "on" and the rest 
*'off". Integrated circuits accept the output of all the comparators, 
count the number "on", calculate the corresponding 5-bit binary number, 
and store it in a 5 bit by 64 word random access memory. When triggered, 
control circuits advance the memory address counter and give write 
commands once every 100 nanoseconds (or selectable, slower rates) for a 
total of 64 cycles. It then switches to "playback" mode, supplying the 
stored numbers, each in turn, to a digital-to-analog converter. This 
analog signal is a reconstructed version of the original signal, and can 
be displayed on an oscilloscope. 

The recorder consists of five such analog-to-digital converters 
and memories, plus two digital-to-analog converters, so 5 signals can be 
recorded, then any two displayed simultaneously. 

If the waveform is acceptable, the investigator may set the twelve 
"fixed data" thumbwheel switches and initiate recording. The shot 
number (incremented each time the device is triggered), the fixed data, 
and the contents of all 5 digital memories are written to a 9-track 
magnetic tape for later computer processing. The waveform recorder then 
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Fig, A-1 Block diagram of waveform recorder 
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reverts to "ready" mode, waiting for the next trigger pulse If the shot 
was unacceptable (due to switch misfire or abnormal time history of a 
monitor signal, for example), recording can be bypassed. 

Details on operation procedures and performance specifications of 
the wavefoTnn recorder appear in the following instruction sheet. 



Digital Data Acquisition System 
Instruction Manual 
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1. General Information 

The Digital Data Acquisition System (DDAS) is a high speed analog 
to digital converter and memory. It can record 64 data samples on each 
of 5 channels, with a sample interval as short as 100 nsec. These 
stored samples can be displayed on an oscilloscope and recorded onto a 
9-track magnetic tape. 


2 . 


Technical Specifications 
sample rate, once every 
internal amplifier risetime 
useful signal range 

maximum signal range 
resolution 
channels 

signal input impedance 
trigger level 
max trigger signal range 
trigger input impedance 
playback sweep output 
analog output 
enabling circuit 

mating input amplifier 


mating digital tape deck 
magnetic tape record 


• 1, 1*> ^*9 10, or 20 []sec« 

80 nsec 
0 to + 32 V 
-1 to +1 V. 

3.1% of full scale 

5 

50 

+1.1 V 
-.6 to +5 V. 

IMfi 

22.7 Hz sawtooth, 0-2.6 V 
0-5 V 

enabled if external circuit resistance 
IS less than 100 

Tektronix type 127 preamp power supply, 
with matching Tektronix oscilloscope 
preamp . 

Cipher model 70M-360, producing 800 BPI, 
9-track, IBM-compatible magnetic tapes. 
329 bytes of 8 bits each. . . 

6 bytes (BCD, 2 digits/byte) fixed data 
from thumbwheel switches 

3 bytes (BCD, 4 low order bits) experi- 
ment count 

320 bytes (binary, 5 low order bits) 
data, grouped by time 
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3. Installation 

For optimum protection against radio frequency interference, the 
unit should be mounted in a shielded 19 in. relay rack Several inches 
clearance below the unit are necessary for ventilation. 

4. Operation 

The Cipher tape deck should never be switched on unless the DDAS is 
on, so the proper logic inputs are provided. 

1. Turn on the DDAS and associated preamplifiers. Allow preamps to 
warm up 

2 If a tape is desired, turn the tape deck on and load a tape 
The ’’RECORDER READY" lamp should light. 

3 Switch the operating mode to "AUTO SEQUENCE", switch to 
"TRIGGER ENABLE INT", and press "tlECORD' BYPASS" 'Jhe "ENABLED" 
lamp should light 

4. Switch "DISPLAY CHANNEL SELECT" to "1". The "A-D DISPLAY" 
lamps are now displaying, in binary digital form, the signal on 
channel 1 . 

5. Ground the channel 1 preamp input. Advance the preamp "vertical 
position" control until all display lamps are lit. If this 
cannot be done, adjust the 127 preamp power supply "DC level" 

(on top of case) . 

6. Back off the "vertical position" control until all lamps just go 
out. The zero level is now adjusted Repeat steps 4-6 for the 
remaining channels now, and frequently during the experiment 

7 Connect the trigger and signal cables If an "enable" circuit 

cable IS to be used, connect it and switch to "TRIGGER ENABLE EXT" 
Set the desired sampling interval. When triggered (by a signal 
or by using the "MAN TRIGGER" button) the unit will record its 
64 samples of each channel and increment the "EXPERIMENT COUNT". 

8. If a visual monitor is desired, connect the "PLAYBACK SWEEP" to 
the "EXT HORIZ IN" jack of an oscilloscope, and one or both of 
the "ANALOG OUTPUT" 's to the vertical amplifier inputs. Set 
"ANALOG CHAN SELECT" to the desired channels. 

9. When a signal is recorded, the unit will automatically switch to 
playback mode, the corresponding mode lamp will light, and the 
stored waveforms will be displayed on the oscilloscope. 
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10. If a recording is desired, set the desired "FIXED DATA", and 
press "RECORD DATA". Otherwise, press "RECORD BYPASS". The 
unit is again ready to record a set of signals The unit may 
be switched to "MAN PLAYBACK" to again display the recorded 
signals 

11. After experiment has been completed, press "EOF" several times, 
and rewind and unload the tape 

12. Turn the tape deck off, then the DDAS and other equipment. 

Alternate operating modes are provided for diagnostic purposes. In 
"SINGLE STEP PLAYBACK" mode, the contents of one word in memory, 
corresponding to the "DISPLAY CHANNEL SELECT" setting and the octal 
address shown under "MEMORY ADDRESS", are displayed under "MEMORY DISPLAY" 
and appear at the "ANALOG OUTPUT" jacks. The associated pushbutton steps 
to the next sample. 

In "MAN SAMPLE" mode, the unit stores samples one at a time, when the 
"SAMPLE STROBE" pushbutton is pressed. The unit must be enabled and 
triggered before sampling can begin. 

In "CAL" mode, the analog to digital converters operate continuously 
and any of them can be displayed on the "A-D DISPLAY" lamps. 



APPENDIX B 


STATISTICS 

In most experxments the investigator assumes a functional form 
governing his data which has several parameters, and the object of his 
experiment is to determine the values of the parameters. If there is only 
one parameter, the quoted result might be 

* 

a = a 1 o , 

where a is the true value (usually unknown), a is the "best" value which 

* 

can be determined using the data, and o indicates the error in a . We 
usually mean by a the mean square deviation of the data from the best 
value 

2 * 2 

o = (x^-a ) , (B-1) 

/ 

where the are the results from several similar experiments. It is 
necessary to extend this to the case of several parameters and specify 
a way of calculating the quoted values. 

Assume the functional form is 


y = f(^,x) , 


(B-2) 


where x is the independent variable, y the dependent variable, and the 

48 

are parameters We define the error function 

M(a) = I 


(B-3) 


k=l 0 (Xj^) 


and let the "best" ^ be that value ^ which minimizes M We find it by 
solving the set of m equations 
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= 0 


(B-4) 


ir^ 


The errors in these parameters are given by the elements of the 
variance-covariance matrix^® 


which can be calculated from 
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(B"5) 


-1 1 9^M(a) 

a, = (H H, = 4 T-- 

ij J-J ij 2 3a^Da^ 


(B-6) 


The variance of one of the parameters is then snd the 


correlation matrix is 
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a 

C 

11 a o 
1 J 

If all the a (xj^) have a common value o, the solution of (B.4) is 
independent of that value. After this least square solution is found 
a can be calculated using 


(B-7) 


2 

a 


1 

n-m 


n 


I 

k=l 




(B-8) 


where we divide by n-m because after the parameters aj^..,a^ have been 
calculated from the data, only n-m degrees of freedom remain. 

If f(a;x) IS linear in its parameters, the calculations are, of course, 
much simpler, since (B~4) is then a linear system which can be solved 
exactly. Failing this, a search must be performed in ^ space for the 


best value. 



APPENDIX C 


PROGRAMS 

The data read from the waveform recorder tapes are processed by 

several programs, each accepting an input file plus control or 

data cards , and producing one or more output files The last programs , 

PROFILE, VPLOT, and THEORY, also print their results. Other programs 

are available to read and list each file for debugging. All mainline 

programs were written in FORTRAN for use on a Univac 1108 computer with 

the EXEC-8 operating system. Intermediate files are "direct-access" 

files on disc or drum storage, like those developed by IBM for their 
50 

computers, but not defined within ANS FORTRAN. Other nonstandard 
features used include PARAMETER statements and FORTRAN procedures. 

The first program, REVERT, uses the assembly-language subroutine 
TREAD to read the 9-track tape produced by the waveform recorder. The 
tape record format is shown in Fig, C-1 REVERT assumes the scale settings 
of the input amplifiers and the sample rate of the recorder were set on 
the "fixed data" thumbwheel switches. The alphanumeric file header (a 
prose description of the run), number of channels used, and wavelength 
for each channel and shot number are read from cards. The file header is 
written into the output file, copied by later programs, and identifies 
all printed output. Specified shots may be dropped at this point. 

Since the waveform recorder stores 6.4 psec of the signal, while 

the plasma lasts only about one psec, REVERT tries to select only the 

useful part of each signal. The first twelve records are read, the average 

time T of the maximum of the monitor signal is found, and the tape is 
max 
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rewound. Each record is then read, and the data for eight samples, 
starting at time T , are scaled and written to the output file, with 

D13X 

format shown in Fig. G-2. An end-of-file marker is written after the 
last record. 

During the experiment, the light is sometxmes attenuated to prevent 
PM tube saturation, and PARAM corrects the measured intensities to 
account for this. Since PROFILE requxres that the monitor sxgnal be 
strictly decreasxng, PARAM also chooses a decreasing portion of each 
signal and discards the rest. The output record format is shown in Fig. C-3. 

BSORT sorts the records, first on wavelength, then on shot number. 
Experimental points can be taken in any order, but in this step all data for 
a given wavelength are collected. The format of the records is unchanged 
by BSORT. 

PROFILE unfolds the data, recorded as intensity as a function of 

time at different wavelengths, into intensity as a function of wavelength 

(a line profile) at different times. Since the ionized helium line 

intensities are sensitive to temperature, all data for one profile must 

be taken under the same plasma conditions. PROFILE does this by taking 

all the data for equal monitor signal (from the total intensity of the 

He II 4686 line) . The time at which the monitor signal decays to this 

level is found, and the shot is discarded if this time is further than 

1.73 standard deviations from the mean. Similarly, any intensities at a 

given wavelength which differ from the mean by more than 1.8 standard 

deviations are discarded. Profiles are then found for successively 

lower monitor intensities (therefore later times) . The means and 

standard deviations of intensities at each wavelength go to one file 

(shown in Fig. C-4), which VPLOT uses to make a printer-plot of the line 

( 
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profile. All undeleted data points are vrritten to a second file (shown 
in Fig. C-5), used for fitting. 

The actual least-squares fit is done by THEORY. As described in 

the section on data reduction, the convolution of the theoretical line 

profile with the instrument response function is done first, in alpha 

space, using an assumed electron density. The instrument function is 

assumed Gaussian, so the convolution integrals are done using the 

52 

Gauss lan-Hermite '3-point quadrature formula. This profile is fit to 

the experimental data and a new electron density is found. The convolution 

and fit are repeated until the electron density converges, usually within 

four Iterations. Each of these fits requires a search for the values 

of the four parameters (line intensity X, background intensity B, line 

center X , and electron density N (line width)) that minimize the 

2 

mean square deviation o of the fitting function from the experimental 

points. The subroutine ZXPOWL, from the International Mathematical 

53 

and Statistical Library (IMSL) uses the function-minimization algorithm 

54,55 

described by Zangwill ’ to find the best-fit values of ^nd N^. 

For each trial values of and it calls the subroutine FUNCT3, 

which in turn calls other subroutines to calculate the best values of 

2 

the two linear parameters, and the corresponding a , using standard 
methods. 

V/hen the best values of all four parameters are found, subroutine 

2 

FUNCT2 finds the second derivative matrix of a numerically, inverts it, 
and normalizes it to get the standard deviations and correlation 
matrix of the best-fit parameters. If the line is He II 4686, it uses 
the line: continuum ratio to calculate the plasma temperature Sub- 
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routine TPLOT plots the average of the experimental points at each 
wavelength, the best-fit theoretical profile, and the background level. 
The entire procedure is repeated for each profile, but since the line 
center and electron density are carried over each time, subsequent 
fits converge rapidly. 


6 bytes 

3 bytes 

320 bytes 

FIXED 

COUNT 

DATA 


FIXED 6 bytes BCD, 2 characters/byte data from thumbwheej 

switches 

COUNT 3 bytes BCD, 1 character /byte shot number 

DATA 320 bytes binary number, 1/byte data, grouped by time 

total, 329 8-bit bytes/record (excluding parity and check bits) 

Fig. C-1 Record format of waveform recorder tape 


2 words 


1 word 

1 word 

16 words 

LABEL 


SCALE 

COUNT 

’’l ^1 ^^2 ^2 • * • "^8 ^8 

LABEL 

2 

words, FIELDATA 

First 5 characters are the wavelength in 
A (decimal point assumed before last 
character). Next 3 characters are the 
shot number 

Amplification on preamplifier (V/div) 

SCALE 

1 

word, real (R) 

COUNT 

1 

word, integer 

Shot number (same as above) . 

^i 

^i 

total: 

1 word , R 
1 word, R 
20 words/ record 

Time of sample (iisec. after trigger pulse) 
Signal amplitude (V) 


Fi-S* C-2 ‘Format of data record written by REVERT 
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2 words 

1 word 

1 word 

1 word 

8 words 

8 words 

LABEL ' 

SCALE 

COUNT 

POINTS 

V 2-'-^8 



LABEL, SCALE, COUNT, T , Y, as before 

' ’ i’ 1 

POINTS 1 word, integer Number of data poxnts (always 8) 
total. 21 words/record 

Fxg. C-3 Format of data record written by PARAM or BSORT 


1 word 

36 words 

36 words 

36 words 

MONITOR 

WAVELENGTH 

AVERAGE 

SIGMA 

MONITOR 

1 word, R Intensity of monitor for 

this profile 


WAVELEKCTH 36 words, R Wavelengths (8.) 


AVERAGE 36 words, R Average of signal intensities at corresponding 

wavelength . 

SIGMA 36 words, R Standard deviation of signal intensities 

total: 109 words/record 


Fig. C-4 Format of plot- file data record written by PROFILE 


1 word 

55 words 

36 words 

36 words 

28 X 36 = 1008 words 

MONITOR 

BLOCK 

NUMBER 

WAVELENGTH 

INTENSITY 


MONITOR 1 word, R 

BLOCK 54 words, integer 

NUMBER 36 words, R 

WAVELENGTH 36 words, R 

INTENSITY 1008 words, R 

total; 1135 words/record 


Intensity of monitor signal for this profile 

(currently not used) 

number of shots at this wavelength 

Wavelengths 

INTENSITY (I,J) is the signal for the Jth 
shot at wavelength WAVELENGTH(I) 


Fig C-5 Format of fit-file data record written by PROFILE. 
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A Note on Program Documentation 

A code IS used to describe the parameters of some subroutines. 

Fo.r example, JLn IIM, 

INT R,1 Given intensity , 

the R indicates INT is real (single precision floating point) and the 
I means it's used only for input (i.e., the subroutine doesn't change 
Its value). Possible parameter modes are* 

F single precision floating point 
DP double precision floating point 
I integer 

S statement number, for alternate return 
L logical 
C complex, 
and possible uses are: 

I input only (unchanged) 

0 output only (changed, contains useful information) 
10 input and output 

W work area (changed, not meaningful on return). 
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Programs Listed 

REVERT 

TREAD*, OPT* 

PARAM 

BSORT 

STORES, START, SADD, SDROP, HADD, HDROP, ADDTO, 
FINDTO, EPUSH, EPOP 
PROFILE 

TIM, INTENS, LOOKUP, YESNO* 

VPLOT 

THEORY 

DBANK, GROUP, FETCHS, FUNCT3, FUNCT2 , TPLOT, 
AXISN**, NEWS, NEWT, NEWU, SIGMA, SYMSLV**, VALUE** 


Programs in UNIVAC Assembly Language. 
These programs may be of general interest 
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KEvtBT 


^Ob373JIM*AOK^SPAC£t(l) .HtVEHT 


1 

2 
3 
H 
5 
o 
7 
a 

9 

to 

11 

12 

li 

l<f 

IS 

lb 

17 

18 

19 

20 
21 
22 
23 
29 
25 
2o 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 
36 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


NA -IE... 

RFVERT 
PURPOSE.. . 

TO ACCEPT A TAPE PRODUCED DY THE DIGITAL DATA ACQUISITION 
system A'jo pro uct A file acceptable to PROORAS *PARA«*. 


usage... 


UXQT RFVERT 

cdata caros> 


{NOT 0, REVERT) 


options; 

•L« PRINTS information FROM DATA CARDS AND FIXED DATA 
FROM TAPE RECORDS (THUMRWHEEI. SWITCHES) 

*N* IGNORES Improper scale or interval from tape 
RECORD header... NO MESSAGES 

•R* OMIT initial rewind (DEFAULT: REWIND TAPE BEFORE 
READING) 

INPUT... 

data tape with name mhtapc*. produced by the digital data 
acquisition system 

data from thum.jwhc.el switches is interpreted a5 Follows* 
digits 1-5... scale tV/DiV) FOR CHANNELS 1-5 
DIGIT 6 ... sample interval (HICPOSEC) 

setting 0 123456769 

MEANING Illegal .005 .01 .02 .05 .1 .2 .5 1 . 2 

.005 >LE'. SCALE .LE. 2. 

.05 .LE. INTERVAL ,l.£. ?. 

interval specified on CARO 3 SUPERCEDES DIgIT 6. 

A WARNI-IG IS printed IF INTERVAL ISN’T .1 HCROSEC. 

the following data cards; 

CAROS 1 and 2... 

(72A1/72A1) alphanumeric FILE HEADER IMAGES 


CARD 3 » . . 


(15) ilUIoER OF CHANNELS BEING USEo 
(Fb.O) SAMPLE interval IN MICROSECONDS 
(Ib) 


(15) 


CARDS 4-N 

(15) 


monitor signal channel (if black* the program 
USES the FIRST CHANNEL WITH BLANKS IN THE 
wavelength specification COLUMNS OF CARO 4.) 
starting sample numbfr (if blank* The program 
reals the FIRST IZ RECORDS AND USES THE 
AVERAGE OF THE MAXIMA OF THE MONITOR. I 
(SORTED BY SHOT tti INCREASING) 

FIRST SHOT HUMBER OF A GpOUP OF SHOTS WITH 
THIS SET OF wavelength SETTINGS. 

(5F5.1) 'WAVELENGTHS IN ANGSTROMS FOR EACH CHANNEL 
BLANK FIELD INDICATES A MONITOR CHANNEL 

CARO N+1... 

•HEOF ♦ IN FIRST FIVE COLUMNS 
CARDS H*2 - H 

(IS) SHOT number WITH INCORRECT SCALE « 

(5F5.4) NEW SCALE #»S (BLANKS FqR CORRECT 0((ES) 

CARD H+1... 


0? Poor quality: 
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htvf’T 


30 

ti9 

60 

61 

62 

63 

09 

Ob 

66 

67 

oS 

69 

70 

71 

72 

73 
/4 
7b 

76 

77 
7d 
79 
60 
01 
02 
63 
04 
dS 
66 
o7 

ao 

d9 

90 

91 

92 

93 

94 

95 

96 

97 
9d 
99 

lUO 

lUl 

102 

103 

104 

105 
lUb 
107 
lua 

109 

110 
lU 
112 
113 


C ‘OtOf • It< FIP5T FI.'E. COUU^'IS 

c 

c ouTf>ur... 

C IN FlLf lOi A FORTRAN RANDOM-ACCESS FILE ACCEPTABLE 

C TO PROSKAH 'PARAM* 

C 

C OATA RECORD FORMAT; LAbEu(2) *SCAUE»‘ISHOT» Tl tTl ,T2» T2 ... T8»Y8 

C total; 20 hORnS 

C 

C SUOPhOSRAMS PeOiJIREO... 

c 

c begin 

c opt 

C TREAD 

c 

parameter PQINTSsA 
PARA lETEK MPTS=SO 
PARAMlTER MBAD=15 
PARA-tETEK FCHANS4 

Parameter HRcn=iono 

PAHAMC.TER N700DS=5fPOINTS*POINTS 
INTEGER H£ADER»CARD 
logical rar.ieo.sgu/cu 
LOGICAL OPT.C'IET.LO 40 

OlMDJblOll SCAlLIV) , ^AVFL(f.PTS^•'CHA^■> .k£COPO(N, ORDS) • 

- IREC(HrfOJDS) .KSHOT{kI“TSI .WIMCHANI .LS«0T{^8AD) t 

- GOOPfMaAD.NCUAtn 

data bCALE/.OnSf .01. .02. .05. •!> «2. .5.1. .2./ 
data f.LRR/O/ 

EOUI VALtfJCE ( record 1 1 1 . 1 RFC ( U 1 < ( HEADER 1 1 ) . RECORD In) 

OATA JPTS. nFILF, ENOhEC . NRO« CARO. WAPfJEO. SOlVEt 

- / It io.'oeof 9. 5i. false. ..false./ 

DATA ; 4CHAr'. /MRTS. M;Kcu. mvoag 

- / lACHAN. Mprs. iRCo. mrAD/ 

COMMO, . NSH0T.M0''tl2>.ISIv[5.64) 

CALL 6EGIM( 'REVERT 1.23 i.') 

<JUieT=OPT(»n>) 

LONG=OPTl»L') 

IFt.NOT.OPTt *pM )CALL Re*II 0 
OCFIRC FILE ..FlLC(^nCJ»ri.iOPDS.Ur.I iCC) 

NREC=1 

C 

c Transfer header inforaatioh 

DO 10 1=1.2 

READ {CARD. aiP , END=1 06 ) header 
aio FORMAT U2A6) 

IFILONGIPRINT 812.HEA0ER 

ai2 FORMAT! IX. 12A6) 

10 wHITfcCNFILE'flBEOHEAOEH 
C 

C 2E ALWAYS C.NTER THE MAXIMUM NUMBER OF DATA POINTS 

IREC(5)=P0IHTS 
C 

c read number of channels in use. sample interval* 

C MONITOR CHAt 4EL «i AMD STARTING CHANNEL « 

READ{CARD.a20.END=X06)NCHAN.SAHPLE.POHlT.KTIME 
820 FORHATU5.F5.0.2IS) 


ORIGINAL PAGE IS 
m POOI^ gUMOT 
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>»LVt.KT 


114 

115 
Ud 
117 



118 


12 

119 


15 

120 


18 

121 

122 


S 15 

123 

C 


124 

12 b 

c 


126 

127 

128 


825 

129 

130 


830 

131 


20 

132 


832 

133 


23 

134 

135 


835 

13 b 

137 

i 3 ri 

139 


25 

140 


836 

141 

142 


28 

143 

144 

145 


. 30 

146 

c 


1^*7 

c 


148 

149 

150 

c 


Ibl 


35 

152 

c 


133 

c 


154 

15 b 

c 


156 

157 


901 

158 


40 

159 

c 


IGO 

c 


161 

c 


162 

lo 3 

164 

165 

166 
167 

c 


168 ' 

169 

170 




IF(.KMA'l.<;T.'irHAIi .OH. m,hAN.L‘='.OICO TO 1.38 
IFtSAf PLE.EO.0.)0O TO 12 

ir <’>*' HLt.LT..O') .OR. SAhPLE.rr.a.lGO TO 110 
SGlVtic.TRlJE, 

IFlFO.Iir.LT.O ,OK. MO.llT.ST.MCHANlGO TO 112 
IFtKTlMt.UT.O .OH. KTir!E.6T.5'»160 TO ll4 
IF « LONG! PRINT 8l6.HCHA»,iA IPLF tHOMlT.KTlME 

F0RKAT(lX//lX,I2f • CHM.iltLS U':EO. ' *F6.?. • USEC SAMPLE INTERVALV 
• CHAHME1-*|I2.* IS MONITOR, STARTING CHANNEL IS».I4> 

HEAR wavelengths FOP EACH SHOT NUMBER 
IF<L0N6.AnD. .NOT.OPTCYMlPHt'IT 125 
FOHMATdX//* SHOT B WAVELENGTHS’/) 

00 23 NPTS=1*:PT5 

READICARO,830,ENO=25)KSHOT(MPTS),(WAVEL(NPTS,M),H=1,NCHAN) 

FORM AT (15, 5 AS) 

DO 20 P=1,NCHAN 
0EC0DE<tt32,PAvEHNPTS*M) )W(M) 

FORMAT <F5.1) 

IFlLOilG.AND. .NOr.OPTC »Y’) )PHINT a35,KSH0T IMPTS) , (V,(H) ,H=1 ,nCHAN) 
F0RHATtlX,I5»2X,5F9.2) 

GO TO U6 
i<SHoT(NPTS)=lnono 
HPTS=npTS-l 
DO 28 rlUAU=l'MDAD 

HEAOlCAKt},a36,ENV=30)LSHUTlNnAU) ,(6Q0D(NBAC*H)>M=1,NCHAN) 
F0«HAT(15.5F5.4) 

CONTII.UE 
GO TO IIT 
i.SHOT{N3AD)=ln000 
N6AD=N3An-l 
IFtMOMT.NE.OjGO TO 40 

WE HEREIJ’T TOLD THE MONITOR CHANNEL B... 

FIGURE IT OUT (PLANKS IN THE WAVELENGTH FIELD) 

00 3b MOMITsl.NCHAM 

IF( <AVELtl,MO lIT).EO. • M60T0 40 

CONTIt-UE 

THERE'S STILL '|0 MONITOR CHANNEL SPECIFIED... 

USE channel 1 

RkItiT 901 

Format (• no monitor specified. ••using chan. i») 
monit=i 

IFtKTIHE.NE.niGO TO 58 

WE WEREN’T TOLD WHICH TIME TO START WITH... 

FIND average TIME FOR PFAK OF MONITOR SIGNAL 
AMONG FIRST 12 RECORDS 

DO 47 L=I»12 

call TREADIHSHOT.KEOF) 

IF(KEOF.NE.O)GO TO 113 
HAX=IS1G(W0,41T.1) 

1=1 

DO 45 K=2>64 

IFtISIG»MO.IIT.K).LE.MAX)G0 TO 45 
I=K 
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««•*•* ]{ElVt;.RT 


17X 

Hi. 

176 
17<* 
175 
17fa 

177 

178 

179 
IdO 
101 
182 
153 
184 
18b 
IBb 
157 
155 
159 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
2U1 
202 

203 

204 
206 
206 

207 

208 

209 

210 
au 
212 

213 

214 

215 

216 

217 

218 

219 

220 
221 
222 

223 

224 

225 

226 
227 


MAX=t5IG(5Jr4lT»K) 

4b CONfl! Ot 
47 KTMt=<Trhb«-l 

KTl HE=!UrjO ( S7,KT19E;/12) 

CAUL liEWINO 
IFtOHTt'2'))STOP 
53 INR£(.=0 
MSMOT=a 

IFtLOfJGJPRiNT 835 

838 F0HMAT41X//' KECOHD SHyT » FIXED DATA*/) 

00 82 L=1>MRC0 
C 

c read a record 

CAUL TRtADlNSHOT*KEOF) 

IFIKE0F.tIE.0160 TO 05 
HSHOT=HbHOT+l 

IFILONS) PRINT B40»HSHOTih,ShOTM4Uw 
B40 F0RHAT(lX.I4.I9.BX.5I2>2Xr J2>2X.6I2) 

II,REL=JNHEC+l 

C 

c FMD TH*^ KAV£LE’'STH INFO FOR THIS SHOT N!lMt,ER 

IF«N5ttOT.LT.XSHOT( JPTS) )JpT9=t 
1ST0P=NPTS-JP*'S+1 
OO 62 I=1.IST0P 

IFtNbliOl.UT.KSHOn JPIb+1) leo to 61 

62 JPTS=JPT5vl 

63 CONTltJUE 

00 82 H=1.NCHAN 

c 

C load DATA PROP O'lE CMAtl'iEL INTO A RECORD 

C 

c DISCARD Data IF OESIRED 

IF(V.4VEL<JPTS»M)-Ea. '9999 ► .OR.WAvELI JPTS»H> . EO- • 99999 * JOO TO 82 

C 

c PICK UP SCALE » l-BO-t PEC"RD HEADER 

15=NUK(M) 

IFIOPTt'n* ni5=riUH(H-l) 

IFUS.6T.0 .AMD. IS.LT.aiGO TO 6S 
IF(.‘«JT.OUlET)PPlf T 902.1HREC.MSHDT,P.*S 
KERR=r ERR+1 

902 FORMmTC TAPE RECORD ' 1 14, • , SHOT' . IS. • *CHAHHEL' » 12. 

- *... SCALE a OUT OF rtA.IGEI'.IS/ 

- • USING ,005 V/DIV*) 

1S=1 

65 K=KTIME 
C 

c SUBSTITUTE CORRECT AMPLIFICATION IF NEEDED 

AHPLFY=SCALECrS) 

66 IF(LbflOTtJ'JAD).EQ.1000qiG0 TO 63 
IFdJbhOT.LT.LSHortJBADDviBADsl 
IFtNbHOT.LT.L5HOT{JBAi>+H)GO TO 67 
UBAD=JBAU+1 

60 TO 66 

67 IFtMSHOT.EQ.LSHOrUBAJ) .AMO. G000( JBAO.M) ,NE. 0.0> 

- AMPLFY=GOODt JOAO.M) 

68 IFtSGIVc.N)60 TO 71 


ORKrCfAB PAGE K 
PB tOOR liUAlffiS 
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I'E-Vt-PT 


22« 
iiM 
2 JO 
2^1 
2J2 

C 


PICK UP sample interval F(?0 I PEfiARL HEADER 

lt=.JU i{o) 

IFtIl.or.3 .AMD. IT.LE.9) oO TO 7fl 
NE«fl=I.E6l *1 

IF< .AOT.i.JlCTlPRlMT 'lOJilVlIEC.IlSHOTtlT 

2J3 

2u4 

23b 

236 


903 

rOHAATt' TAPE RECnnO'.l4f'. SHOT*. IS. 

•...sample interval b out of range:»*i5/ 
- • USING .1 MICROSCC) 

1T=2 

237 

236 

C 

70 

SAMPLt=5CALEtIT) 

2j9 

C 


GIVE ONE WARNING IF INTERVAL ISN'T .1 US 

240 

241 

242 

243 


71 

IFISAPPLE.E0..1 .OR, VARnEOIGO TO 72 

HEHR=NEKR*1 

WAHNEj=.TkUE. 

PRINT 9U4.lNRrC*MSH0T,SA.^PLE 

244 


9u4 

FORMAT! • TAPE RECORD* »Uf' SHOl • , J5, < , . .SApPLr INTERVAL I.OW'i 

245 

246 

c 


- F5.2.* MICROSeCONUS'l 

247 

c 


TRANSFER Data POINTS 

246 

249 

250 

251 

c 

72 

TIMt=K«5AvPLE 

00 75 hK=7..w'('HPS,2 
RECOHO!KK-l)=Tr'£ 

252 

c 


WE FAKE A (.ORRFCIIOM OF A FACTOF OF 10 6ECAUSC 

233 

254 

255 

c 


, TMtPt ape Tiio EX IRA A"PLIF1EPS IN THt DATA SYSTEM 

REC0rtC(KK) = ,l*ISIG(Mr<O*ANPLFY 
K=K41 

2jt> 

257 

c 

75 

TME=TIME+SA SOLE 

250 

259 

c 


SET UP output REC0R9 HEADER 
ENCOUt !650,1RECI*AV£L[JPTS,M) fNSHoT 

2o0 

261 

262 

2o3 


350 

FORMAT (AS, 1 3) 

REC0hD(3)=AHPLFY 

IREC!4)=;JSH0T 

WB I TE ( NF ILE • r IP EC ) PECO,- D 

264 

265 

266 

c 

62 

CONTl.NUt 
GO TO 124 

267 

c 


END OF FILC...OUITTIMS TIME 

2u6 


85 

PRINT SbO.PSHOTrNOEC 

269 

270 


BbO 

F0RrlAT(lX,I6,» SHOTS PROCESSFO'/IX, I6» ’ RECORDS WRITTEN') 
IF(i‘Ci P.GT.OIrRIi.T 90S>f,_Rrt 

271 

272 

273 

c 

905 

FORMATdX.Io,' ERRORS CR aAo.JIiJGS') 
GO TO 199 

274 

c 


complain 

275 


1C6 

PRINT 906 

276 

277 


906 

Format (• header cards arc missing*) 

GO TO 199 

278 


108 

print 91.6,1 chan, whan ' 

279 

260 


908 

FORHmTC NU'IHEH of CH.\(,NfcLS,',I5,'» IS OUT OF RANGE 1 TO* >12) 
GO TO 199 

261 


110 

PRINT oiO,SA*'PLE 

262 


910 

FORMAT!* sample INTERVAL 0F*,Gl5.5,* IS BAD*/ 

263 

264 


- • USING TAPE RECOPD HEADERS*) 

GO TO 12 
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'<LVt«T 


20b 

112 

PhlNt 912.M0N1T 

20b 

912 

FORiWrC* 0IVC9 f'OMlTOR CHA.J'JFL IS "ATI 

207 


10NIT=0 

2o0 


GO TO 15 

2d9 


PRItlT 914,KTI‘'E 

2sa 

9m 

F0H9hT(‘ 01v£ 9 STIRTMO tHANlEL IS oAD') 

291 


KTli-(t=0 

292 


GO TO 10 

293 

116 

print 916.FMPTS 

299 

916 

FORXATl* TtiERF ARE MORfc, THAN*tl3»* WAVELENGTH CAHJS'J 

295 


GO TO 199 

296 

117 

PRINT 917^^.V3AD 

297 

917 

FOR.VIATC THERf=’ ARE MORE THAf'*»l3»* CORRECTION CAROb’l 

290 


CO TO 199 

299 

118 

PRINT 918 

3U0 

918 

Fg 9'’AT(' FF.REP Tinv 12 S16.IALS ARE PfiRSENTM 

301 


GO TO 199 

302 

129 

PRINT 924..4MPC0 

303 

929 

FORMAT (• MORE THAH»»Ijt* TAOE I'ECOROS. . .QUITTING' > 

3U<f 

C 


3U5 

c 

CLOSE OUTPUT file 

30d 

199 

WR ITl ( NF IlE • nr EC ) ENDREC 

307 


STOP 

306 


•lno 


tJWORK.LJCT 


DBIGINAIJ PAGE IS 
OP POOR 



****** revert 


66 


(Sample data) ****** 

Lo.ng printout .is wanted. 

F^R':‘'.SURE='0'j' =.L>0T=/,7V» Rt FLCC FOR =.B 
6, AL II 16VC ... FILTdR RUf! 101 

3 channels were used. 

First shot number at these wavelengths 

Channel 1 set at 3280.0 % 

Channel 2 set at 4688.0 8 

Discard data for channel 2, shots 84-92. 


Channel 3 is for monitor signal. 


Shot 76 had wrong amplifier setting recorded 

Channel 1 amplification was .005 V/div 
Channel 2 amplification was correctly 
recorded 

Some signals were attenuated 
Attenuation factor was .151 . 

Signals for 4670 A through 4696 A were 
attenuated 
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'’EVcKT [TIEA'-i» 


20b373JlN*^Or, 


A*3!‘At.CiU) .TtliAJ 


1 

i 

3 

5 

6 
7 
6 
9 

i.0 

U 

12 

13 

14 

15 
lt> 
17 

la 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
29 
JO 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 
46 

49 

50 

51 

52 

53 

54 

55 
5b 


rU'iE... 

DA TAT APE 
PU4P0SE... 

TO READ TilP HIOH SPEED DATA ACOUISITIOm SYSTEM TAkE. 

CALLING SEQUENCE... 

call tread tBUFKEHiEOF) 

BUFFER 333 WORD DATA INFUT DUFFER, IMTrQER 
BUFFER<1) IS THE EXPERIMENT COUNT 
B'JFFER(2-13) ARE HISC FIXED DATA 

BUFFER( 14-333) AKE DATA POINTS, DIVENSIONEU (5,64) 
EOF E*10 OF File. -.ARfLR, SET ‘I0N-2EP0 IF £0F 

WAS DETECTED 4HEN THIS BEAD WAS ATTEhPTEOi 
ZERO otherwise. 

ONE TAPE RECORD IS REAP ANQ UNPACKED INTO THE USER'S flUFFtR 
call PEWl'lO 

THE data tape is rEitOUNU. 


input... 

tape record fop*'at is: 

6 BYTES (3cn. 2/BYTE) FIXED DATA 

3 BYTES (BCD) EXptRlHENl COUNT 

320 BYTES (dlflAHY) DATA, GROUPED BY TIME 


s(l) 

AXR*. 




• 

% 

MAIM ENTRY POINT 



o 

TREAD* 

LA 

AO *EnO. 

ARE WE AT THE END OF THE FILE? 


JHZ 

AO, 3, XII. 

IF SO RETURN AT ONCr 


ANOTHER 

LA,U 

A0,PKTT. 




ER 

lows. 

FETCH A BUFFER FULL 



LA 

A0»p><TT+3. 

PICK UP STATUS Ann CHECK. 



TE 

AO, STATUS. 

SHOULD HAVE AN AFC OF 5 



J 

EOF. 



OKAY 

SX 

XltSAVEX. 




LXM 

XlfOfXll. 

XI WILL POINT TO the NEXT 

WORD IN THE 


LXI.U 

Xlfl. 

USER’S PUFFEP ALL THROUGH 

THE program 

• 

• 

PICK UP EXPT COUNT 




LA 

ao.inbuf-h. . 




AND.U 

A0.0I7. 

WIPE OUT EXTRANEOUS STUFF 



LA 

A2>A1. 




SSL 

AO, 8. 




ANOrU 

A0,Ql7. 




MSI,U 

Al,10. 




AA 

A2,A1. 




SSL 

AO, 8. 




[DRIGEJAE PAGE IS 
PS POOI^ ^UAUT^ 
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REVLfiT ITREVO) 


b7 


ARDfU 

A0»ni7, 




M5I»U 

AUlOO. 


M 


AA 

A2*A1. 

EXPT COUNT IS IN A2 

bO 


SA 

A2»0f*Xt. 

STORE EXPT COUNT IN DUFFER 

«»1 

* 




o2 

• 

PICK UP MISC DATA 


bi 

• 




64 


DL 

A0»I«3UF. 

SIX BYTES OF TbO CUARACTERS/BTTE 

oS 


LRrU 

K2»5. 


6b 


DSC 

A0>4. 

BACK UP ONE ni3IT 

67 

LOOPl 

LOSC 

A0>12. 

finVF FORWARD 3 niBiTS. PUT IN Al(3-0> 

68 


ANDrU 

Alf0l7. 

MOVF LAST dibit TO A? (3-0) 

o9 


SA 

A2fO»*Xl. 

STORE DIGIT 

70 


OSC 

AQf4. 

MOVE SACK 1 digit 

7i 


AND tU 

A1»017. 

PUT ONE digit into A2(3-0) 

72 


SA 

A2»0i«Xl, 

STORE BCD DIGIT 

73 


JGC 

R2*L''0P1. 

OECPEVE)JT COUNTER 

74 

• 




75 

• 

6ET REAL DATA 


7b 

• 




77 


SX 

X2.SAVEX+1, 

X2 POINTS TP NEXT PAIR 

73 


LX 

X2»(?»INHUF+2) 

OF WORDS IN INPUT puffer 

79 


LR»U 

R2»34. 

HOVE 35 pairs of words (315 BYTES) 

dl) 

100P3 

OL 

AO>0»«X2. 

PICK UP TYO packed WORDS 

8i 


L«»U 

Rl>8. 

NINF WRCS AT A SHOT 

toZ 

LD0P2 

LOSC 

A0*8. 


a3 


AND'U 

A1»0377. 

TRANSFER ONE BYTE To A2(7-0) 

<»4 


SA 

A2.0r*Xl, 


bs 


OGO 

R1.L00P2. 

shift until "oinLE ford is finished 

Ob 


JGD 

H2rLOOP3. 

finish 35 WORD PAIRS 

37 


OL 

AOfOrX2. 

PICK UP LAST TWO words 

68 


LrtfU 

Kl>4. 

process ONLY S BYTES 

89 

LO0P4 

LOSC 

AO. a. 


90 


Ai40>U 

Al.0^77. 

transfer OHS' BYTE To A2(7-0) 

91 


SA 

A2.0.*X1. 


92 


JCD 

R1.L00P4, 

finish the last 5 dytes 

93 


52 

•l.Xll. 

CLEAR EOF Flag 

94 

RETORU 

LX 

Xl.SAVEX. 


95 


LX 

X2.SAVEX+1. 


96 


J 

3*X11. 


97 

• 




98 

• 

END OF FILE OR TAPE 

ERROR 

99 

• 




100 

EOF 

LArSl 

A0»PkTT+3. 

TEST FIRST FOR EOF 

101 


TE»U 

AO. 01. 

01 r^EAtJS END OF FILE. 

102 


O 

BO. 


103 


SA 

AO.E^'D. 

hotf that EOF reached 

104 


SA 

AO.*lfXll, 

TELL calling PROGRAM ABOUT EOF. 

lOS 


J 

RETURN. 


106 

BO 

LAfH2 

AO.PKTT*!. 

, TEST FOR WRONG RECORD LENGTH 

107 


TNEfU 

AO . 74 . 


108 


J 

Bl. 


109 


SLJ 

PRINT* 

length wrong... PRINT STATUS 

110 


ER 

ABORTS. 

AND QUIT. 

111 

B1 

LAfSl 

A0.PKTT*3. 

GEt I/O COMPLETION STATUS CODE AGAIN. 

112 


TtlE»U 

AO. 00. 

00 MEANS NORMAL COf^pLETION, (MUST BE 

113 


0 

OKAY. 

COPIED TAPE, SINCE FRAME CT IS NORMAL) 
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‘?EVlRT (tkea[i> 


114 


TtiU 

A0>Q4. 

04 MEAN" ABNORMAL Fftl.it CODf.T. 

Uj 


J 

82 • 


116 


61.J 

PRINT. 

AFC MOT 5. SO wr PRINT STATUS AND 

li7 


J 

another. 

GET NEXT RECORD. 

Hi) 


SLJ 

PRINT. 

S0« P OTIIFB t/O fRROh. . .phINT STATUS 

114 


ER 

AUORTS. 

Ann ouiT, 

IcO 

• 




121 

t 

SUBROUTINE FOR PRINTING OUT THE STATUS WORn 

122 





123 

TRIUT 

-f 

$-$ . 


124 


LA 

A0rP<TT+3. 

GET the status WORD 

125 


LH.U 

R2»l. 


126 

£2 

SA 

A2.I.inuF+l. 

(USEFUL only on 2ND ITERATION OF LOOP) 

127 


LR»U 

Rl»5. 

t2R 

£1 

ANLfU 

A0r07. 

AOVF ONE octal digit TO AU2-0) 

124 


AA>U 

AltUSO. 

CONVERT TO FIELDATA DIGIT 

130 


DSC 

A1.6, 

MOVE INTO A2( 35-30) 

131 


SSL 

AC>3. 

MOVE next octal digit to AO<2-0) 

132 


JGD 

HlrEt. 

133 


UuO 

R2<E2. 


134 


SA 

A2. InUUF. 


135 


LA 

AOf (0104.1NBUF 

-2). 

13b 


El< 

PRINTS. 


lo7 


J 

.PRINT. 


l30 

• 




139 

• 

REwiNo entry point 


140 

• 




141 

HZ/iim* 

LrU 

AO.HPD. 


142 


Eti 

lows. 

REWIND THE TAPE 

143 


LAfSi 

A0.R-'D+3. 

CHECK FOR BAD STATUS 

144 


TNZ 

AO. 


146 


J 

RE2. 


146 


SLJ 

print. 

STATUS bad... PRINT THE STATUS 

147 


EH 

ABORTS. 

AND OUIT 

148 

RE2 

S2 

END. 

'JOTE WE AREN’T AT AN END OF FILE NOW 

149 


J 

l.Xl] . 


130 

• 




151 

• 

storage area 


152 

» 




1S3 

»(0) >5AVEX K£S z. 


154 

PKTT 

I SOT 

'INTAPE’ »RS 74. 

INBUF. 

155 


+ 

• 1/0 STATUS ’ 

• 

156 

INOUP 

RES 

74, 


157 

STATUS 

t0420[)s00nil2. 


158 

EtJO 

t 

0 . 

SET NOmzEPO “HEiI EOF FOUND 

159 

RWD 

I SOT 

•INTAPE’.REWS. 


IbO 


£.J0 




U.EJCT 


ORIGINAL PAO^ a 
OP POOR gUAlilSJ 
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HEVcWr 


jl^»i,ORKSPAtEi 1 1 > .OHT 


1 

o 

• 

)Cn ^C. * 9 * 



< 

3 

'» 

• 

« 

• 

OPT 



5 

• 

PURPOSE. . , 



& 

7 

• 

■ 

TO OdTftIM FOR THE USER PROGRAM THE OPTIONS SPECIFIED OH 

8 

• 

the 

EXECUTIMG STATE'!EHT. 

9 

10 

• 

% 

CALLIHG SEQUENCE... 


11 

• 




12 

* 

logical switch. TEST rOPT 

13 

• 

• 



!<* 

a 

* 



lb 

• 

SWITCHzOPTC'Un 


16 

• 

TEST=OPTt 'T* ) 


17 

• 

• 


-r 

IB 

• 

SWITCH WILL have THE 

VALVE .7RUF. IF ThF OPTION WAS 

19 

• 

specified C"J The «.XCT ok QfILE.prOGRAw CAhOf A!0 •FALS'^. 

20 

• 

OTHERWISE. TEST WILL 

SI"1LARLT indicate the presence of 

21 

» 

' THE 

*t* option. 


22 

SCO) 

LIT. 



23 

SCI) 

axpa. 



24 

OPT* 

T2 

HAI/E 

go get options if we HAVEN'T ALREADY 

25 


J 

GET 


26 

GO 

LA 

AI,*Q.X11 . 


27 


SSL 

A1.30 . 


28 


A VA.U 

A1.06 . 


29 


LA 

AU.WOrO . 


30 


LSSL 

AO.lO.Al . 

SHIFT CORRECT GIT TO HIGH-ORDER BIT 

31 


SSL 

A0.3S . 

SHIFT TO LOW-ORDER CHESULi IS 1 OP 0) 

32 


J ' 

2.XI1 . 

RETOR'! 

33 

get 

SX 

XU.WORD . 


34 


£K 

OPTS 

GFT OPTION WORD 

3b 


LX 

Xll.rfOPD . 


36 


SA 

AO. WORD . 


37 


S3 

HAVE 

NOIE OPTION WORD IS PRESENT 

38 


J 

GO 


39 

SCO) . 




40 

RETUR.) Rilb 

1 


4l 

HAVE 

4 

1 


42 

w02O 

HtS 

I 


43 


END 

* 



W.EJCT 
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(•AHA^' *«*•«« 


2OS>i73JlM*<«0hKil'ACCl < 1 } .PARA« 

1 C 

2 C NA^E>.. 

3 C PARAM 

4 C 


6 

7 

a 

9 

10 

11 

12 

IS 

14 

15 

u 

17 

la 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 
90 

51 

52 

53 

54 

55 

56 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


PUHPObE... 

ro CREATE RECOKUS WITH bTRICTLY DECREASING RONlTOh SIGNAUb 


usage... 

U.PAAAH OR OlOT .PARA4 
OPTIONS! 

A loXOT 0 (LY) amplify SOMF SIGNALS WHICH WERE 

ATTENUATED WITH A MEUTTIAL OEnSITY FILTER 

<WL0W»WHI6H> {FREE) SIGNALS WITH WAVELENGTHS IN ThE 

RAMGE (ViLwWi JHIGU) WILL IE A»<ML1FIE0 HT 1/.294 


IHPUI... 

ACCEPTS FIUS CREATED Rl PROGRAMS *RCOVER* OR 'REVERT*. 

RECORDS I A 2! 

(72A1/72A1I file HEADER 

RECORDS 3-H! 

WORDS 1.2; {F5.KI3) WAVFLE'IGTH.SHOT » 

WORD 3: IR) SCALE. VOLTS/DIV 
WORD 4: tIJ SHOT ft 
WORD S: I II ft POINTS IN TplS OECORO 
HORPS 0-13: (R) tI«CS (USEC) 

WORDS 14-21J tR) signals (V) 

RECORD N»l: 

*BEOF • IN FlfiST (ORD 


OUTPUT... 

file ADCEPTABLE to program 'BSORT* 

RECORDS 1.2: 
t72Al/72Al) file HEADER 
RECORDS 3-tI. 

WORDS 1.2; (FS.1.I3) WAVELE’IgTH.SHOT « 

WORD 3: IR) SCALE (V/DI'M 
WORD 4: (1) SHOT U 

WORDS 5-20! (R> pairs OP TI^E (USEO, SIGNAL (VI 
RECORD N+l: 

'heof * In first «ohd 

SUBPROGRAMS REOUIREO... 

BEGIN. OPT 

PARAMETER HRCDPIOOO 
integer 0UTF1L.0UTREC.TEST.END1 
Integer head ( 201 .high. vallev. peak 
logical opt 

DIMENSION A(8) .B(8) .FACTOR(B) .WL0W(8> .WKIGK(S) 

EOUl valence t a 1 1 ) .HEAD ( 51) . ( B tl J .HEAD (13)1 
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param 


&7 



DATA ENDl.INULE.OUTFIL. I'lREC .OUTKEC 



- 10. 15. 0. 0/ 

t>9 



OXMENSIOH X{A) »Y(P) .<(8) 

60 



E6DIV..Le.‘'Ct (XUI.MID 

6i 



LOGIC. L S«ME 

oZ 



CALL EEGINI’PARAH 2.U C) 

63 



AATT£l.=l 

69 



lFt.UOT.OPTt'A')»GO To 5 

6b 



PSlflT ooo 

66 


ooo 

FOHMAK* EtJTEn ATTFN AND WAVELENGTH RANGE*) 

67 



DO 3 KATT£.«1.8 

68 



read eo5.END=5.FACTOH(KATT£Hl . aLO <tKATTEN) .WH igH(KATTEN) 

69 


805 

FORMAT ( > 

70 



WLOW tkATTEM =^LOW ( K ATTEN) 05 

71 



WHlGHCKATTEN)=WHlGH<KATTLN)+-nt> 

72 



IF«FACT0R(KATTEN>.GT.1.)PRINT 901 

73 


901 

FORMAT (• >1?*) 

79 


3 

CCNTI,.Ut 

7S 



KATTEl4=8 

76 


5 

WLOWIKATTEN)5ROOO. 

77- 

, ^ 

WHIGH(kATTCN)=799«». 

78 ' 

* 


define file IflFILElHRCO.il. U.IN rec) 

79 



OEFIIJE file 0UTFIL(MRCD,20,U,0UTREC> 

80 

C 



81 

c 


TPAt’SFER FILE HFAOER 

82 



00 3 J=1.2 

83 



READtINrILE'J) (HEAD! I ). 1=1 . 12) 

89 


8 

WRlTElOUTFIL*J)(HFAD(n. 1=1.1?) 

85 

c 



86 

c 


REAr A RECuRO 

87 


10 

READ(IMFlL£*IflHEC) (HSAO( I) . 1=1 .9) ,H, (X( I) . Y( I) . I=1.N) 

88 

c 



o9 

c 


QUIT AT EHu OF FILE 

90 



IFtHEAOtl) .EO.EHDl )G0 TO 90 

91 



akplfv=i. 

92 

c 


REMOVE BAD RECORDS 9 CORRECT AMPLIFICATION 

93 



DEC0DtI830.HEAa)nAVE 

99 


830 

F0HMAT(F5.l) 

95 



DO 12 KK=X.KAtTE.) 

9b 



IF(WAVE.LT.WLO*(kK) )G0 Tu 12 

97 



IFtWAVE.LE.WHIGH(KK))GO TO 19 

98 


12 

CONTIt.UE 

99 



GO TO 15 

lUO 


19 

akplkt=factorckk) 

101 

c 



102 

c 


transfer record to output area 

103 


15 

N=MXi46<N.8> 

109 



DO 20 1=1. N 

105 



A(N*1-I)=xa) 

106 


20 

B t N+l-I ) = Y 1 1 ) / AMPlF Y 

1U7 

c 


• 

108 

c 


SORT POINTS ON TIME 

109 



IFIN.LE.DGO to 90 

110 



DO 28 J=M.2.-.1 

lU 



same=.true. 

112 



DO 26 1=2. J 

113 



IF(A(II.GE.A(I-1))GQ to 25 



73 


PflWAN 


114 



5=A(I> 

11& 



ACl)=A<t-lJ 

116 



A« r->ii=5 

U7 



s=dct 

no 



dtDsHi 1->1) 

IIP 



8{ 1-11=5 

li>0 



SA AES. FALSE. 

121 


25 

coutx^ue 

122 



IF < SAME! 60 TO 30 

123 


£d 

CONTI'JUt 

124 

e 

\ 


12S 

c 


IF THIS IS A MONITOR SlpNAL. ENSURE 1T*S 

120 

c 


MO.JOTONICAi.UY ofcreasing 

127 


30 

DECODE 1 54Q .HEAD) Tf ST 

120 


a'to 

FOKhATCAS) 

129 



IFJTEST.Eu.* »>G0 To 31 

130 



IQSN 

131 



toO TO 44 

132 


31 

DROPSO. 

133 



HIGH=l 

134 



2ALLEY=1 

135 



PCAK=1 

136 



IB30 

137 

c 



13d 

c 


find a LOC*tU yAYI'Vi 

139 


32 

Ih=X3*-l 

140 



IFIIO.QT.nJGO to •'A 

141 



IFtdU8).LE.D(M-APS{IB+l-H)) )G0 TO 32 

142 

c 



143 

c 


note location of local maximum 

144 



HIGH=I0 

145 

c 



140 

c 


FI.IO A LOCAL HlUirU.'. 

147 


34 


14S 



IF( IB.GT.NJGO TO 36 

149 



IFtimU).CE.PtM-ABStID+l-NH JGO TO 34 

160 

c 



151 

c 


IS THIS Drop bipger Than previous biggest? 

132 



rFfBtl)IGH)-B(IB).Le.DHCpjGO TO 32 

153 

c 



154 

c 


NOTE THIS IS PIGOEST DROP 

155 



DROP=L (HICH>-o(lO) 

150 



PEAHSI.IGII 

157 



VALUEYrlH 

150 



60 TO 32 

159 


36 

IFtathlGHl-BtlO-lJ.LE.DROPlGO TO 3fl 

160 



PEAKshIGH 

161 



VALUtY=I9-l 

162 

c 



103 

c 


HAVE Fourjo ’ biggest DR0P.,,GET RIO 

lo4 

c 


OF POINTS BEFORE PEAK 

lOS 


3S 

13=0 

166 


40 

I3=IU*1 

167 



IF<IB.GE.PEAK)G0 to 42 

168 



A(IB}=A(PEAK) 

169 



B(lB)=a(FEAKl 

170 



GO TO 40 



**•*«« 


pAa^'i 


***•*« 


in 

iiz 

1^3 

174 

l7b 

I/O 

177 

l7<i 

179 

lao 

lai 

lt)2 

103 

104 

las 

106 

107 

108 
109 
190 


c 

c GET RID OF POINTS APTFR VALUft 

40 IC=VALLfeY 

44 in=w*i 
lF<Iu.GT.O)Cn TO 45 
AnjJ=A(VALu:V) 

B ( IB) =U( VALUEV) 

GO TO 44 
C 

c WRITE A RECORD 

45 WH1 Tc( 0UTFIL*0UTREC)HCA0 
GO TO 10 

c 

c CLOSE THE FILE ANO EXIT 

90 WRITEIOUTFIL'OUTREOEiOl 
I=IHHEC-4 
PRIHT 690,1 

890 format (I5»» CURVES PROCESSED*) 

STOP 

END 


U.EJCT 
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'jbcur •«**»♦ 


2Ub37jJl»««aOhK<ot AL£S( U.n.i'll.T 


I 

? 

i 

It 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 
1» 
18 
17 
10 

19 

20 
21 
22 

23 

24 

25 
2E> 
27 
20 

29 

30 

31 

32 

33 

34 

35 
3t> 
37 
30 

39 

40 

41 

42 

43 

44 

45 
4b 
47 
40 
49 
90 

51 

52 

53 

54 

55 
Sb 


C 

C NAME... 

c asoST 

c 

C PUHP05E... 

C TO SORT THE RECORjS. FIRST BT WAVEtENGTHr THEN «JT SHOT «. 

C 

C USAGE..* 

C 

C M.fiSOHT Oft aXQT .OSORT (NO OPTIOfJSJ 

C 

C SUOPrtCOftARS CAELEn... 

c 

C 0E6lN.STAKT.SADD.5DROP.HAOOiHDROP,AORTO.FINDTO»EPUSH.EPOP 

C 

C METHOO... 

C on A FC^R'iAftO PASS» EACH RECORD IS REAII. IF IT RANKS HIGH 

C ENOUGH (HIGHER Th«N LEAST ONE Ct'RREfjTLy SAVED* IF THE 

C BUFFERS AHE FULL) IT IS SAVED Ann. IF THE BUFFERS ARE FULL* 

c the least one currejtlt Saved is *<rittfh in its place. «hlh 

C the end of the UflbORTED PORTION OF THE FILE IS REACHED* ALL 

C RECORDS REINS HElJ APE EXC (ANGFR *(ITH STORED RECORDS. THE 

C PROCESS IS REPEATlD, IN ALTERNATING DIRECTIONS* UNTIL THE 

C entire PIL®’ is S0HTS3. 

c 

INCLUDE ST0RES*U5T 
LOGICAL FOUND 
DATA Ei/*8E0F */ 

DATA NFlUt/15/ 

CALL FEGInI'BSOHT 2.01 O*) 

DEFINE file N=^ILE(100n*20.U*NR£C) 

call STAHT 

MEXT=2 

POUNO=.FAl.SE. 

LAST=0 

IRECS2 

INCKSl 

NPASi=0 

N»RITE=0 

NflEAJ=0 

KNOm.sO 

SURES.TRUC. 

C BEGIN NEB SBEEP 

10 CHAN5E=.FALSE. 

WPASS=HPA5S+t 

C CALCULATE HE^ PPCORC T(UMf»ER 

12 1RECsIREC«^INCR 

IF(I<<EC.EO.LAST)60 to 40 
IF((KNOBN.LE.O).AriD.(.(IOT.SURE))GO TO SO 
lF(HJASE)lb*lS*14 

14 IFtlSLC.EQ.FROHtHnASE) >80 TO 24 

C read neb ftkCOWJ 

15 CONTINUE 

REAOCNFILE* IRFO (R (I .BUFFER) tIsl.kOH> 

NREAU=NftEAD>l 
IF (FOUND >G0 TO 20 


OKIGIN-ab page J3 
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t SOnT 


i>7 

C 


STOP SitEP IP T‘US IS Lt’” OP FILE 


t>8 



IF(H(l»UUFPEO).fJL.El)»jO TO 20 


b9 



FOUfIt)=.THUt. 


bO 



LASrsIREC 


ol 



RECRii=IREC-3 





GO To 9(1 



C 


AOD THIS RlCORO to SORT STACK 




20 

CONTI. Ufc. 


bS 



FROM { PUFFER 1 SIREC 





Call 5rtOO(ilUFFCR) 


67 



OUFF£K=OCMAND H89 1 


66 



60 TO 12 


69 

c 


WE ARE HOLUINO A RECORD WE WANT TO 

1 WRITE here... 

7o 

c 


00 SO. THEt, ADD IT TO THE SORT STACK 

71 


29 

CONTINUE 


72 



wrtlTEUIFlLE'lREC) (RII.HpASFJ .tsi.LGH) 


73 



N»HITL=N<RITE+1 


79 



CALL 5AD0IHDR0P(SB9>) 


75 



IF<SJ>E>GO TO 12 


7b 



CALL CRUSH lSDR0P(?09n 


77 



bO TO 12 


70 

c 


THERE'S NOT enough ROOM To 00 AnT 

sorting LURING 

79 

c 


THIS SCrtN..,nTPASS ALL RFADING 


dO 


30 

IFCHoASE.LE.OlGO TO 32 


tli 



lREC-FP0.tliDA5EJ 


02 



WR ITfc. ( NF I LE ' IrEC ) ( R ( I . H0 ASE > # 1=1 f LGH ) 


83 



NWRITE=NWRITC:+1 


89 



CALL tPUSHtH0p0P(S891) 


as 



60 TO 30 


86 


32 

IREC=LAST 


07 

c 


CALCULATE NEW SWEEP LIMITS 


as 


90 

CONTIr.UE 


69 



LAST=riEXT 


90 



NEXT=lREC-rbCP*KNOWN 


91 



1F( (.NOT. CHANGE) .AND.IKriOJkJ.GT.O)] 60 TO SB 


92 



KNOWN=0 


93 

c 


RECORD 3ESTI lATIOHS OF RECORDS TO 

3E HELD 

99 



IsSSASE 


95 



NTOSO 


9b 


51 

IF(I)55.S5.52 


97 


52 

CUNTlUUb 


9b 



CALL ADDTO(l) 


99 



IsUPlI) 


100 



GO TO 51 


101 

c 


SET POINTERS 


102 


55 

CONTIf.ue 


103 



J=LAST-INCR*MTO 


109 



JT05NT0 


105 



CALL FINOTOIJ) 


106 



INCR=-INCR 


107 



SURE=.TRUE. 


108 



TTOP=STOP 


109 



STOP=t) 


110 



SBASe=0 


111 

c 


WRITE SORTED RECORDS IN PROPER ORDER 

112 


58 

IFITTOP)79.79f60 


113 


60 

IRECSIREC+INCR 






tiSOHT 




114 



If lliltC.E0.F9''>'1( TTOP) ir.O TO 64 

U5 


bl 

IFt(ru(JT0)-l9CC)*INCit)b2i63t6U 

lib 


62 

JT0=JT0-l 

U7 



GO TO ol 

118 

C 


FILE RECORD THAT WAS HERE 15 IH THE SORT 

119 

c 


STACK AiiD CAH RE OVERWRITTEN 

120 


63 

contd.ue 

121 



WRITUiJFILC'IoEC) (RIIiTTOP) . Irl.LCHl 

122 



NWR1TE=NWRITE+1 

123 


64 

IsDOJtMTTOF) 

124 



CAUL EPUSh(TTOP) 

125 



TTOPSI 

12b 



00 TO 58 

127 

c 


PILE RECORD THAT WAS HERE HOST RE HELD 

128 

c 


AND WRITTEN BACK LATER THIS SWEEP 

129 


68 

CONTI) U£ 

130 



READINFILE'IHEC) (P { I .I’UFFER) . t=l .LGH) 

131 



NREA0=NREA041 

132 



UftlTE(NFlLE'ItlECHRtI.TTOP).I=l.LeH) 

133 



NWRirtS'tWlTE+1 

134 



CALL HA0D( DUFFER) 

IJb 



I=DO«MlTrOP) 

136 



buffer=ttop 

137 



TT9P=1 

13a 



GO TO 58 

139 

c 


ASSIGN nESrit'ATIONS TO THE RECORDS HELD 

140 


74 

IsHTOP 

l4l 



170=0 

142 


75 

IF(I)it2.C2.70 

143 


78 

ITO=ITO+l 

144 



FROH(l)=TOtITO) 

145 



I=00aN(I) 

146 



GO TO 75 

14T 

c 


GO FOR REST OF THIS SWEEP 

148 


82 

CONTI’IUE 

149 



IF(AcJStLAST-MEXn .GT.DGO TO 1U 

IbO 

c 


SORTING COMPLETE.., report a OUIT 

131 


88 

CONTINUE 

132 



Print aoa.HECRO.URAbS.NWRITE.NKEAO 

153 


800 

F0RHATtl5.* RECOrfns SORTED*/ 

154 



Ib,* PASSES*/ 

IbS 


- 15.* WRITES*/ 

15b 


- 15.' READS') 

157 



STOP 

IbS 


89 

PRlNl 905.NPASS 

159 


905 

FORMAT t’ STACKS FOULED UP BT PASS*. 14) 

160 



STOP 

lot 



ENO 


o»ejcr 





78 


*•*«** nSOKT tSSOOl 


2Ub37^tHo«0()KSPAC£«( D .5AD|) 


1 



SUiKUUTlNE SA'IDC WHERE) 

a 

C 



3 

C 


PURPOSE... 

«» 

c 


TO ADD A RECORD TO THE * SORTED* STACK 

s 

c 



6 

c 


usage... 

7 

c 


CALLED BY PR0GRA4 toSORI* 

6 

c 



9 



INCLUDE ST0HE5.LIST 

10 



in SURE 1 KNO «INSKM0KN4 1 

11 



XFtSTOPIS.S.lO 

12 


5 

STOP=»HERE 

13 



SDASE=«04EHE 

14 



UP<HI<£RE)=0 

IS 



UOWNt WHERE )=0 

lb 



RETUHU 

17 


10 

IF ICMPARE (HHEPE. STOP > ) 20 r 15. 15 

lb 

c 


ENTER RECORD AT TOP OF STACK 

19 


IS 

UPIST0P»=2H£RE 

20 



UP 1 -HERE) =0 

21 



OOWN(WIIERE)=STOP 

22 



STOP=wHERE 

23 



return 

24 


20 

CHANGEz.TROE. 

25 



IF ( CHP ARE t rtHERE . SP ASt ) ) 2*. 22 .25 

2b 

c 


ENTER RECORD AT BASE OF STACK 

27 


22 

D0iiNCSi3AS£)=*HERE 

2b 



UPtWHt.RE)=SBASE 

29 



D0j)N(hHER£)=0 

30 



SDASe=RHERE 

31 



RETUR) 

32 


25 

I=STOF 

33 


30 

ISDOwfKl) 

34 



IF (CHPARE ( WHERE. 1 ) ) 30 . 35. 35 

35 

c 


ENTER RECORD ABOVE PECORD I 

36 


35 

J=UPCI> 

37 



DORII(J)=VfHERE 

33 



UPt WHERE )=UP(T) 

39 



DOWNC~HEHE)=I 

40 



UPa)=WHERE 

41 



RETURN 

42 



END 


U.EJCT 
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*••**« jSOkT (STOREO) •*«**« 


20b3 /3J1 H*NO}<Kb>' AcCf ( 1 ) . S rOKES 


X 

2 

3 

4 

5 

6 
7 

a 

4 

ID 

U 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 


STONEb PROC 
C 

c proccouke name««* 
c STORES 

c 

C PURPubE. . . 

C TO COMIIWICATE THu OECOKO STACK INPORHATJON FOR SoRTlNO 

C PROGRAM »8S0Rr*. 

C 

C USED BY PROGRAMS... 

C 

C BSORT.BE®Ir»»START*SAOD»SOROP.HADO»HOROP.ADOTO.FINOTO»EPOSHf 

C EPOP 

C 

C LL IS the WMIER of RECORDS HELO XfJ MEMORY, 

C AMO MUST BE AT LEAST 3 

PARAMETER LLS24 
PARAIETtR LGHsaO 
IMPLICIT INTEGER (A-2J 

C DEFUSE KeC0R0{lL0C)=<PtJoJ*IL0CJ#JJJsl.LGH) 

COMHO. • I HCR . St«e . CHANGE » aNOMH . ITO. MTO i PUFFER. 

- STOP.MTOP#FTOP.SBrtSE»HPAS£, 

- OOWMLL) »UP«LL) .FROMtLLi .To(LL) .R<L6H.LL» 
logical SURE, change 

END 


****•« CSORT (START) ♦•*♦♦♦ 


20b373wlH«i*0RKSPACES 1 1 > . START 

1 subroutine start 

2 C 

3 C PURPOSE. . . 

If c TO INITIALIZE storage FOR PROGRAM *BS0RT* 

5 C 

b iNCLUuE STORES. LIST 

7 3UFFEn=i 

6 ET0P=2 

9 DO 10 1=3, LL 

10 10 D0WN(I..1)=1 

11 OOWN(LU=0 

12 STOP=0 

13 58ASE=0 

14 HTOP=U 

15 HBASE=0 

16 INCR=1 

17 RETURN 

18 END 


O.EJCT 
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JSortT (SOrtOP) ♦•••♦* 

eois3;3jiM*rfQKKsrAC«:*ii) .so icp 

1 WTEoLll FuHCTtOiJ *;3K0P(») 

2 C 

3 C PO«pOiE.., 

H C TO DHOP * RECORD FRQV THE *SOHTEO* STACK 

5 C 

6 C USASta... 

7 c called bt pnoerA -1 »ssort* 

0 c 

9 I^CLOCE STORES. list 

10 IFt5UAS£)90rS0,5 

11 5 5DKOP=S0ASE 

12 - S3ASE:UP(SUA5F) 

13 IFCSOASEIlOaO.lS 

t<* 10 STopsa 

15 BETORt, 

15 15 DOHM(Si)ASE)=l) 

17 RE run.. 

la 90 PRINT 901 

19 901 POR'ImH* SORT STACK OVERuROPPEO* ) 

20 RETURN 1 

21 ■' END 


y.CUCT 
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i.sdkT CH'faj 


<!u;i3r3JXM«,aORKS)>ACi:t ( 1 ) .i|AI)a 


1 



SOimOUTl'It llArD(.iitt,Ri) 

2 

C 



i 

c 


Purpose... 

•* 

c 


TO AJO A HFCORJ TO THr 'tlOLO* <;TACK 

5 

c 



6 

c 


USAGE... 

7 

c 


CAtLEO BT PROGKAvl 'aSOPT* 

s 




9 



iNCLUOe STORES. LIST 

10 



1F(HTOP)5.S*10 

11 


5 

riTOP=WH£RE 

12 



HdASE=WHERE 

13 



UP 1 where 1=0 

14 



OOWM WHERE 1=0 

IS 



RETUHt, 

16 


10 

IF ICi-iPARE t where. riHASE ) U5. 20 . ?0 

17 

c 


EMTfR TO bottom OF STACK 

13 


IS 

OOWN t hBASE ) =WHER£ 

19 



UP (WHERE) =HBASE 

20 



DOWN (aHERE 1:^0 

21 



H8ASE=WHERE 

22 



RETURN 

23 


20 

IF(CMP«RE(.tHERe.ilT0P))25.72r22 

24 

c 


enter to top 

as 


22 

UP(HT0P)=WHEP£ 

26 



OOWN(v.H£HEI=HTOP 

27 



up(wheh£)=u 

28 



HT0P=„HEHE 

29 



RETUrtfJ 

30 


25 

I=HTOP 

31 


30 

1=00«N(I) 

32 



IFICMPARECViHERE. I > 130.35.35 

33 

c 


enter ABOVc. ELE'^ENT I 

34 


35 

J=JP(I) 

35 



DOWt.(J)=riHERE 

36 



UP(WHER£1=UP(I) 

37 



□OWN(.;hERE1 = I 

38 



UPtIJ=WHERE 

39 



return 

40 



ENO 

O.EJCT 
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tiSOKT (HPHOP) ***••« 


2Ub37JJIM*«'Jhi;5PACe»( U .HUHOP 


1 



1 ,TEutt{ FUNCTtOM H3ROP(%) 

2 

C 



3 

c 


purpose... 

■» 

c 


TO OKOr> A RECO’'0 FPO ' THE 

5 

c 



6 

c 


usage... 

7 

c 


called by program *bSORT' 

6 

c 



9 



INCLUDE STORES* LIST 

10 



IF(Ht3A5E)90.90.S 

11 


5 

HOrtOP=hHASE 

12 



HBASEsUPtHUASE) 

13 



IFtHrfAsE)10.10.lS 

!<♦ 


10 

HTOP=0 

15 



RETU-^lj 

lo 


15 

DOWNtHdASE)=0 

17 



RETUrtfi 

1& 


90 

PRiNf 902 

19 


902 

FORMAT (' hold stack OVERUROPPPUM 

20 



Re-TURN 1 

21 



END 


stack 


*«••*« oSOHT (AD3T0> **»♦*♦ 


E0b373JXM*rt0KKSPACE5tl) .AOUTO 


1 



subroutine addtoiii 

2 

C 



3 

C 


purpose... 

9 

c 


TO ADD A RECORD TO THE 'TO* STACK 

S 

c 



6 

c 


USAGE... 

7 

c 


called by program 'BSORT* 

8 

c 



9 



include 5T0RES*LIST 

10 

c 


stop if stack is full 

11 



NTO=i4TO + 1 

12 



IF(NT0.5T,LUGO to 90 

13 

c 


FIND PROPER PLACE TO INSERT RECORD 

14 



CALL FINDTOd) 

IS 

c 


HOVE REST records A INSERT NPW ( 

16 



UsNTO 

17 


16 

U=J-1 

la 



IFIU.LE.ITOIGO TO 18 

19 



T0lUn)=T0(U) 

20 



GO TO 16 

21 


18 

T0(ITC+1)=FR0M(I) 

22 



return 

23 

c 


ERROR... STACK OVERFLOW 

24 


90 

PRINT 708 

2S 


708 

FORMAT!* NTO .GT. LL*) 

26 



STOP 

21 



END 


S.EJCT 
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iSOrtT <FI''orC'l 


<!Ub373JlH*i*0KKSPA(,£5< I» .FMdTO 


1 

2 

3 

5 

6 
7 

a 

9 

la 

11 

12 

13 

19 

15 

lb 


C 

C 

c 

c 

c 

c 

c 

c 


bUJKoi'TiM: FI ii)ro<ij 

plhpose... 

SET POINTER MTO' TO POINT IM *TO* STACK WHERE 
« GIWEN NLW BECORO <;H0ItLD »Jf JuSERTEO 


usage... 

CAULLO BY PHOG.^A'^ 'bSORf* 

INClOL'E STORESftlST 
1T0=HT0 
5 ITOsITO-l 

IFIlTo)lb,l5.lO 

10 IF((FR0M(I)-T0tIT0)J*lNCK)5.lS.15 
15 RtrUKH 
END 


**«««* BSORT lEPUSH) 


Z0b373JlM«-0RKSPACEStl> .EPUSH 


1 

2 

3 

9 

5 

6 
7 
a 
9 

10 

11 

12 

13 


C 

C 

c 

c 

c 

c 

c 


SUBROUTINE EP'JSlULOC) 

PURPOSE... 

10 STORE A LOCATION I'l THE »EMPTT* STACK 
USAGE. .. 

CALLED BY PROOHAK 'TSSORT' 

INCLUDE 5T0RE5.LIST 

OOWH1LOC)=ETOP 

ETOP=LOC 

RETURN 

END 


SSORT (EPOP) ♦♦•••• 


20b373JlH*iiORKSPACc«U> .EPOP 


1 

2 

3 

9 

5 

G 

7 

a 

9 

10 

11 

12 

13 

19 

19 


INTEGER FUNCTION EPOP(S) 

PURPOSE... 

TO GET A location FROM THE 'EMPTY* STACK 


USAGE. 


Called by program *bsort» 


INCLUDE STORES .LIST 
lF(ErOP)90.90.5‘ 

5 EP0P=ET0P 

ETOP=DOWN(ETOP) 
RETURN 
90 return 1 
END 


»raeiH4E PAGE B 
W fOOB JSDALinj 
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mhwile 


^Uaj7jJlM**0r«K5PACe5l 1> .PK''FILE 


1 

2 

3 

< 1 - 

Si 

6 

T 

8 

9 

10 

u 

12 

13 

19 

15 

16 
17 
IS 

19 

20 
21 
22 
23 
2*» 

25 

26 
27 
26 

29 

30 

31 

32 

33 
39 

35 

36 

37 
30 
39 

90 

91 

92 

93 
99 

95 

96 

97 

98 

99 
30 

51 

52 

53 
59 

55 

56 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

L 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


NAV£... 

PROFILE 


PUHPObE. , . 

TO UNFOLO the f'ATA INTO IMTEMSTTY VS WAVELENGTH CURVES 


■JSACE..» 


M.pROFtLE OR CXOT .PROFILE 
OPTION?: 

A PRINT AND .(RITE TO OUTPUT FILE POINTS FOR mLL 

MONITOR Intensities (otherwise* if < so* of the 
POINTS a«E ACCEPIEO. thf whole profile is JISCAROEOI 
c DISCARD Points with tihfs for this monitor intensity 
further Than i.6 sicaa {Default 1.73 ) from the nean* 
THUS deleting AROUT lU (DEFAULT 8*) 

E DFLETE 10 UATA POI ITS (OVER<»nrS 'C* option) 

L provide F'JLL LISIINO (OTHEKJISE* only a SUFf*ARYJ 

R report NLSjLTS for all wavelengths (EVEN IP most 
points are oiscahdeo) 

S SKIP SO*’E profiles 

T GENERATE FILES FOR PROGRAM •THEORY* 

SUSPHOGRAMS PEQUIPEO... 

DEGlN*OPT*TIM.iriTLNS*LOOKUP.YE5NO 

MAXIMUM numder oe wavelengths 
parameter LIMEXPS36 

MAXIMUM NUMBER OF POINTS PER WAVELENGTH 
PARAMETER LIMSHF8 

OUTPUT RECORD LENGTH 

PARAhtTER C5I7EF360 
PARAMETER CSI:E2=CSI2 t:+3 

DIMENSION NU«H(HHEXP)*1NTEXP1LI 1EXP*LIMSH) *LENGTH<LIMEXP) 
DIMENSION OLOCK(SS) *EVf*0TE(3) 

EQUIVALENCE ( IHTKEF, ^UOC^ (!))*( TAVER . BLOCK (?))*( J • BLOCK 1 3) ) 
integer TEMP (30)* GRP 

OIMEUSIO.I HEAf!ER(24)*SIG(LI‘-*£Kf’)*ERR0R(LI lEXP) *L«BEL(20) 
eOUIVALENCe(TEMP (9)*HSH0T)*(TEMP ( 5) * AREF( 1 ) ) * (TEMP (13) .BREF(l) ) 
equivalence (LABEL (9) *HSHOT) * <LABEL(SJ . ALIN( 1 ( ) . (LAUELU3) *BLINtI) ) 
DIMENSION AREF(fl)*BREF(8)*ALIfH8)*nLIN(8)*L0CATE(CSIZE)» 
LI5T(C5I7E2) 

LOGICAL SKIP (lO) *PETOFM*PHlNT*OPT. CLOSE* ALL* theory, PRSkIP. every 
LOGICAL teSNO.REPORT 

REAL INTEHS. IUT( 10) *IHTAV(LTMEXP>, LAHDA, length* INTREF 

REAL IBEGIN.ImTEXP 

EXTERijAL IHTEmS 

IHTEGLH poor* early, FOUND *pcent 

InTECCR 6 H0UP(6).F0EQIN{a) ,ENO‘S) ,CR£C*PFILErPREC/l/,CP*TRY*OUT 

integer RSHOT* REFFND *REFI • CEUO'CFILE, El * POINTS * BAD 
integer SHOKLIPEYP.IOJ *TAG(Ll’''EXP.10)*NTAGtl0)*HMTA6(LlMEXP) 
integer TFXLE.TREC.TRIT 
DATA (3LOCK(I),I=9*21)/18*0./ 
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57 

5d 

59 

bQ 

*>1 

02 

03 

04 

05 
66 
07 
oa 
69 
7a 

71 

72 

73 

74 

75 

76 

77 

78 

79 
ao 
ai 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
tuo 
181 
102 

103 

104 
195 
106 

107 

108 

109 

110 
111 
112 
113 


rlfOFlLS 

, 0 .U .--.U. ,, „„„ 

oarA utiTii). ALTER.SPRtiD.Evnorc 
"" ^ 132< .aSf 3«f ' ,, 

call LE0IN( ’profile Z.Ol Jt) ' 

^''f^“-£‘lOl’0,20,lM8RCcJ 
0E^16t FILL CFILE(CS12£|^0,(J,CKEC1 

ALL =OPT(.A OUT- 

CLOSE =OPT(»C') 

EVEHT =OPT(*E*) 

Print soptc’l*) 

REPOr<r=OPT( 'R«I 
PRSKIP20PT( ‘S'J 
THE0RT=0PT('T»> 

IF(CU0SE>SPREAD=2.S6 
1PI.H0T.£VERY)S0 TO 5 
eVNOTtfl»=> I'lC 
£VN0TC(2)=»LU0E0 • 

5 IPt.HLT.THEORTJTPrLEsO 
iDEotr =0.0 
L=0 
K=0 

<• SAVE FILE HE.AnE» 

KEAO UNFItE • 1 > 1 1 TEA DER (I ) , 1 =1 , 1 Z J 
H£ADUNFILE*2> (HEAOEfidj ,I=is,Z4) 

10 READllNFIUE-rNfiEQCAU^"'"'''® 

IF [LAoELt I ) ,E0.E1>60 TO £5 
DECOOL («Oo»LArELKH»P 
000 F0RMAT(A5j 

DECODE I BO 3 * LAPEL) OP 
803 FORMAT ( 16 j 
eP=6P*100 

IFIGrtP.JE.* .)50 Tj) ,2 

t8E6Iri=MAXtt6EeIN.BLlM{i )> 

K-K-t-i ^ 

LOCATttK)=NSHOT 
12 1=0 
IS I=I41 

IFTI.6T.L)bO TO 2(J 

TO IS 

20 IFIUOE.0,00 rS^io? 

LaUi 

GROUP IL)=GP/ 1 Q 0 
F8ESIN(U=I)JHEC-1 
PRINT 805,GHOUPiL).FaEGrN(L) 

805 FORMAT (• FOUND GROUP! *,15,* 

IF(L.ME.l)END(L-l)=m,i£c-2 
60 TO 10 

25 lFfL.N£.0jEND(L)=INREC-2 

1NREC=1 eeOlNNlNS A END OF REFEirtWCE SIGNAL GROUP 

K=0 

00 30 KK=1.L 


AT RECORD' r 14) 
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U4 

lib 

lib 

117 

llH 

119 

1^0 

lai 

uz 

1<!4 

1<!S 

U6 

U7 

126 

129 

IJO 

IJl 

132 

133 
1J4 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 
14b 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 
156 
159 
IbU 
Ibl 
162 

163 

164 
lo5 

lob 

167 

166 

169 

170 


•«»*•* pHOrILE 

30 lrtGI»-U? lrtiO.ro. 0)60 TO 32 
GO TO 191 

32 HtPl^Ft'tGINlKK) 

«EFE.)0=LfWtKK) 

Krt=0 

Kl.IM=L-l 

IneAJ=0 

IFl.i)&T.PRIIJT)PnlMT 306 
306 F01MATO0FILES IMTENSITY POIHTSV 
' (VOLTS) found BADO 

C PRODUCE )!£«(■ FILE FOR EACM LINE PRESENT 

00 100 K=lrKLIH 

c INITIALIZE THE FILE 

3S KKSKK+1 

1F(KK.GT.KLIH)G0 to 190 
lF(GHOUPtrtK).EO.' OGO TO 35 

PFIL£=CHLE4K 
CREC— 1 

JASFUEGINIKK) 

Je=EN0(KK) 

IF(Jj-JA+1.GT.C5IZE) 30 TO 130 
00 30 J=JA.JO 
R£ADU»jFIlE*J)LA8EL 

33 rtRITElCFlLE'CPEOLAQEL 
CENO=lHCC-1 

OD^KJE file PFILE114»3*LI 1EXP4lrU.PREC) 
<RITE1 PFIlE' 1) HEADER 
lMTR£f- = .93*IDEGI). 

TRY=0 

NT=0 

LINES=100 

IF (.NOT. theory) GO TO 40 

OEFIHE file TFILE(13*55+(2FLIMSH)*LIHEXP.U,THrC) 
rfRlTtlTFlLE*!) TFMT . HE AUEK 

C START NEXT MOHITOft SIGNAL LEVEL 

40 lNTREF=ALTER*It)TREF 
TaY=TRY*l 

lF(TRY.6T.20)r.O TO 97 
IF(.NOT.PftSKIP)GO TO 41 
PRINT 013 ,PFIuE.TFILE.X’ITREF 
IFty£5N0(*WA'lTE0?l5*.$97))G0 TO 41 
NT=0 
GO TO 40 

41 AVERH=0.0 

CREC=i 

,JSHor=o 

j=i 

length a )=-2. 

P01NTS=U 

N0M0N=0 

founo=o 

P0OR=G 

SA0=0 

OUT=0 

EARLT=0 

LATE=0 

N=0 



C-X 
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171 

172 

17«» 

175 

176 

177 

i/a 

179 

160 

IBl 

162 

163 

16(f 

16S 

ISu 

167 

168 
169 

190 

191 

192 

193 
199 

195 

196 

197 
19a 
199 
2J0 
261 
202 
203 

209 
205 
20b 
207 

20 a 

2U9 

210 
211 
212 
213 
219 

215 

216 

217 

218 

219 

220 
221 
222 

223 

224 

225 

226 
227 


«•**«• fHOhlL'’ 

1 FI <n -ItA.i « ’jMl-JAhJ ItVlATlO'J OF Tlt-ILU FOf, THIS 

c •M.ITO'I irjrtr.SiTT. 

IFlNl.GT.DoO TO 44 

C •!£ I'UST 60 PEAD THE FILE THE FIrST Tl'lt THkOUGM 

7=0. 

TS=0. 

00 42 L=KtFl<PEFt'lD 
REAOCItjFILE'DTF P 
rA=TM(lrirR£F.AREF,BnEF»»42.1i42) 

NT=Nl+l 
TsTlTA 
TS=TSiTA*TA 
42 CONTINUE 
GO TO 50 

C RECORD TI9£S FOP AU. 'lOMtTOR SIGNALS NOT USED 

C FOK THIS PkOFILF 

44 00 45 L=HEF£rin.PEFl(-l 

CmLL UC0/<UP(L0CATEtL*l-RLFl|,I*LISTr$4S) 

SEADllnFILE'LlT£«-H 

TA=TI'(INTKEF.Ab£P,HREF,145»S45) 

T=T+TA 

NT=HT+1 

TS=T5+TA*TA 

45 CONTI JUE 

50 IFtNT.LE.UOO TO 192 

call EHPTyiLISn 

TAVEF(=T/FLOAT(riT) 

T5I6F.A=50RT UT5-TAVFH*T) /FLOAT (HT-l)J 
TSIG=I (TS-TAW:R«T>/FL0AT<HT-1) l^SPREAn 
T=0. 

TS=C, 

N7=0 

InREC=REF1 

RSH0T=-2 

C 5FT ANOTHER CHRVF U.'LESS WC'BE AT THt, END OF THE FILE 

52 IFICHEC.LE.CE JOJGO TO 54 

c IF we have data for this wavelength, find averages 

IFiN.GT.OjGO TO 70 

C IF WE have data for this TI’-ie Go record it 

IFtJ.OT.lJGO TO 61 
C GO TO next time 

GO TO 9b 

C SET NEK CURVE 

54 READtCFILE’CneClLABEL 
DECODE (812.LAPEL>LA‘10A 
812 F0RH4T(F5.1) 

IFCAJS(LAMDA-LEHGTH(jn ,LT. .021G0 TO 58 
IFdi.GT.OJGO TO 70 
56 UtNGTHC‘JJ=LA*'DA 
N=0 
NF=0 

C GET REFERENCE CURVE FOR THIS SHqT 

58 NF=N1nO(NF+1. JO) 

SHOT(JmNF)=NSHOT 
TAG(0.NF)=> > 

IFMNReC.LE.RFFENO)GO TQ 59 
I=REFl 




niiOKIUE 


2^S 



00 TO 60 

ZiO 


69 

IFILOCArEtlNP«'C + l-«£Pl) ,tO.N5MOT)r.O TO 63 

ZiO 



IrlJltTfC 

Zil 


oO 

00 61 m?EC=hCFl •«£?£: IC 

ZSZ 


61 

IFCLOCArttlNHFC^l-RCFU.K.Q.N'iHOTlCO TO 63 

Zi3 

C 


NO f'OKlTOR SISrnL FOUflO FOR THIS SHOT llU-tOcK* 

Z^<* 



lhREC=I 

2.15 



NO')O.J=HO^OH+1 

2Jp 



TA6{JitJFJ = *N' 

2J7 



GO TO 52 

2ja 


63 

READ C INF ILE • IriREC J TEMP 

239 

C 


RECORD TIME AT f'EXT INTF^SITY 

240 



TA=TI<‘{lNTREF«ALTFR»AR£F.BREFfS64rS69> 

241 



nt=nt+i 

292 



T=T+TA 

293 



T5=TS+TA*TA 

299 

C 


IF TIKE for TMI5 HOUlTOP INTENSITY IS MORE THAN 

295 

C 


1.73 (OP l,6i UNDER «C» OPTION) SIGMA FROM THE '^EAN 

246 

C 


THE monitor SIGNAL IS 8A0 AND WF DISCARD THIS POINT 

297 

C 


THIS SHOULLii assuming A GAUSSIAN DISTRIdUTlON. 

298 

C 


DELETE AUOUT B* (OR 111) OF THE POINTS, 

249 


69 

TA=T 1 • ( INTREF , AREF , bRt F i >66 . S 1Y3 ) 

250 



1F( (TA-TAtfErt)**2.LT.TSIG)G0 TO 65 

2bl 



P00R=P00H+1 

252 



TAG(J.NF)='P' 

253 



IF{.N0T.EVERY)G0 TO 52 

259 

c 


FOUND one mORE curve FOP THIS LEVEL. .RECORD 

255 


65 

N=N+1 

2Sb 



IF(fl,GT.LIii5H)00 TO 1»6 

257 



INT(N)=INTEfJS(TA,ALIN,BLIN.«67.So0) 

250 



SKIP(I,) = . FALSE. 

259 



NTAG(N)=NF 

2o0 



GO TO 52 

2t>l 

c 


INTENSITY IS OUT OF RANGE OF MONITOR VALUES 

2b2 


66 

0UT=0UT*1 

2b3 



TAG(J»NF)=*0* 

2o9 



SO TO 69 

2bS 

c 


TME IS TOO early For signal CURVE 

2o6 


67 

early=early+i 

2o7 



TAG(J.NF)=»E' 

2b8 



GO TO 69 

2b9 

c 


TINF IS TOO LATE FOR 5IG WL CURVE 

270 


o8 

LATE=LArE+l 

2/t 



TAG(J,Nr)=*L* 

272 



IF(LATE.GT. 10)60 TO 97 

273 


69 

H=N-1 

279 



GO TO 52 

275 

c 


found all data for one wavelength... 

276 

c 


CALCULATE AVERAGE INTENSITY 

277 


70 

returns. false. 

278 


71 

Nll=0 

279 



SUM=U 

280 



SUMSU=O.D 

281 



DO 72 1=1, N 

202 



IF(SKIP(I))60 TO 72 

283 



nn=nn+-i 

289 



SUM=SUM*INT(I) 



89 


<•*»•* i>kafiLE 


2ttS 

idb 

ial 

zaa 

^A9 

2V0 

2Vl 

2<#2 

293 

299 

295 

296 

297 

298 

299 
3U0 
301 
3U2 

303 

304 
3oS> 

306 

307 
306 

309 

310 

311 

312 

313 

314 

315 

316 

317 

318 

319 

320 

321 

322 

323 

324 

325 

326 

327 
326 

329 

330 

331 

332 

333 

334 

335 

336 

337 

338 
539 

340 

341 


5U'‘bi-:SU ibOfliiTHl*! UU) 

I,JrtA^ {j*fjN»=iNru> 

72 com’iue 
tlUHSMl J)=NN 

Ik tN.i-l)79,75.73 

73 ir TAVtJ)=SUkVP|N 

SJOCJ) = lSUM50-SUM»lNTAO(J>>/ClN-l) 

c DISCARD Points only oncf 

IFIRLTU4N)00 TO 76 

C TMROU OUT all points FURTHER AWAY THAN l.B SI6MA 

C THIS SHOULD DELETE ABOUT 6St OF THE P0lt4TS. 

DO 74 1=1, H 
IFISKlPtrj )60 TO 74 

IF! (I,jTlIJ-INTAV!jn**2,Lr,SI6IJ)*3.24)60 TO 74 

BAOsHAOkl 

IlFUTAStl) 

TAI5(J, II) = 'ii’ 

IF1EVLRT!gO to 7L 
REfURN=.TRUE. 

SKtPtn=.THUE. 

74 CONTINUE 

C recalculate. A'EAN a standard deviation if 

C WF'VE THRO.v^f OUT A.JY ‘10P=" POINTS 

IFJRETUR.OGO TO 71 
GO TO 7n 

75 INTAV(J)=SUH 
6IS(U)=0. 

76 SI3tJ)=S0ttT(SIG!Jn 

IFItlf.TAV(J) .LT. l.E-6).0o.(N.E0.0))6O TO 194 
ERROR (J)=100.*SIfa(JI/INTAV(J) 

AVERH=AVERH+FRROR ( J| 

LIMTAg! J)=MF 
POlNTS=POINTS+W‘l 
J=J+1 
GO TO 80 

79 IF(RLP0RT)J=J+1 

80 found=foijnd+n 
IF(CREC,GT.CEND)60 to B1 
IFtJ.LE»LINtXi>)GO TO 56 
J=J-1 

print 960, J 

980 FORHAT! 'O**** more THAN*, 14,’ WAVELE'JGTHS, REST OlSCAROfcO *»♦•«/'> 

c 

c DON’T dOTHEH RECORDING OP PRINTING IF MORE 

c than 40« OF THE POIl'TS WFRE THHOWM OUT 

81 FOUNQxFOUtJO’rio.nN+POOR+OljT+FAPLY+LATe 
IF<F0UND.LE.0)GD to 19S 

PCENT=( 100* (NOMOr,+POOR+OUT+SAPLT+UATE+HAOn /FOUND 
J=J-1 

IFIALDGO to R2 
IFlPCENT .GT. 40 ) go TO 92 

C IF wavelength is in UV OIVIGE PY TWO 

C TO COMPENSATE for li'CORRFCT HCPHEARSON SCALE 

82 IF(LENGTHd) .GT. 3300. )G0 TO 86 
DO 84 1=1 #J 

84 L&JGTh{1)=.5*LENGTH(I) 

86 IF(PR1NT)G0 To 87 
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PHOKILE 


3<*2 

34b 

3<»e> 

347 

348 

349 

360 

351 

352 
3a3 

354 

355 

356 
35^ 

358 

359 
3aO 
3t>l 
3o2 
3b3 
3b4 
365 
3o6 

367 

368 
3o9 
3/0 

371 

372 

373 

374 

375 

376 

377 
3,78 
379 

3ilQ 

361 

382 

383 

384 

385 

386 

387 

388 

389 
330 
391 
332 

393 

394 

395 

396 

397 

398 


HHINT 813»PFILE» TFILE* lMl.7ePtF0a«J0tPCCNT 

813 FOHMATIlXf I2> I3»F8.3« 19i 1 5> I"*') 

GO TO 91 

C PRINT PROFILE IP REOUEST^'D To 

87 LINES=LIHlS*J* 17 
lF(I8EAO.t.T.0)6O TO 88 
IF<LlrJES.LT.57)G0 TO B9 

88 PRINI 814rHEAnER 

814 F0RMATUHJ,12A6/IH fl2A6) 

L1NE$=J*'19 

09 print 815»TAVtR*TSIC'4A.INT6EP«PFILE*TFlLe»FOUnO,UOHON» 
- POOR • EVNOTC*OUT » EARLY , t ATE f 8AD . EVNOTE . “O IMTS 


815 FORMAT! IX//' TIME = 

•,F6,J,*+-',FS.3i' MSEC 

MONITOR S*, 


• V FILES',15,* <',I3// 


— 

T1P,I3,' 

POINTS found*/ 


- 

T1D»I3,' 

ri„D wO UONITOP signals 

!N) '/ 

• 

T10»I3,* 

HAD POOR MONITOR SIGNALS 

!P) ',2A6/ 

- 

T10*I3,» 

OUT OF MONITOR RANGE 

!0)»/ 

4» 

T10*I3. * 

Til £S TOO EARLY 

!E)'/ 

•• 

TIP, 13,' 

TIfES TOO LATP 

!L) '/ 

- 

T10,I3,» 

signals too far from mean 

!B) *,2A6/ 

- 

T10,I3, ' 

ACTUALLY llbEO*//) 


print 816 




816 FORMAT!' wavelength 

1 JTENbITY ERROR', 

T50, •POINTS FOUND'/ 


• (ANG5TR0«S> (VOLISJM 


IHEAD=1 

DO 9U I=1»J 

m=LIMTAG(U 

90 PRINT BlB.LLNfiTHtl) •I'lTAVtl) »5lG(I>*EPR0Rf I)* 

«SHOTa»II) •TA&lI,in»Il=l.ITTJ 
313 FORMAT (F10.2»*^9. 5.' (F7,5tFb.0.T3S* 1H»* 

2XflO<I>i.lX.All ) 

A VEIW=AV£RR/FlQA T t J ) 

PRINT SalfAVERR 

821 FORMAT (T17»* AVERAGE SISFA* >F6,C»T35rlH*) 

C 

C WRITE l-jFOrt'-'ATION TO FILFS 

91 IFU.LE*29)LEnGTM{J+1)=0. 

WR I T£ ( PF ILE • PPEC ) I'JTREF , i_e;i jGTH r I ' JT A V » G I Q 
IF ( THEORY ) rtrt I TE ( TFl LE • TRtC ) OLOCk » lU 'S'< ' LENGTH, INTEXP 
60 TO 96 
C 

C REPORT ONE profile OISCARDEO 

92 NOPROF=NOpHOF*l 
IFf.lJOT.PRlHTlGO TO 96 

c GO TO iJExt page if this Page has a profile already 

ifuheao.gt.oiprint ei4 

PRINT 825,1 MTpEF , PCENT « HOMON, POOR, OUT , FARLY , LATE , BAD 
825 FORMAT!' MONITOR =' »F3.-3, 15, 'X DISCARDFO:', 

- I«,' N»,r4,' P',I4»' O', 14,' e*,I4,' O') 

IHEAOS-1 
GO TO 96 

c 

C REPORT ERRORS OCCUpING WITHIH MAJOR LOOP 

c 

180 PRINT gUi,GB0UPlKK>,JA»Jb»CFILE,CSI7E 

901 FORMAT!' group: ',I5»* USES RFCOROS' ,14,' THROUGH', I4» 
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r'HOMLE 


39<J 

4UO 

4U1 

4U2 

403 

404 

405 
40b 

407 

408 

409 

410 

411 

412 

413 

414 
416 

416 

417 

415 

419 

420 

421 

422 
42J 
424 
426 

426 

427 
42B 

429 

430 

431 

432 

433 

434 

435 

436 

437 
436 

439 

440 

441 

442 

443 

444 

445 

446 

447 

448 

449 

450 

451 

452 

453 

454 

455 


'...TOO wAfjr F0>) nur'tiM/ 

- ' WHICH HOLOS* f l4» • RECOFUb*) 

GO TO 97 

192 Pn ZHT 9U2>1HT0EF 

902 FORMAT!* MOHITOR =**F5.3.'t TOO FfW VALID HONITOK SIGNALS*) 

GO TO 40 

193 I=INREC-1 

PKINF 903.1* !TE-1P(I) *I = li4) . <TEMP(I) p Irii7,P) . (TtMP{ I) p 1=2,6 p 2 ) 

903 FORMAT!* INVALID OATA It) RECORD* ,13, * p GROUP *,2A6,*. SCALE*, F6. 3. 

- *, SHOT*, 14/ 

- 1X,6F8.5/U,6F«.5) 

194 FOUNOsFOUND+NOHON.POOR+OUT+EAPLY+LATE 

195 print 904,J.LFNGrnlJl rlNTAVIJ)fN.NM 

904 FORMTI* error OJ P 0IHT*,13,* ( * ,F6. 2 , * ) . . .AV INTENS=* iFB. 5/ 

- * FOUND*. 13 .* points. USING THF FOLLOWING* . 13. IH' ) 

print S 16 ,LEHoTH!J) . INTAv (J) .SIGU) .e'H’OR! J) . 

- ISH0T!J.II) .TAG! J.Ii). 11=1 . NF) 
print 905 

905 FORMAT! IX/* PREVIOUS POIuTS...*) 

PRINT 8l5,TAVER.TSIGHi\.I iTREF’pHFiLE.TFiLE. FOUND, NOMON. 

- I’OORpEVNOTE. OUT. early. LA Tt.CAO.EV'IOTE. POINTS 
GO TO 9o 

196 KKK=LIM5H 

print 906,KKK.J,1 hTAV[J), (SHOTIJ.n .TAr,(j,i) . I=1.,JF( 

906 FORi'AT!* >*.13.* CURVES FOUND FOR POINT ',13, 

- • OF Profile. ..AV intensity =*,F 8 .s/ 

• INCLUDING these: ', 10(15. IX. All) 

60 TO 97 

c 

C TRY TO FINo 10 PROFILES IN ALL 

96 IFIPREC.LT. 13 .AND. (PCENT.LT.40 .OR. PREC.LT.4) )G0 TO 4o 
GO TO 97 

C 

c finish fill 

c 

97 IF(P«EC.LE. 14 ) WRITE (PFILL'PREOEI 
I=FRlC -3 

print 830.1. IOPROF 

030 FORMAT ( T16. 'TOTAL 0F*'l4i* PROFILES ENTERED ArO'.IS.* DISCARDED') 
if! .HOT.T,ic.ORY)GO TO IGO 
IF!TKLC.LT.41 )WRITE(TFILc.'TREC )£1 
TFlLL=TFILfc,+ l 
too CONTII.US 
STOP 

c 

C REPORT errors OCCURING PEFORE )’aJ 0R LOOP 

c 

190 PRINT 910 

910 FORMAT!' NO CURVES FOUNp FOR MAIN SIGNAL') 

STOP 

191 print 911 

911 FORMAT! » NO CURVES FOUND FOR REFERENCE SIGNAL') 

STOP 

197 print 917,U.!GR0UP(r).I=lrL).GP 

917 FORMhT!' more THAN' , 13.* GROUPS FOUND;'/ 

- (bis)) 

STOP 


oeigwaceagbb 

0® POOR gOAlBM 
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PKOHLP 


PROFILE (TIM) 


>.OrtASPACE»(l).Tl^ 

FUNCTION TI.M(INTrA,B* jt$) 

C 

C PURPOSE... 

C TO FIND the time CORRESPOiJOING TO A GIVEN INTENSITY 


T::TIM(INr.A.S.SF0.5O0) 

TIM R.O CAUCULATFO TIME 

INT R.I GlVtN INTENSITY 

A R.I ARRI.Y OF TISES 

H R.I ARRhT of CORRESPOMOiriG INTENSITIES 

t50 S.I EXIT USEH IF DESIRED INTENSITY IS OUTSID*^ THE 

range of the INTENSITIES IN ARRAY B 
SV»0 5*1 EXIT USED IF given qATA IS INVALID 

(I.t. TI»ES NOT INCRCASUG OH Il.TENsITIES HOT 
DECREASING) 

method. . « 

linear interpolation, 
real int 

DIME ISIOU A(8) .n(fl» 

IF((A{ai.LT.A(in.OR.(B(»).GT.Btl)))RFTUR I 5 
IF( (INT.LT.aCBl ) .OR. t INT.GT.Btl) ) jRETlJPN A 
DO 10 1=2*8 

IF(IHT.GE.DU) )G0 TO 12 
10 CONTii.UL 

12 IFtBtU.NE.B(l-U)GO TO 15 
riM=A(I) 

RETUTr 

IS TIM=A{I)-(Ain-A(I-l» )*tlj(I)-IflT)/(HtI)-8(I-l>) 

return 

END 


M.EJCT 
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PUOFtLE IlMirNS} 


. IMC.NS 


1 



REAL FUMCTIOM riTniS(TUi.i A*6>«r*l 




C 





3 

c 


PURPOSE. . . 



*♦ 

c 


TO FIND THE IMTENSITY pOR A OIVEN TIME 



5 

c 





b 

c 


USACc,* < • 



7 

c 





li 

c 


T = lNTEHS(TIHEtA*9.*aO#$90) 



9 

c 





lU 

c 


TIME R»I OiVtN TIME 



u 

c 


A R.I sample times 



12 

c 


8 R«I samples {IMTENSITIESI 



^3 

c 


560 S»0 ERRua EXIT TAKEN IF SIVEN TIME 

TOO 

early 

14 

c 


590 S.O ERROR EXIT TAKEN TF GIVEN TIH£ 

TOO 

LATE 

15 

c 





15 

c 


method . . . 



17 

c 


LINEAR INTERPOLATION 



16 

c 





19 



DIMENSION A16J»4(0) 



<>0 



IFITI' E.LT.AIinRETURM 4 



21 



IF CTII'e.GT.AinM-^ETURil 5 



22 



□0 10 1=2 1 a 



23 



^PITIME.LE.A(I) )C0 TO 12 



24 


10 

cc.mi-uE 



25 


12 

IFIAtll .NE.Atl-inGO TO 15 



26 



INTE 4S=dU) 



27 



RETURl. 



26 


15 

INTENS=a(I)-{ptI)-BtI-in*(Am-TI«EI/(Ani-A(I'll) 



29 



RETURN 



30 



EMU 




U.EJCT 


ORIGMAB PAGE IS 
OE POOR gUALOT 
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««■*•* r>toflLE ILOOKUPI *»*«•♦ 



4C£i(l> .lookup 

1 


SUHhOO rific L00KUP(IJU4'«ER»I»LIST**| 

, 2 

c 


3 

c 

NAME... 

4 

c 

TABLE LOOKUP 

Si 

c 


6 

c 

PUftP05E«.. 

7 

c 

TO CftEATEt ANO later FI ID ENTRIES IN* A SORTED TAsLE OP 

0 

c 

INTE6ERS. 

9 

c 


10 

c 

USAGE. •• ' 

11 

c 


12 

c 

call empty tLlST) HARK LIST AS EMPTY 

IJ 

c 


14 

c 

LIST 1*10 WORK array (SAVED BETWEEN CALLS) 

15 

c 

LISTU) should be set TO THE LENGTH OP THE 

lb 

c 

AHr„Y. 

17 

c 


18 

c 

CALL ENTER (NUMBER* LIST* ») ADD NUMBER TO LIST 

19 

c 


^0 

c 

NUhBER I* I NUMBER TO 8E ADOFO 

21 

c 

S S*I THIS FETURN TAKEN IF LIST IS ALREADY FULL 

22 

c 


23 

c 

call lookup (NUMBER, I *LIST.S) SEARCH FOR NUMBER iN LIST 

24 

c 


25 

c 

number l.r NUFoER TO DE FOUND 

28 

c 

1 1*0 relative location in list* if found 

27 

c 

ZERO* IF not found 

28 

c 

S S*I THIS RETURN IS TA^EN IF NUftiiER IS FOUND 

29 

c 

(Otherwise* normal retupm) 

30 

c 


31 


DIMENSION LIST(50) 

32 


1=0 

33 


IF t LIST ( 2 ) .LT , 3 ) RETURN 

34 


DO 20 I=3*LIST(2) 

35 


20 IFtLlST(l).EQ.NUHRE«)REToRN 3 

36 


1=0 

37 


RETUrt , 

33 


ENTRY ENTERINIIMBER.LIST.s,) 

39 


LSHAX0(LIST(PH-Ii3) 

40 


IF(L.r,E.LIST[lJ»KETUfri 3 

41 


L1ST(2)=L 

42 


LlST(L)=NUMtiEP 

43 


RETURN 

44 


ENTRY EHPTT(LIST) 

45 


LIST(2)=2 

46 


RETUR 1 

47 


END 


((•EJCT 



95 


PKOFTU'' 


dU>>J7jjlM*«OKKbl'A(.ClIU .YtS> lO 


1 

2 

3 

4 

5 

6 
7 
a 
9 

II) 

u 

32 

13 

14 

15 
lt> 
17 
lU 
19 
2U 
21 
22 

23 

24 

25 

26 
27 
2d 

29 

30 

31 

32 

33 

34 

35 
3b 

37 

38 

39 

40 

41 

42 
45 

44 

45 
4b 

47 

48 

49 
bO 

51 

52 

53 

54 

55 

56 


NA't*.. 

YC5N-* 

PDRHOSE... 

TO OrT A TEG OR MO AMbbtH TO A QUFSTION FOR 
AN Interactive proarum. 

calling seouence... 


LOOtCAL TESNO 


IF(VESN0(*SHALL I SKIP OPTIONAL PART5I3»,S90) )G0 TO 40 
20 continue 


<OPTIO.JAL PART OF PROSRA' > 


IF (YESMOI 'SMALL I OEPEAT OPTIONAL PAHT?t')J60 TO 20 
40 CREST OF PR0GRAM> 


90 CSPECIAL section FOR END-OF-FILC RETI)RN> 

THE QUESTION MUST ENO WITH THP STOP CHARACTER 8 {WHICH 
ISN’T PRINTED) » 04 THE QUESTION WILL 0E FOLLOWED BT A 
WHOLE STRING OF TRASH* AND THF PPOORAM F.AY BLOW UP. THE 
QUESTIO') WILL nt PPINTEOt AND THE TELETYPE wILL WAIT OH 
THE SA{£ LIME FOR THE USER TO TT°E HJS ANSWER. THAT 
answer .AAT start IN ant column AMD CONSIST oF At.T NOPHAL 
AFFIRMATIVE OH NEGATIVE WORD CI.C. r ANT ONE I COUUO THIi.K 
OF HlEJ I wrote IT). the USER SUPPLY AN 5.E0F RETURN 
ADDRESS AS THE SECOMO ARGUHENT, PUT THIS IS OPTIONAL 
(The PROGRAM will FIUU IT’S VAT 3ACF EITHER WAT). 

CONTENTS OF ALL REGISTERS EXCEPT AO ARE SAVED. 

ERFOR COHDITIOMS... 

IF the USER'S AliSWER IS A PLANK LINF OR A CHARACTER 
string the PROGHAv DOESN’T RECOGNIZE (IT CAN BE FOOLED) 

IT TdIES again, printing OfILY ’WHAT?’. IT WILL REPEAT THE 
ORIGINAL OUESTION O'IE TIME IH FOUR. 

IF THE RFSPO.JSE IS 3E0F* THE PROGRAM USES THE ALTERNATE 
return address* IF SUPPLIED. THIS CAN BE USED* FOB 
INSTANCE* TO terminate THE PROSHAH OR RE-ASK A 
PREVIOUS QUESTION (THIS JAY YOU DON’T HAVE TO RE-RUN THE 
WHOLE P*<0GRA4), IF THE A ISWEP IS ANY OTHER CHARACTER 
STRING BEGlNJllMG WITH 0 THE SYSTEM THINKS THE PROGRAM 
IS TRYING TO READ A CONTPOL CARD, AND ALLOWS NO FURTHER 
READS. IN THIS CASE* OR IF OEOF IS ENCOUNTERED AND NO 
OEOF RETURN ADDRESS HAS BEEN SUPPLIED* THE PROGRAM PRINTS 
A SHORT MESSAGE AND FXITS. 
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I'ROHLE < yes 10 1 


or 

/ 





58 

s(0) 

LIT . 




59 

SU) 

, 




60 

ftSNO* 

SA 

XllfPETURH 

» 




DS 

Al.SAVF . 



bZ 


06 

A3»SAVF3 . 



bS 


Sr< 

RltSAVERl , 



64 


S'< 

R14*SAVER1A 

♦ 


65 


ESDIT 

PACKET . 



6b 


LX 

Xll.RETURN 

• 


67 


L«H1 

AOflfXll 

« 

GET Ft^T WCRO FOLLOMING QUESTION POINTER 

68 


AN.U 

A0» 0742000 

• 

Compahc to JUMP instruction 

69 


S 

AO. PRESENT 


THIS INOICATES PRESENCE OF OEOF ADDRFSS 

70 


L 

AO.O.Xll 

• 

«£T AUOSgSS OP USER'S QUESTION 

71 


L>4J 

XlUEMSfiS 

• 

*N0 insert it into PrlNT LINE 

72 


LMJ 

X11.E0ITX9 

• 

END fU|T MODE S RESTORE REGJSTEkS 

7J 


LiU 

A4.3 


IWITlALIrt OUERT COUNTER 

74 

ASK 

LtU 

AO.REAOPXT 

* 


75 

B£A0 

L>U 

ai.rlamks 

f 

FILL ItIP.fT LINE KITH blanks SO THE 

7b 

* 

LXI>U 

A2.1 


PREVIOUS INPUT UOESfHT CONFUSE THE iSSUw 

77 


lx.<i»u 

A2. INPUT . 


78 


LH»U 

R1.14 . 



79 


at 

A2.0.«A1 . 



80 


EH 

TPEAOJ 

• 

DO A TreaO* to print OUtSTION S GET ANS, 

81 


LSSL 

AO. 4 

• 

15 REPLY ACTUALLY IN 'INFOR' FORMAT? 

82 


JH 

AO.EXA'tiNE 

• 

REPLY IS PRES*^NT 

83 


LA 

AO.RFAnPKT.-l 

. INP'JT is in 'INFOR' FORMAT... REREAD 

84 


EK 

REAOS 

• 


85 

EXAvitjE 

LX UU 

A2.0 

• 

INITIALIZE index REGISTERS 

o6 


U»S1 

A3. INPUT 

• 

IF FIRST CHARACTER IS NON-BLANK, SKIP 

87 


L 

AO. INPUT 

♦ 

TO SEARCHINS POR A MATCHING hORQ 

as 


T£.U 

A3. 5 

• 

OS IS A rLANK 

89 


J 

HATCH . 



90 


L-<»U 

Rl.lA 

• 

0K...WE'LL have to do IT THF LONG WAY 

9l 


LA 

Al.BLAMnS • 



92 


Si4E 

Al.INPllT.*A2 

, FIND FI?ST WORD NOT ENTIRELY BLANKS 

93 


J 

AGAIN 

* 

OOPS.., THE WHOLE line WAS BLANK' 

94 


L 

AO,tNPUT-l. 

A2 

. GET THAT FIRST NON-BLANK CHARACTER 

95 


u^u 

A3.5 

• 

FIGURE OUT WHICH CHARACTER IT IS 

96 

LOOP 

ALU 

AO. (07-7(J00000i 

uOOO) . WAqK 0"T LAST S CHARACTERS 

97 


T£ 

Al. (005000000UOCO) . IS IT A BLANKS 

93 


J 

LOAD 

• 

NO... 00 look for a I'ATCMING «0RL 

99 


LSSL 

AO. 6 

• 

TES...TRr next character 

lUO 


JSD 

A3. LOOP . 



101 

LOAD 

ot 

AO.INPlJT-1. 

A2 

, GET THE MESSAGE 

102 


EX 

MOVE. A3 

* 

SHIFT the blanks AWAY 

103 


LXM»U 

A2.0 . 



104 

HATCH 

LR«U 

Rl.WnCOUMT 


search for matching word in table 

105 


SE 

AO. table. *A2 , 


106 


J 

AGAIN 

• 

'no hatch... ASK again 

107 


L 

A0.A2 

« 

FOUND A natch. Get index 

lua 


LSSL 

AO. 35 

• 

LAST dIT IS 1 IF 'YES'. 0 IF »N0* 

109 


SSL 

A0.3S 

» 

{SINCE A2 HAS BEEN INCREMENTED) 

110 


SLJ 

RESTORE 

f 

RESTORE registers A RETURfl 

111 


U 

PRESENT , 



112 


J 

2. XU 

» 

RETURN TO CALL *2 IF HO OEOF AODR SUPPLIED 

113 


J 

3.X11 

• 

RETURN TO CALL+3. SINCE 8£0F AODR PRESENT 
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PHOPtL'T (V£!)iO> 


114 

AGAIN 

L.U 

lAA 


Jo J 

Uo 


L.U 

117 


J 

ita 

EOF 

SLJ 

119 


T,»t 

120 


JP 

121 


L 

122 


EK 

123 


C« 

124 

RESTORE 

HES 

125 


OL 

126 


DL 

127 


lr 

128 


LR 

129 


LX 

130 


J 

131 

/ 


132 

sto> . 


133 

WBCK 

•YESNO. 

134 

RETJH!) 

+ 

135 

SAVE 

RES 

136 

SAVE3 

RL5 

137 

SAVER! 

RES 

13S 

SAVER! 4 

RtS 

139 

PRESE’ir 

+ 

140 

PACKET 

EIPKT 

141 

• 


142 

. 


143 

PRTLlNE 

•control 

14+ 

prtcnt 

tOU ; 

145 

PF 

FORM 

146 

reaopkt 

PF 

147 


+ ( 

148 

LINE 

RES J 

149 

AOAINPKT 

Pr ] 

ISO 


+ I 

iSi 

QUESTION 

••HAT?* , 

152 

INPUT 

RES j 

153 

ULAT'KS 

, • 

154 

. 


155 

HOVE 

LDSL t 

15tt 


LUSL / 

157 


LOSL t 

158 


LOSL t 

159 


LUSL P 

IbO 


NOP . 

Ibl 

table 

•YES' 

162 


'HO' . 

163 


•Y* , 

164 


• N» , 

165 


•OUI' . 

166 


•NON* . 

167 


'ALRICH' 

16S 


.NOPE* . 

169 


•YEP* . 

170 


•NAW' . 


AlUft'lAlflPKT 
rtt.REAO 
A4.3 . 

ASK . 
RFSTOKr 
PHfSCMT 
AOi*l#Xil 


A0,(PF ;fP'?TCuT,PRTLI»'C) 


p-'T' r nr f.jTint. q'p ;tiou oit-t one Tiotr 
Out nv Four 


CJ3-CF-FILF SFTUWI.., RELOAD HEGlSTEH'i 
R:TII0 . QUir tF i«t.0F HOOK SUPPLltO 
RETURN TO OCOF address 


pniNTi 
EXIT'S 
1 

AlrSAVe . 
A3rSAVE3 . 
RlrSAVERl . 
Rl4»SAVEBl*f 
XI If RETURN 

•restore . 


0 . 

2 . 

2 . 

1 , 

1 . 

-1 . 

24, LINE 


, NOfl""EOF CONTROL CARo ENCOUNTERtO.., 

. RE'RE FORCFO TO OlItT 

.•routine for RESTORiriG REei'STERS 


I-PRTLINE 
12,6, la . 
1.22, LINE 


. JF THE USER FOoGETS THE STOP 

character, a 0 is the character host 

LIKELT To BE found By ACCIDENT. 

■'"S PR06RA-! EXIT' . 

, hu^ser of words in PeSSAGE 


24 . 

1,1. question 


14 


this ensures the SECaUO .oRO lOhOEO "AS 

BLAHPS IF THr LAST WORD Ifi INPUT HAS INSTR. 

Instructions for shifting out blanks 


even relative locations for 'YES*. 

ODD FOR 'NO' 


OSIGINAB'^AGB IS 
OF POOR 
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17l 
1 7a 
17J 
17<* 
17b 
l7o 
177 
17b 
17-J 
IbU 
Ibl 
182 

183 

184 

185 
Ibo 
187 


•jihO.lj* . 
•HIGriT' . 

*JA* . 

’■'.YET* . 

•DA» . 

•JUNT’ . 

•DO* . 
tJOI4**T* . 

»OK» . 

•HOT* . 

•okay* . 

•N* . 

*51* . 

wOCODM EOU S-TABLE 

E.^0 . 


(’ROFtLF (YESnO) 


f T>'E MUMBEH OF WORDS IN THE TAOLt 


M.EJCT 
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Vf»tOT 


auMrajiH^^oiwspAce* < i > . vj’lot 


I 

C 



2 

c 


I4A IE... 


c 


vplot 

•» 

c 



S 

c 


PURPlibE--. 

b 

c 


TO PRODUCE A PRlNTe«-PLOT OF THE EXPERIMENTAL LINE PROFILE 

7 

c 


PROOUCFD DY program 'PROFILE'. 

a 

c 



9 

c 


USAGE... 

la 

c 



11 

c 


OXQT .VPLOT OR J. VPLOT 

12 

c 


<FILE «> IFREE) NtJllBER OF 1ST FILE TO PLOT 

13 

c 


<FILE tt> (FRtEJ number of 2nO FILE TO PLOT 

19 

c 


BEOF 

15 

c 



16 

c 


options: 

17 

c 


T PLOT IS TO RE OM A TERMINAL* SO IT USES ONLY 

18 

c 


columns 1-72. 

19 

c 



20 

c 


Input ... 

21 

c 


OATA FILE AS PRODUCED OT PROGRAM 'PROFILE' 

22 

c 



23 

c 


subprograms R"^ouipeo... 

2<f 

c 


BEGIN. OPT.NUMBtH 

as 

c 



26 



PARAMETER HMEXP=?6 

27 



LOGICAL OPT.RaTCH 

28 



REAL iNTf 0:*LrN5rH*lMTAV.lMAX*lNCR.NL 

29 



aiMENSIOfJ LINE(IOS) .HEAD£H(24) *LER0TH{LIMEXP) * INTAV (LlMEXP ) . 

30 


- SIG(lWEXP1 

31 



OATA tI/''OEOF'/.HARX/*!'/.lZERO/*0'/,I ONE/ * 1 ' / 

32 



JCOLlXl=l+IFlX(X*IUCR) 

33 



call t.E6lNt*VPL0T 1.12 O'J 

34 



IDIF=I0NE-12EP0 

3S 

c 


find file »ND READ FILE HEADER 

36 


5 

NFrLE=NU IbERCFILF NU bCrt’T .7*29.599) 

37 



OEFINC file f<FILE(14.1+3«ui'!EXP.Uptl'>EC) 

38 



READ tuFiLE' 1 ) HEADFR 

39 

c 


OECIOE flHETHE" S-'ALL CT' OPTIOf) OR LARgE GRAPH 

40 



IC0L=4l 

41 



BATCHi.NOT.OPTCT*) 

42 



IF<8ATCH)IC0L=lfll 

43 



COL=FLO AT (I COL-1) 

44 



LASTS' * 

45 



I F < B ATCH ) LA5T=M AR R 

46 

c 


read a PROFILE 

47 


a 

1FINREC.LE.14)G0 TO 11 

48 



1=MAX0|U.NKEC-2) 

49 


10 

PRINT 805,1 

50 


ao5 

F0RMATUX//I3.' PROFILES') 

51 



PRINT 808 

52 


008 

FORMA T(lHl) 

S3 



GO TO 5 

54 


11 

READNJFILE'NREOlHTMOtl.LENGTH.lNTAV.SIG 

35 



IFUhtmon.eq.edgo to 10 

56 

c 


SET UP SCALING PARAMETERS 



100 


vflvjt 


tif 

59 
bU 
61 
62 

03 

04 

65 

06 

07 

08 
p9 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

60 
61 
62 

63 

64 

as 

66 
87 
6B 
69 

90 

91 

92 

93 

94 

95 

96 

97 
96 
99 

too 
101 
102 
' 103 

104 

105 

106 
107 
106 

109 

110 
111 
112 
113 


ir'AXSb. 

.JL=1U. 

'I— J 

12 lFi(LeN»TH(H*l) ,LT. lC-5) .OR. tM.GE.LI"CXP) JGO TO 15 
ilzHi-l 

IF«fi.oE..2)«L=ArUi^UNL.LEi.0TM{ I)-LE\3Tl I tfJ-1 1 ) 

IMAX=A ^AXl ( 1«AX, IMTAVIN) +SIS (4) ) 

GO TO 12 
15 NL=1./NL 

IFM.’ 0 .*tt-ENGTH{rg)-LEHGTHll)).LE;. ?4.).AND. BATCHJNL=NL*2 

C 

C NU IS NO-g THE SCALE IN PRINT LINES PER AN55TR0M 

T0TAL=NL« ( LEHRTH ( N > -LENGTH 1 1 ) ( 

IFITOTAL.GT.200.JGO TO 115 
POWERsl. 

1JP=0 

20 IFtINAX*POWER.GT.l.)GO TO 25 
POWEH=POHErt*10. 

NP=NP-1 
GO TO 20 

25 lF(INAX»POf(ER.LC.10.)00 TO 30 
POWEH:PO^ER*.l 
,gp=NP+l 
GO TO 25 

30 IFtIMflX*HOhER.LE.2JG0 To 3o 
lFtHAX*POWER,LE.G.)GO TO 34 
UT=10 
GO TO 40 
34 NT=5 

GO TO 40 

36 NT=2 

C PRINT graph HEAPER 

40 PRINT dIO, header* INTMON. uFlLE 

810 FORMAT UH1.12A6/1H .12A6//' MONITOR INTENSITY =»»F6.4.» VOLTS'* 

ax, 'FILE' *13/) 

INCR=COL»POWEP/FLOAT < rJT ) 

PHEV=LENGTH(1) 

C PRINT TOP UOROER 

00 42 I=1.IC0L 
42 LINEU) = '“' 

L=IC0L/'1T 
00 44 1=3*1C0L*L 
44 LINEU4i) = ' + * 

PRINT Sll*(LItl£(J).J=l*ICOL) 

811 FORMAT C27x*10SAl) 

DO 4o I=1*1C0L 

46 LINE{I)=‘ ' 

C 

c INSERT NEE6E0 BLANK LINFS 

DO 60 1=1, N 

IF<SIG(I) .LT. O.OJGO TO 60 
LINES=lFIX(NL*tLEhGTH(I)-paEV)+.S) 

PREV=LCNGTH(I) 

IFILINE5.GT,25)G0 TO 50 
48 IFfLIllES.LE.DCO TO 55 
HHES=LINES-I 
print 812.MARK.LAST 





VPLuT 


U.4 

tlb 

llD 

117 

llu 

119 

1<>0 

121 

122 

12J 

129 

125 

12b 

127 

128 
129 

i:so 

tJl 

l.}2 

133 

139 

135 

13tj 

137 

13B 

139 

190 

191 

192 

193 
199 
145 

196 

197 
196 
199 

150 

151 

152 

153 

154 

155 

156 

157 
156 
159 
IbO 
161 
Ib2 
163 
lo9 


Hi2 roKVAT(27X*Alt99AfAl» 

5C ro 98 

50 PRINT B12» (rtAUh.LASTi Ilslr9) 

PRINT 613.UA5T 

813 F0ltiATt25X»*/////*iTl?8iAU 

PrtiNl 812* (MARK»UAST.Xl5li9) 

c 

c PRitiT ONE ul.JE OF GRAPH 

55 LINE<U=HARK 
U JEUCOUSLAOT 
Jl=JC0l.(INTAV(n*SI6(I) ) 

J25JCCL I AMAX 1 [ 0 • INTA V (I > -S lG(t ) > 1 
00 57 JSJ2*J1 
57 LIH£tJ)='*» 

PRINT 815. INTSV r I ) .SIGI 1 1 .UENrTKC I ) , (LlNE( J) * J=l* ICOL) 
315 FO‘IMATtlX.F7.5* ' +-' *F7,s*F9.2* lOlAU 
00 59 J=J2*J1 

59 LINE(dl=* • 

60 CONTINUE 

c print LOKER BOPnEP 

00 62 I=l.ICOL 
62 LIH£tH = »-* 

00 69 1=0*IC0L»L 
69 LI JEU4l) = ' + * 

PhINT 811. (LI JE(J) *J=1*1C0L) 

DO 66 I=1,IC0L 
66 LIM£(I)=* » 

C PRINT HORIZONTAL SCALE 

00 68 1=0 .NT 

68 LINE(U*I+1) = IZER0M»I0IF 
IF(NT.NE.10)G0 to 69 
LINE(ICOL)=IZrRO 
L1NEUCOL-1)=ION£ 

69 IFI.'JP.EO.PIGO TO 70 
L1N£(IC0L+1)=' * 

LINE(IC0L+2)=*E' 

LINEUC0L+3)=' • 

IF(NP.LT.O)LI <E(lC0L+3)=*-* 

LINetIC0L+4)=l2EH0+IAdStliP)*lnlF 

NCOL=ICOL+9 

print 811. lLlNEtJ).J=l.NLOU 

GO TO a 

70 PrUT all. <L1'|E(J>. J=l.ICOL) 

GO TO 8 

C OUIT 

99 STOP 

c error messages 

115 print 915. total 

915 FORMAT!' GRAPH LENGTH 0F.'»F6.0>’ LINES IS TOO LONGJ'J 
GO TO 8 
ENO 


8.EJCT 





102 


1 

2 

3 

5 

6 

7 

8 
9 

10 

U 

12 

18 

14 

15 

16 
17 
10 

19 

20 
21 
22 
28 

24 
2b 
26 
27 

25 
29 
80 
31 
82 
83 

34 
8b 
86 
87 

35 
89 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 
So 


• **»«* THEutiY 


• «OHKbI'ACC»(l> .THEOKT 


C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


riAME... 

theory 


purpose.. . 

ro Fi'io the values of the parameters (line position, 

LINE INTENSITY, LINE HIUTH OR ELECTRON OENSITY, A|,0 
BACKGROUNU INTENSITY) WHICH BEST FIT THE THEORETICAL 

profile to THE experimental data. 


USAGE... 


B.THEORY or OXOT .THEORY 
<DATA> 

options; 

•U» DELETE SOME POINTS 

•N* PRODUCE Plots mow 

•T* PRODUCE files FOR PLOTS LATER 

(FEATURE UhTESTEO) 

»S* SKIP SO‘E profiles 

*Y» PRINT DIAGNOSTIC INFO FOR EVERY POINT REOUtSTEO 
BY 'ZXPOWL* 

'Z* PRINT diagnostic INFO EVERY TIMF *ZXPOwL' 

CHANGES SEmRCH DIRECTION 


DATA I I AGES; 
card l;tFHEE) 
CARD 2: (FREE) 
CARD 3; (FREE) 
CARD 4: (FREE) 

1 

2 

3 

4 

5 

CARD 5; (FREE) 
CARD 6. (FREE) 
(REPEAT cards 


INPUT FILE NtlPDE'’ 

output file NU^UFR, (ONLY WITH OPTION) 
WAVELENGTHS TO DELETE (OfLY WITH *0* OPTION) 

LInl number 

HE II 46S6 (density 1.E17) 

HE II 1640 
HE II 1215 

HE II 1025 

HE II 4606 (DENSITY I.E10) 

•yes* if fit to THIS PROFILE IS WANTED 
(Only with * 5 * option) 

•YES* IE PLOT BESIRED, »f'9» IF NOT 
(Only with *n* option) 

5 S 6 F0“ REGAINING PROFILES AS NEEDED) 


SU5PK0GRAMS REQUIRED... 


fetchs obtains The nor ulized line profile and theoretical 
line position FOP the lime. 

ZXPOWL (FROM IHSL) GENEPAL FUNCTION MlfJlMIZER, UStO TO 

FIND the riOf.LINEAR PARAPETERS WHICH MINIMIZE THE SUH 
OF THE SQUARES OF THE DEVIATIONS OF THE EXP£RIMEHTAL 
points from the THEORETIEAL profile 
FUNCT 3 HAS THE THtORETICAL PROFILE CALCULATED FOR THE 
CURRENT PAKAmETFHS, AND O0TAINS THE MEAN SQUARE 
deviation of THE DATA FPOH THE MEW PROFILE, 

FUNCT2 finds THE CORRELATION MATRIX OF THE BEST-FIT 

PARAMETERS, USING THE SECOND PARTIAL DERIVATIVES 
OF THE error function. 

TPLOT PLOTS THE EXPERIMENTAL DATA AND THEORETICAL PROFILE 
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TMCoiiir ♦*•»♦♦ 


b7 

C 


r,.E oFF-LI IE PLOTTi-ft. 


ba 

C 


AXISN»REbI'l»FETCHb, JF*«;tUE VT.Ne'WU.0Pr,PLmc*Sll>HA,SY'15LV» 

b9 

C 


TP|_0T. VALIJP 


oO 

C 




d1 

c 


method... 


62 

c 


I'’SL HOUTIMF <ZXpo.JL' FltlOS THF MIHIMltn OF THt ERKOO 


o3 

c 


function (THE SUM OF ThE S3UAKFS OP THF OEVIATlONb OF 

THE 

6l» 

c 


fitted FUlCTIO'J FK01 THE DATA POINTS) AS A FUNCTION 

OF 

6S 

c 


ITS T»0 NO l-LlNEAh PARA'*ETFRS. THE COVARIANCE MATRIX 

IS 

66 

c 


COMPUTED AS THE IiiVERSE OF THE MATRIX OF SECOND PARTIAL 

o7 

c 


DERIVATIVES OF THt ERROP FUNCTION NEAR ITS MINIMUM. 


68 

c 




69 



INCLUDE DQANK.LIST 


70 



INCLUDE GROUP .LIST 


71 



COMMON /F3C0.M/ ARG2(2( 


72 



ExTEKt AL FUNCT3 


73 



heal max 


74 



LOGICAL TESNO , KEEP , OUT»OPT» LATER , ALL 


75 



DI XENS I ON HE APER 1 ?4 > . H 1 3u ) . ARGO ( 4 J . DAO 1 25 ) 


76 



D1 lEHSIUN LARFL(2n) 


77 



INTEGER 0UTFIL»OJTFMT.0UTrEC.''FILE.DRFr,TF-T 


78 



data EST. CLOSF»CLOStP. E 1D» FPSNE. EPSI.QUTFMT, tFMT 

79 


- / l.E-4, l.E-4» 3.E-6**(3=’0F ** .2* .02i 700» 

503/ 

80 



DATA LLviEP 


ai 


~ / 0/ 


82 



DATA (LABEUI) .I=1.10)/*mE II ’*'4686 »»*HE II '.*1640 *, 


83 



*HE II »,»I215 »i»rt£ II '.'1025 ','HE II 'r*4bP6 •/ 


64 



IPCeNT(00) = IFIX(IOO.*AtSSlQO)*.b) 


65 



CALL UEGIl.t 'THEORY 1.74 jt ) 


86 

c 


GET options FR0« CO ITPOL CARP 


67 



OUT=QPTt 'T' ) 


88 



LAIEH=.NOT.OPTt 'N* ) 


89 



ALL=.NOT.CPT('S') 


90 

c 


SET UP PARA"£TERS 


91 


10 

IFtLINE?.N£.0)PRIMT 820 


92 



DF1LE=NU'!B£R( • INPUT FILE?9' .7.29.$96) 


93 



OE'^INE FILE 0FILEtl3.ISlEC.'J.ohEC) 


94 



IFIOUT)OUTFIL=NU’ RERl'OUTfUT FILE’D' »7>29ii96) 


9b 



IF(OUT)DEFINE FILE OUTHL 1 13. ISI2E.U, OUTREC) 


9b 



KEEPb.MOT.OPTI'O'I 


97 



IF(Ki.EP)GO TO 15 


98 



PRINT 300 


99 


000 

F0R.<ATC ENTER WAVELEIGThS TO DELETE •) 


lOQ 



REAOt3.810.END = 13) taAOCfldAOI r«'bAD=X.25) 


lOl 


310 

FORMA TO 


102 



NBAD=26 


103 


13 

NBAD=t.DAD-l 


104 



IF t NBAD . EO. 0 > KEEP= . TRUE , 


lOS 


IS CALL FETCHStABGOMJArLlJJE) 


106 



READ(DF1LE»1)hFPT.HEA0EA 


107 



IF(NF.1T.NE.TFmT)G0 to 97 


108 



IF I OUT ) »>R I TE ( OUTF IL » 1 ) OUTFIT » HEADER 


109 



print 320. header 


110 


820 

format (1H1.12A6/1X.12A6) 


m 



LIN£P=2 


112 

c 


read data for NFXT monitor intensity 


113 


20 

READ ( DFILE* DREC ) BLOCK. NSHOT .LENGTH. INTEXP 
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theory 


U4 

US 

U6 

U7 

118 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 
129 
13U 

131 

132 

133 

134 

135 

136 

137 
13» 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 
IsO 

151 

152 
lt>3 

154 

155 

156 

157 

158 

159 

160 
161 
162 
lo3 

164 

165 

166 

167 

168 

169 

170 


,JAR6=NA 
00 22 1=1»MAK6 
22 AR4(I>=AKG0(T> 

LINE5tl>=LINF 
IFtlNTi-l0N.EO.rNO)4O TO 10 

C CHECK FOR acceptable WAVFLEfJGTH RANGE 

IF(AKG(2)'f9ASE«LT.2«LENGTH(l) -LENGTH ( rjU ) . OR . 

- ars(2u-ba5E.gt.2»lehgth<nu)-lengtm(ii) go to 92 
NU1=1 
NU2=NU 

IF (KEEP) GO TO 50 
NUl=C 
NU2=HU+’1 
38 NUl=NJl+l 

1F(NJ1.6E.NU2)60 TO 49 
00 40 l=ltl'3A3 

1F(AUS(LENGTH(NU1)-BAU'(I) ).LT. .001)60 TO 38 
40 CONTIfWE 
42 NU2=HU2-1 

IF(N01.v>E.NU2)GO to 4" 

00 44 I=l,N6AP 

lF[AbS(LEN3TM(N02)-BA0(I) ).LT. .OODGO TO 45 

44 CONTINUE 
GO TO 42 

45 HOUD=LEHGTri(MHl) 
length [ .JUl ) =LENGTH ( NU2 ) 
length (NU2)=H0l-D 

N=HAXO ( NSriOT ( MUl ) . MSriOT ( NU2 > ) 

00 48 KSH0T=1.N 
H0LD=1NTEXP (Mill rKSHOT) 
lNTEXF(NUl.K5H0T) = I(lTEXpiN0?rKSH0T) 

48 UTEXP(NU2»K^H0T)=ri0L0 
N=NSi«0T(NU1) 

NSH0T(NUl)=NSH0T{flU2) 

NSH0T(HU2)= J 

60 TO 3U 

49 NU2=NU 

c Print labeling information 

50 TOTAL=0 
90TT3II=LENGTH(NU1) 

TOP=LENGTh(N'J?) 

AVG=0. 

DO 51 KU=NUl.flU2 
0OTTO”=A'IirJltOOrTnH»LENGIH(KO)) 

TOP=Ai!AXl( TOP. LENGTH (KL) ) 

N=NSH0T(KUJ 

total=total+h 

00 51 KbH0T=l,N 

51 AVG=AV6*INTEXP(KU.KSH0T) 

AV6=AVG/T0TAL 
IF(AVG.LT.1.E-S)G0 TO 95 
LlHEP=LlNEP+n 
IF(LINEP.LT,57-3JG0 TO 53 
print aaO.hEAOER 
L1N£P=15 

53 PRINT 83Q.TAVEB.IflTMON.TOTAL.flOTTOft.TOP 
830 forhatux/' tihe;*.f6.3» • monitor: '.F e. 4/ 
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****** TMtvi'ilf 


17i 


- 

1X»I3** ■’OliiTSi #A»ELE .GlHS: '.FP.2,* TOifFd.21 

172 



IFlAULIGO TO S5 

173 



lF<rt’j40l»FtT’-')49b))60 TO 55 

174 



LINEH=L1N£P-U - ^ 

175 



uO TO 20 O 

176 

C 


CHECK TiIAT the system I5N*T oegeherate 

177 


55 

I=NU2-NUl«'l 

176 



IFU«LE«'IARG*2)G0 to 93 

179 

c 


SET UP for solution 

IdU 



meval=o 

131 



DO 6J 1TER=1i6 

132 



A1=ANG(1) 

183 



AfiS2U)S9.9 

134 



ARG2t2)=9.9 

135 



CALL hEFSTARG.NAKGI 

136 



ul <IT=iO 

l«7 

c 



188 

c 


HAVE MINIMUM FOUND (USING KOuTItlE FROM IMSl) 

139 



CALL ZXP0liLlFliNCT3fCL0SEiNA30,flR(j,v4LUE»LlMTT,Hi XER) 

190 



IF(dPT( »f)» > )PrI>ll 760. AR0( 11 .flH6(2J ■ VALUE* IE'» 

191 


760 

Fo.MATt* results: ».3Gl3,u. 141 

192 

c 



193 

c 


IF ELECTRO,. DENSITY CHANGED MUCH, REPEAT SOLUTION 

194 



IF(AdS(ARG<l)-Al) .LT. ,1*ARG(1)>G0 TO 65 

195 


60 

CO«'JTIf,U£ 

19o 



1TEP=d 

197 



A1=Ai(G(1U1.E16 

196 



PhINI 838.ITFP.A1 

199 


338 

FORMAT!* DENSITY DIDN**T CONVERGE AFTER*. I4»* ITERATIONS, HAVE'i 

200 


~ lPEll.3) 

201 

c 



2U2 

c 


REPEAT ,(ITH closer TOLERANCES 

203 


65 

lTeR=lTER-H 

204 



CALL f'EWSIARG.NARG) 

205 



ARG2(l)=9.9 

206 



CALL 2XP0hL(Fl(NC 13 , closer, MAPG.ARG.VAU'E. LIMIT. H. lER) 

207 



IF (VALUE.lT..01*ESTJG0 TO 94 

208 

c 



209 

c 


UPDATE estimate of MfII*")M VALUE 

210 



ESTsVALUE 

211 

c 


FIND ERRORS IN CALCULATED PlPAnETERS 

212 



CALL FUNCT2 

213 

c 



214 

c 


PRINT results 

215 

c 



216 



PRINT 840.1 TEO . HEVAL , lER , H .lOTH , Ft> . VAL' IE 

217 


840 

FORMAT!* ITEKATIO 'SS ■ .IS, ■ EVALUATIONS. *. 1 7, • ERROR CODE!*. 

218 


- • HALF wrCrTH;*.P5.2,' FIELD STREMGT4:*iF5.0r» SI6 HA:*iF9 

219 



II=IPCENT!SIG(11/ARG!1J J 

220 



print 84b,ARGIl).SIGUJ,II,C0VAR(l,l) 

221 


84S 

FOR.MAT!lH0»Tl3.*aFST FIT VALU'’S» ,T40, *C0RRELATI0M MATRIX*/ 

222 



- * ELECTRON JriSITY=*,lPEU. 3,* ,0PE».3. 14, , T47.F6.2) 

223 



A=ARG!2)+faASE 

224 



PRINT 84&,LABEL(2*LINES-1.) .LAHEL(2*LINES) .A.SIGI2) . 

225 



- C0VAR(1.2) ,C0VAR!2,2) 

226 


846 

F0RMAT!T?,2A6, 'CENTERS*, F9.a. • +-' ,F4.2,T47,2F6.2) 

227 



H=IPCENTt£ll)/INT!l) > 
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228 

229 


847 

200 

201 

202 


848 

203 

204 
206 
20o 
237 


850 

208 

C 


209 

240 

241 

242 


855 

240 

244 

C 


245 


70 

246 

247 

C 


248 

249 

250 

251 

252 
250 


860 

254 

255 

256 

C 


257 


92 

258 

239 


902 

2ol) 

261 

2o2 


93 

260 


903 

204 

205 


94 

266 


904 

2o7 

268 

2o9 

270 

271 

272 

273 


9S 

274 


905 

275 

276 


96 

277 

278 


97 

279 


907 

280 

281 




^ **•*•* rilEOBT «•♦••• 

p;?l,ir 3‘*7,rjT<l»itH)»IX» <COV4tUI,3),I=l*3J 

POH IAT(Tl3r*U'TF4*‘.ITV=*,h7.‘**' +-• .Fo.**. In, »«• .T‘*7,3Ft>.2| 

XI=It*l.£'iT{ei'’l/I.,T(2)) 

FHtMl O^O.IMta) rri2)rll, (C0VAK(I,4),I=i,4i 

FO'l’tAK* OACKf.ROUf'U I JTENS1TY=' ,F7.4, • ♦-* ,F6.‘I, lt>, ,T<*7, 

IF<RATI0 ,L7. I.t-5»o0 TO 70 
I r=IHC£MT ( E.HAT10 1 1 J /RAT lot U > 

PRirtT 3bO,RATIO(l) •EtTATlO(l>*n 

FORMAT tlH0,T13,»LlriE;COnT=',F7. 2,* +-* .F6.2* 16. '*» ) 

IF LINE IS HE II 4686. FIND « PRINT TEMPERATURE 
IFfAdStAR6(2»+BASF-46o5,7S» .GT, lO.JGO TO 70 
IIsIPCtMTtET£‘<P/TEMP) 

PRINT ab5,T£MP.ETEMP.II 

FORMAT ITU,' TEMPERATUHEr . . F7. 2 » ♦ +- » , F6.2. 1 6 . • S * > 
LINEP=LINEP+1 

WRITE results to file IF DESIRED 
IFtOUTJWRITEtOUTFIL* OUTRtC TRLOCK.KSMOT, LENGTH . I NTEXP 
ARo(1J=I.E-16*AR(,J1) 

PLOT NOT IF DESIRED 
IF t LATER) GO TO 20 
LINEP=LIMEP-*l 

IFt,N0T.TES/J0(»PL0T?3*.S'i6))G0 TO 20 
CALL TPl.OTtIAAX,(JoTTO".ToP) 

PpiNT 860. INAX. BOTTOM. TOP 

FORMATt* AXIS LABELS; *.F7,‘».2F6, I) 

LINEP=Llr.EP+l 
GO TO 20 

notify of lrrors 

A=BASE>ARG(2) 

PRINT 002.A.LENGTHtl). LENGTH ( PU ) 

FORMATt' line center OF.,F9.2.' TOO FAR FROM fXP RANGE'»F9.2. 
« * T0».F9,2) 

GO TO 9o 
print 903.1 

F0RHAr(T6.*tO>lLY'.r2.» POINTS PRESENT) •) 

GO TO 9b 
PRINT 904 

FORMATt' SIGF'A SOUAREO IS TOO SMALL*) 

A=1.E16*ARG(1) 

s=o. 

print 640, iter. . lEVAL.IER. value, a. S.S 
A=AR6<2)+OASE 

print g4S,LAMELt2*LIMES-l) ,LA"ELt2*LINES) .A.S.lMTtD.S 
GO TO 20 
print 906, AVG 

FORMATt* ISIGMALS AVERAGE ONLY* .E13.6. ' > * ) 

GO TO 20 

IFtOUT .and. OUTREC.lt. 13)WRlTEtOUTFIL*0UTREC)E)J0 
STOP 

PRINT 907,NF''T.TFMT 

FORMATt* INPUT DATA FILE FORMAT IS*. 17.* RATHER THAN*. 14) 

STOP 

END 



107 




.£05J7jJI 

1 

Z 

4 

5 
b 

7 

8 
9 

10 

11 

12 

13 

14 

15 
lb 
17 
la 

19 

20 
21 
22 
23 

e 4 

25 

26 

27 

28 


»,<0.iKM 1C 
UBAI<K 
C 
C 
C 

c 

c 

c 

c 

c 

c 

c 

c 

. c 
c 


r!» ( 1 1 .OJAr'K, 

PUOC 

PrtoctSjRE; NA’i-:... 
ORAiJK 


PURPOSE., , 

TO T^lAilSMtT THEORtTICAL' EXPEfttME UAL. AdU CALCULATED 
INFOH'IATIOM PAST THE LIBRARY RmJTTJE 'FMCG' 


programs using this PRCCt-DURE.. . 

FETCHS.FUi-ICTl.FUNLT2rNEWS,NEWT,NEWU*SlGPA,SUMF*THtORY.TPLOT 

Paha JETER LI>JS=2o 
PARANLTER LIMFXPS36 
parameter LI'1SH:b 
INTEGER TOTAL 
LOGICAL OPT 
real 1.JTEXP*LCN6TH 

common KS . KU , • lU 1 f N’J2 r iCShOT > TOTAL » ' 'EVAL . I TER * lER » tPSI i EPSNE # HHM , 

- NSl.NS2.F0 1 DASEp 

ALPH1ILI'‘S).SALPH1 [lI ’S) .ALP^2tUr'S) .SALPH2(LIMS) . 

CtULHS) .TOELILILG) . 

ISHOTIH •EXPJ.LEJCTHtU'FTPJ.UTfXPlLIMFXP.LINSH). 

- SUM(4,51 pERRORCM) .pi?0) pOUO) 

EQUIVALLHCE {f3RIGHT.SUH{l,31 ) . IBACKGD.SUMt 2 1 3) ) 


END 


theory (GROUPJ **♦*♦* 


205373UlNp«*iOKKSPACes< 1) .GROUP 


1 

2 

S 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 


GROUP 

C 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 


PhOC 

P^lOCtCURE pJAPF... 

GROUP 

PURPOSE... 

TO communicate PARAMETEP ERROR ESTIMATES. 
programs using THIS PROCEO'JPE .. . 

THEORY . FUNCT2. 1 JEaU 
PARAMETER ISI7E=55* ( 2+LIMSH) «LlHEXP 

COMMOl I/uROUP/ INT^'ON . T A VLR. NU .*<1 . W2 . NARG . TEMP. ETEMp • 

* LINES (10). 

- VALUE.HWIOTH. 

■< ARG(4J»SIG(4) .INTt4)(E(4).RATIO(3)»ERAT10<3) .C0VAR(4,4) 

DIMERS ICM UL0CK(55) 

EQUIVALENCE ( BLOCK t 1 ). INTPON) 
real IHTMON.IriT 


EKO 


PACB2S 
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rhEv;,1Y (rCTCHSI 


1 
i. 
:> 
4 
b 
6 
7 
6 
9 


..’ALL»U».ftrCHi, 

iOnaOUT I , J£ FFTCHS ( ARO . flARGi L I'lES ) 


NAME. 


FtrCHS 


punposE... 

TO DETFR'II'IE WHICn LI JE HA9 ^EE'! SCANHEO' flLAO IN THE 
COPRESOOHOYNO THEORETlrAL PROFILE. AND INITIAUI2E THF 
NONLINFA'’ PARAf'ETc-Rb. 


LO 

C 



11 

C 

CALLIrtS SEQUetlCE... 


12 

c 



13 

c 

call FETCHS 

(ARG.HARG.LINES) 

14 

c 



15 

c 

ARC RrO 

ARRAY 'IF HONLIMlAR PAPAW£TERS 

lb 

c 


ARG(l) IS THE ELECTRON DENSITY TIMES l.E-16 

17 

c 


ARG(2) IS THE distance FROM THEORETICAL 

IS 

c 


TO actual line CFNTER. 

19 

c 

NARO 1.0 

nU^bEP OP NO.lLItlEAR '’ARA^'ETERS 

20 

c 

LINES I.O 

LiNt CODF NO'-REK* 

21 

c 


1 hE II 4646 (DENSITY 1.E17J 

22 

c 


2 HE II 1640 

23 

c 


7 HE II 1215 

24 

c 


4 he II 1025 

25 

c 


5 HE II 4636 {DENSITY J.EIBJ 

26 

c 



27 

c 

variables in 

3LAM common: 

2S 

c 

NARG I.O 

P NOtJLIHEAR AR3UMENTS 

29 

c 

N51.M52 


30 

c 

I.O 

« ENTRIES IN arrays ALPHI. SALPHI. AND 

31 

c 


ALPm2. SALPH2i RESPECTIVELY 

32 

d 


(PptSET OOSITIONS) 

33 

c 

HHW R.O 

half MALP width PF I JSTRUMEr T FUMCTlON 

34 

c 

BASE R.O 

THEORETICAL LINE CENTER MINUS 10 ANGSTROMS 

35 

c 

ACPHl.SALPHl 


36 

c 

R.O 

AHRmYS of values of alpha « S( ALPHA} 

37 

c 



38 


iNCLUfE DBANK.LrST 


39 


DUEhSION ARGI41 


40 


INTEOlR STHEO 


41 


LOGICAL OPT 


42 


data 5THEO/29/.lFl>T/oOO/ 

43 


OEFIHE file STHE0t:0.2*LIMS+2.Li.NrEC) 

44 


read < STHEO 'UNFHT 


45 


IFlNFi.T.NE.IPHTJGO ' 

TO 90 

46 


NARG=2 


47 

c 

START 

with an electron OFNSITY of 10**17 

48 


AhG(1)=1.£-16*1.EI7 


49 


LINE5=!W<UER{*LItlF NU IBER’O* . 1.9. *30 ) 

50 


READ t STHEO 'LlnES+l JBASE . tN ALPHl . SALPHl 

51 


N51=n 


52 


N52=HS1+1 


53 


AR6C2)=10. 


54 


base=base-io. 


Sb 


ifidase-lt.sooojgo ■ 

ro 25 

56 


HHM=.3 
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ffiEviRY (FtTCHSI 


b7 lF(opn'c»))HHn=.rooi 

bu RETunr. 

b9 25 dASE=oAbE-.68 

bO 

61 iFConri 'CM )H(u=.nooi 

oZ RETURN 

63 30 STOP 

61* OO PRINT 906.tlFyT.lFMT 

o5 906 FORMATt* SlALPHA) FILE FORI AT IS*. 15.' RATHER TUAN'. I'M 

66 STOP 

67 ENO 


•I.EJCT 


teGINAE PAefS’B 
^ POOR 



ThEoHr tFl' JCT31 


^U?>J73Jlf».OkKS<*ACe,tl) .FO ICTJ 

1 FUMCTIOI FWJCTS (flHOI 

2 C 

3 C PUnpOSEf.. 

4 C TO PIK*) Try Mg'Vti.ili'AXF OCVlATt^N OF The DATA FHOh THE 

s C THEOKY U‘5IN6 CW»P£NT -W*l*MeTf«S. 

t» c 


7 

8 
9 

iU 

XI 

12 

13 

14 

15 
a6 

17 

18 

19 

20 
21 
22 

23 

24 

25 

2b 

27 

20 

29 

40 

31 

32 

33 

34 
45 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 


INCLOSE 0BANK(LIST 
IfiTEutR tRKCT/0/ 

OI.-«Ei<5XOtI ARCtlJ 

logical opt 

COMM0N/F3e0'</A (2) rOLOVAL. OLOANG 
C NEEfiM'T CALL NE'ill UNLESS DENSITY CHANGED 

IF(All).NE«ARC(l) .0R.AI2) ,NE.AHG(2)> 

- AHGLE=S7.206»ATAi.2(ARG(a)-A(a>»ABG(l)-'Ul ) 
IF<AHG(A(a)-10.) .GT,20.)tRf!CT=ERRCT+l 
IF{AI»S(AWgU)-A( 1)) .LT, l,E-7»ARG(tl)G0 TO Jr 
IF<A(11 .LT.l. .OH. A(l) ,GT.10n.)EPRCTS'’RRCT+l 

C IFfEKRCT.Gt.BJGO TO 90 

IFtBRlGHT .LT. -I.E-5)G0 TO 90 
CALL I E<TtAHG.2) 

A(1)=ARG(1) 

10 CALL t EaU(AR6.2) 

F(JNCT3=blG'.A{ARG) 

IFtAHCLt.LE.O. ) A JGLE=A’ CLErieO. 

IFtOPK'Z*) .AND. ABS(AN(,LE-OLDA’J3).OT..4) 

- PRINT 701.A.0LOVAL.OLOAHG 
1F<0PT('Y»)) 

- PRUT 701»ARGtl» rARtoJE) .FUNCT3.AM''LErDRlGHT»BACKGD 
701 F0RMAT{3Gl5.a.FS.tt»lP2G9.3) 

AI2I=ARGC2) 

0LDV«L=FUNCT3 
OUOAiJG=A JCUE 
NEVAL=N£VAL+1 
RETUh . 

90 PRlnr 703 

703 FORMAT! » ' 'FtlrlCTS' > FILob UNRFASOHAPLE ARGU'lruTS !6Ti« TtPE).') 
PRINT 701,A«0LDVALrOLD-Nb 

PHlijT 7a7,ARGU) rl'ASEr AR(,{2 I ,FO,dPIGHT.BACKGr 
707 FORMhK* NOa: ML=1.E16.'.IPG0.3.' CFN=« pG^-Sr » + ’ »G13.6, » F" 

- 09.3/* XlRIGItTriJACKoDz' »2G9.3) 

00 9£ KU=nUI.>)U? 

*X=LL..bIH(KU) 

TT=V«LUE { del p TOELp NSa r Xx-d ASE-AR g ( 2 ) p 0 ) 

Y2=0. 

HsflSHoriKU) 

00 91 KbHOT=J.H 

91 Y2=r2+INTtXP(KUpK5H0T) 


48 

49 

50 

51 


Y2=Y2/N 

92 PRINT 724,XXpY2pTT 
724 F0RMAT(1PG14.6»2G10.2J 
STOP 
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JHFuRr (FUNCT2) 


ZU!>J a: 

:£i(I).FUf CT2 


1 

f* 

SU’IHUUTI IE FU ICT2 


3 

u 

c 

NAME. . , 


«* 

c 

FWICT2 


i> 

c 



b 

c 

PUHP05E... 


7 

U 

c 

TO find THF EXoecTEO EP10HS IN The BFSt-FIT PARAMtTERS. 

O 

*1 

c 

c 

CAt-LlNG SEOUFflCE... 


10 

c 



11 

c 

call FIJNCT2 


12 

c 



19 

c 

variables in 

blank common; 

14 

c 



15 

c 

TOTAL I»1 

U CaTA POINTS 

16 

c 

EPSI R»1 

step in wavelength 

17 

c 

EPSNE R'I 

step in electron density = ARG(l) 

16 

c 

FO PfO 

HOlTSMARK field strength 

x9 

c 

3A5E Rfl 

theoretical lime center minus lu angstroms 

^0 

c 



21 

c 

VARIABLES I.‘ 

COM.*oN /GROUP/' 

22 

c 



23 

c 

NAR6 I > 1 

HUboER OF elemfnts in apg 

24 

c 

A»G Rrl 

arrat of nonlineap parameters 

25 

c 

VALUE RiIO 

input: sum of SQUARF5 OF DEVIATIONS 

26 

c 


FOR THE REST-FIT PARAMETERS. OUTPUT: SQUARE 

27 

c 


ROOT OF SUM OF SQUARES. 

26 

c 

516 R>0 

ARRAT OF expected ERRORS IN THE NONLINEAR 

29 

c 


PARh'IETEFS 

30 

c 

I. IT R>0 

array of linear parameters 

31 

c 

£ RrO 

array OT expected ER'’C°S of LIHtAK PARA IETERS 

32 

c 

COVAR R<0 

upper TRIA.>1GLE is set to a f.ORMALIZEO 

33 

c 


VARIA«E-C0VAR1AHCE "ATPIX. 

34 

c 

RATIO HtO 

line to (loo ANGSTROM! CONTINUUH RATIO 

3S 

c 

ERATIO R.O 

expected Error i»i ratio 

36 

c 

TEMP RrO 

TeMpFRATUWE FROM THE U'lFtCONTINUuM RATIO 

37 

c 

ETEMP RtO 

expected expfrihemtal error In temp 

36 

c 

HMIOTH RrO 

half itALF XIDTH OF fcXRFRIMENTAL LINE 

39 

c 


(ANGSTROMS) 

40 

c 



4l 

c 

method. . . 


42 

c 



43 

c 

THE VARIAfICE-COVAKIANCE MATRIX IS NORMALIZED bT DiVIOIMS 

44 

c 

EACH ROW AND 

EACH COLUMN BY THE SQUARE ROOT OF THe ORIGINAL 

45 

c 

0IA50UAL ELEMENT, THF CALCULATIONS OF FRATIO AND ETEHP HAKE 

46 

c 

USE OF THE APPKOXIHATELT KNOWN COVARIANCE MATRIX ELEMENTS, 

47 

c 

THE TEMPERATURE IS FOU'lU USING THE THEOPETICAL RESULTS OF 

40 

c 

DELCROIX AMD 

VOLONTE. 

49 

c 



50 

c 

SUBPROGRAMS REOUlPCD... 

51 

c 



52 

c 

NEWTrNEWLNSUMFrSleMArSTMlNV 

53 

c 



54 


OIMENStON AR€UM(4I 


55 


LOGICAL OPT 


56 


OtMENSlON TATTLEI6) 



fOOfii 
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rHEORV (PUNCT2) •*«*»• 


b7 

ae 

59 

60 
ol 

63 
03 
6<» 
dS 
o6 
67 

oa 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 
83 

64 

65 
86 
67 
83 
69 

90 

91 

92 

93 

94 

95 

96 

97 

98 
9tf 

100 

101 

102 

1U3 

104 

IDS 

106 

107 

108 

109 

110 
111 
112 
113 


DATA oI5I/10./r91CNE/X0,/ 
real INT 

IfjCLULiE Ji>A. 4'(>LIST 
INCLUDE GROUPfUST 

c GET THP ft, til. VALUES OF THE SUV OF SoUAHES d ERROR 

c hatrix elements 

CALL f!E*»T(ARGtNAKr,> 

CALL NE»‘J(ARG«NARG> 

V1=S10MA(ARG> 

XF(OPT< '0*» (PRINT 701»AR6(1>.AHG(2( »V1 

701 FORMATt*0 *r2Gi3.6.» FlNn»»6l2,6J 
Vl=FU.XT3(AO;) 

IF(0PT<*0'( (PRINT 701>AR3a)>Ah6(2) »V1 

MARKLRsO 

BEPSI=8IGI*EPSI 

8EP5nE=8I0Nc,*'rPSHI' 

00 lU 1=1>NARS 
10 ARSUM(I>=ARG(T) 

call UC»r(ARGUM>MARG) 

CALL tiEAU(ARr,'iM.UARG) 

CALL bUVF ( ARGI W t ,-IARG* 5F , SF3 . VALUE) 

TATTuc.t41=VALUE 
lMTtl)=SU,vUrTI 
IMTt2)=SU( 12.3) 

V=VALUE 

AR60H ( 2 ) =ARG ( 2 ) -SEPSI 
313SIGM(ARGU'() 

IF(Sl.LC>VALUe) -*arker=i 

TATTLL(3)=61 

ARGU'i ( 2 ) =AK6 ( 2 ) +8EPS1 

CALL 9UMF(AaG’W.NAnG.SF2*SF02.S2) 

IF < S2 .LE . VmLUE ) "<APKEa= 1 

TATTLE(5)=S2 

COVAR{2f2) = .5*(TOTAL-<*>*(51-2»VALU£+S2>/0'CPSI**2 
ARGUS ( 1 > =ARG ( 1 > +UEPSNC 
CALL I4E*T(ARG'JM»NARG) 

Si=SlC.,^A{AH6ll/) 

IFtSl.LE,VALUE>MAPKEH=l 
TATTLE (2) =S1 
AKGU'i(a)=AKG(2> 

call bUMF ( a R 61J4, t JA RG*SF3,SF')3»S3) 

IF t S3 . LE . VALUE » MAOnER= 1 
TATTLE tl)=S3 

COVAK{1»2) = .5*(TOTAL-4)*( V ALUE+S 1-S 2-53 ) / < OEPS I *BEPSNE ) 

ARGUM U ) = ARG ( I ( -UC^PSNt 
call I.E*(T(ARG'IM.NARG) 

IF(SI.le.valuei4arker= 1 ■‘■CALL NEWT ( ARG,. lARG) 

TATTLE(6)=SI 
IFlMAKKER.eO.DPHINT 702 

702 format (' »*«** LEAST SQUARE SOLUTION NOT FOUND. . .NEAR0T VALUES', 

- • FOLLOW **••*•) 

IFlHAhKER.EQ.l .OR. OPT( 'O' )) PRINT 707, TATTLE 
707 FORMAT (13X, 2013, 6/3G13.6/13X,Gl3.6/) 

C0VAH{1*1)=.5*( TOTAL-4) *tSl-2.*VALUF+S3)/8EPSNE**2 
C0VA«(1'3)=( SFD3-SF0 ) /BEPSNF 
CO V A R { 1 » 4 ) = C SF3 -SF ) /HEPSnE 
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T-HEoRY ('•I'NCT?) 


1 14 

llC) 

U7 

uu 

119 

lai 

122 

12} 

124 

125 
I2b 

127 

128 
129 
IJO 
lil 
1}2 
1}3 
U4 

135 

136 

137 
135 

139 

140 

141 

142 
14} 

144 

145 

148 
l47 

146 

149 

150 

151 

152 

153 

154 

155 
15b 
157 

156 
159 
IbO 
161 
162 

163 

164 

165 
186 
lo7 
loB 

169 

170 


C 

c 

c 

c 


c 

c 

c 


c 

c 


c 


c 


6 


c 


COVflrt (■>*})= I ^rO->-SF3)/^EKSl 
CoVAIi C2*4 1 = 13F2-JF) / Jc.J SI 

C01/AR13.3), C0VA«I3,4)» AND C0yAR<4»4) ARE SET RT 
PROGRAM •Nt'IU* 

iNi^FRr TO FING the vARI ANCE-COWARIANCE MThIX 

00 15 J=1.4 

15 IFCOPTt'Q' ) JPrMT 715f ICoVAR(lrJ)*I=l»J) 

7l5 FC«NAmx»401}.61 

CALL 5TNIflV(COVAR»4»4»0»ib01 
00 lU J=lr4 

la IFtOHTl'OMJPRlNT 715f [COVAR(I»J)tI=l»JJ 
00 20 JSli4 
00 20 1-1>J 

20 COVAK(I*J)=COVAR{IrJ)*V 

RESCALE ELtCTROH DE4SITV 
ARu(1>=1.E16«ARG(1) 

CORRECT RFFERFNrFS TO DFflSITY 

DENSITY actually VARIES LIKE (LINE »I0TH1«*1,2 

R=l. 

lF<Ab5{ARGf2)+BASF-46G5,7S} .GT. lO.JGO TO 75 
R=T 1.^/1. 3) »tARGln /I. El 7) *»( 142/1.5-1.) 

Afl5(l) = l.E17*( «APG(l)/l,tl7)»«<].2/1.5> ) 

WE AUST recalculate THE HOLTSMARK FIELD STKENCTH 
WITH THE lltW rENSITY 
F0=1 .w 503E“9*ARGU ) *♦ .66o66b6687 
25 COVAHH»l) = l.E32*COVARU.n»R**2 
COVAn(l«2l = l.El6»COVAf (lr2)»<’ 

C0VAR(1^3)=1.E16*C0VAR<1,3)*R 
COVAR 1 1 »4)=1 .E18*C0VA« ( 1 ,4 ) *R 

EXTRACT The variances of THF ITtnrviDUAL parameters 
SIG(l)=5I6MtSnRT(AaS(C0VAR(l'l 1 1 ) .COVARTlt] )) 

SIG(2)=SIGN{SoRTt ADStCCVAR12.2) I ) ,C0VARt2»2)) 

£{!)= 5IGN(S0RT(ABS(CCVAR{3>3J ) ) .C0VAR(3»3) ) 

E<2)= 5rcN{SORT(A3S(C0VrtPl4.u) ) ) .COVAR t 4.4 1 ) 

FI JO THE LIVE'CO’ITI’lUUM RATIO < ERROR 
aATio=iNTa)/(iNTt2)«ioa.) 
fcRATIiy=RATIO*SQRT(AaS( C0VART3*3 )/IWTU)a*2 

- +COVAR(4.4)/irJTt2I«*7 

- -2.»C0VARC5p4)/ni.T(l)*I'II(2)) 1) 

IF line is HF II 4645* FIMD TF' pFRATUHE a ERROR 
IF(RATIO.LT..001 .or. ABStARS(2l+nASE-4fea5.75) .GT.lO.lGO TO 30 
TEMP=.208»ALOG(ARG[1) )+,209*ALOG(?ATIO(1) )-5.10 
TT =.208/ARG( 1) 

TL=.209/INTT1> 

TC=-»«-09/If!T(2) 

£TE«P=SQRT TABS! TIJ*T’l*COVARtl*l) 

- ♦TL*TL*C0VAR(3*3I 

- +TC*TC*C0VAR(4.4) 

- +2.*TN*TL*C0VARUt3) 

- +2.*TI.*TC*C0VARtl>4) 

- +2.*TL*TC*C0VART3.4) )) 

N0RHALI2E the VARIANCE-COVARIANCE MATRIX 
30 C0VAR(l*2)=C0VAfi(li2)/JSI6<l)»SlGT2) > 

C0VAR(l*3l=C0VAR<1.3)/(S16t 1)*E<1) ) 

C0VAR(1*4]=C0vART 1.4)/(SlG(t)*E(2>) 
C0VAR(2»3)=C0VART2*3)/(SI6(2)*ET1) ) 




114 


171 

17 ^ 

173 

1/4 

175 

176 
1/7 

1/e 

179 

160 

161 

162 

163 

164 
185 
166 
187 
168 

189 

190 

191 

192 

193 

194 
19b 

196 

197 


• •••** IFUUCT2) 

COVArt(2*4>=COvA'U?,4»/tSI612J*bl2l) 
C0VAR(3»4)=C'’\/AB(3*4)/tEU)*CI2J 1 
00 40 J=t»4 
40 COVAt<( J)=l. 

C FlflO half IHTENSITY POI'lT OF EXPERIHtMTAL PROFILE 

HALF=T0EL{2)/2. 

00 4b Kb=2«NS? 

IF<TJLLtKS).LT.HALF)00 TO 4S 
45 CONTINUE 
KS=Nb2 

48 X1=DLL(KS) 

Tl=TJEL<KS) 

X2=0£L(KS-1> 

T2STU£L(KS-1» 

00 5U US1»4 

XS(<T2-rlALF)*Xl-(Tl-f<ALF>*X2>/tYa-Tl) 

X2SX1 

T2=Ti 

X1=X 

so Yl=VALUEtOeL.TDEL,NS2.X,uJ 

HW IDT6= ( I Y2-H4LF J *X 1- ( 1 1-HALF ) *Xa ) / ( Y2-Y1 > 

VAUUt=SQRTtV) 

RETURN 
60 print 908 

903 FORMAT I • IilVEflSION FAILURE FIhUING COVARIANCE MATRIX*) 

RETURN 

END 


U.EJLT 
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1 

2 

5 

6 

7 

8 
9 

XO 

il 

X 2 

13 

14 

15 
lt> 
17 
la 
19 
2 U 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 
3d 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 
46 

49 

50 

51 

52 

53 

54 

55 

56 


TIIEO^r (T>"LOT) ****** 


.OiM*jrACtKl>»r»»coT 

subroutine TBlOT t VBAX f POTTOM, TOf> I 

c 

C 

c tplot 

c 

C PURPOSE... 

c ro PLOT THE EXPEni'lENTAU DATA A*1D THEOHETICAL PESt- 

C FIT PROFILE. 

C 

C USAGE... 

C 

C call tplot (AMAX. bottom. topi 

C 

c AMAX fi.O LABtL FOP EMO OF Y AXIS 

C bottom fi.O LAStL FO® ORIGIN OF X AXIS 

C TOP R.O LAatL FOR END OF X AXIS 

c 

C variables in common; AL'-’OST EVFFTTHINO 

c 

c subroutines used... 

c 

C PLOTC standard plot suhroutime. used to position the '■EN 

c NSCALE FI, JDS PLOT SCALIf'G PARAMETERS FOR EASILY INTERPRETED 

c AXIS labels 

C AXISN draws All AXIS 

c 'STMBOL' draws a jYMROL AT THji hesIRSO POSITION 

c PAGEUP COMPLETES THE PLOT AND P^'ROSITIOHS THE PEN ONTO 

C the NEXT page 

c 

c method... 

c note That the axes used in this routine are rotaTlD so 

C DEGREES CCH FROH THOSE USED BY The SYSTEM ROUTINES. THUS. 

C HY t.Y* OIRECTION IS THEIR .-X'. AND MY >+X* IS ThEIR •+Y*. 

c 

DIMENSION ARG(4( .INTI4I 
INCLUDE OBANK.LIST 
INCLUDE GROUP. LIST 

real l.int.lbound 

data UlDTH.HEIGHTr NY( NX 
" / S.» 5. Sr 5i 15/ ^ 

C 

C BEST FIT CURVE 

FIT(Z)=IffT(l)* VALUE (0ELrTDELrMS2.Z.0) + lNnNAPG> 

C 

c find largest S SMALLEST VALUES ALONG EACH AXIS 

aOTTONsLEiTGTHTll 

T0P=LEN6TH(NU) 

TmAXsFITTO.) 

CENT£R=AR6t2)+BASE 
UBOUND=CEtjTeR+X O..HWIDTH 

lbound=center-io.*hwidth 

DO 15 KU=1.JJU 

C DON'T PLOT POINTS FURTHER THAN lO HHW TO EITHER SIDE 

IF ( LENGTH IKUJ.GE.LOOUND) aOTTOM= AM IN1(80TT0M»LENGTHIKU)> 
IF(LENGTH{KU),LE.UOOUNO|rOP=AMAXI(TOP.L£NGTH{KU)> 

N:NSH0T(KU> 
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• **••• TMEu»^Y ^►'LOT) •««*«* 


t>7 

b<} 

OQ 

61 

62 

63 

04 

05 
66 

67 

68 

69 

70 

71 

72 

73 

74 

75 
7o 

77 

78 

79 
OO 
81 
02 

03 

04 

85 

86 
87 
aa 

89 

90 

91 

92 

93 

94 

95 

96 

97 
90 
99 

100 

101 

102 

103 

104 

105 

106 
lOT 

lua 

109 

110 
111 
112 
113 


00 16 k3ffOT=l.H 

15 YMAXSAMAXI IT'WXt IWTEXt’fKU.rtSHnrn 

C PLOT AT LEhST 4 Hhw TO EACH SIDE 

BOTTO =A IJHl tnOTTC'ltCLdTtH-A.^MWIjTH) 
TCPsAPAXl(TOP.CEHTER*4.»riXIOTHJ 
C 

c PIHO SCALING PARAMETERS 

NTICX=NX 

call flSCALE(80TTOM.T0P»NTICXrWI0TH.DX) 

DX=1./0X 

NTICT=ny 

ZER0=0. 

CALLNSCALE(ZeR0.VMAX,NTICY.HEI6HT,0Y) 

OYs-l./DY 

C 

CALL PL0TC(7.5»1.0*-31 
C 

C DRAW AXES 

CALL AXISN{0.,0.»HTICY*Hc16HT.“, 1,180. ) 
call AXlSN(0.,0.,NTICX>WiaTH,.l,9n.) 

C 

C DRAW A COTTED LINE FOR THE CONTINUUM LEVEL 

C IIF POSITIVE) 

IFUHT(NARG).lT.O.)GO to 20 
UxOYAiNTtfJABG) 

L=DX» < TOP-BOTTOM ) 

DU=-L/69, 

DO 18 KU=1,35 
CALL PL0TC(U,t.,3) 

CALL PL0TCIUrL+DL,2) 

16 L=L+UL+DL 

CALL PL0TC(U.L'3> 

C 

C PLOT THE EXPEBI’-'ENTAL POINTS 

20 NSYMi1=4 

DO 30 KU=1»NU 
1F<KU.EO.NU1)nSYMP=0 

IFJLtMGTHIKU) ,UT . bottom ,0P. LENGTH tKU) .6T, TOp) GO TO 30 
U=0. 

N=NSH0T(KU) 

DO 26 KSHOT=l,N 
25 U=U+It.TEXPlKU,K5H0T) 

CAi4.SYMO0L(Oy*U/N,DX*(LEMGTH(KU) -BOTTOM), .l.nSTPB, 90. •-!) 
30 IF(KU.CQ.NU2)MSYMB=4 
C 

C DRAW THE THEORETICAL CURVE 

L=TOP 

DL=- < TOP-BOTTOM) /250 . 

NPENX3 

00 40 KU=l,2Sl 
0=FIT(CE71TCR-L) 

CALL PLOTC (DYaU.OX* IL-BOTTOH) ,HP£M) 

NPEN=2 
40 LSL+DL 

call PAGEUP 

RETURN 

£ND 
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THtOPY (AiaS'l 



.Orf^ 

GPACti.* 1) .AXIStl 

1 


SUlHObTriE AXISN (XX»YTiHTIC«AlHTri»TICiANSLe» 

2 

c 


3 

c 

NAME.>> 

4 

c 

AXIS * 

5 

c 


o 

c 

PU«PObE.. . 

7 

c 

TO URA9 ONE AXIS fOH A 2-OIME()SIONAL GRAPH 

a 

c 


9 

c 

US ACL* *• 

10 

c 


U 

c 

call AXISN <XrTiNnc»ALNTH.TIC»ANGLE! 

12 

13 

c 

c 

X»Y R»I POSITION OF START OF AXIS (INCHES FROM PAPER 

1<( 

c 

ORIOIN) 

IS 

c 

^ NTIC If I ’.UltoES OF LI'E S<^GME’ TS I pTaEEN TIC MARKS 

16 

c 

ALNTH Rfl LEtJoTH OF AXIS (1HC1IES» 

17 

c 

TIC R»I length of TIC MARKS (I’JChESJ 

IS 

c 

TIC.GI.O FOR NARKS OJ CLOCKWISE SiDt 

19 

c 

TIC.LT.O FOR harks ON COUMTLP-CL0CK*ISE SIDE 

20 

c 

ANGLE R*I angle OF AXIS FROM X~AXIS (DEGREES) 

21 

c 


22 


a=alhth/ntic 

23 


x=xx 

29 


T=rY 

25 


C=C0Sl.0l79533*AilGLE) 

2b 


S=SIN( .017953T*AI1GLE) 

27 


ax=c»a 

28 


AY=S»h 

29 


TX=-S»T1C 

30 


TY=C*T1C 

31 


CALL rLOTCtX+TXfY+TY/3) 

32 


DO 10 I=1*NTIC 

33 


CA>-t PL0TC(XfYt2) 

39 


X=Xf^AX 

35 


VsYfAY 

36 


CALL PL0TC(XrYf2) 

37 


10 call PLOTCtX»TXfY+TYra) 

38 


CALL f LOTC(X*TXfY*TYf 3) 

39 


RETUKtJ 

40 


END 


H*£dCT 
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THEu>?t 


iJUb37jolH*,OHKil’Al.E,i < I > V5 


1 

2 

3 

4 

5 
b 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 
10 

19 

20 
21 
£2 

23 

24 

25 


C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUrtHOtrl 111 H£'fS( APOrilAHOJ 


4FrfS 


FOHPObCo. 

TO CONVOLtf*: THE THEORETICAL Ain THE INSTRUMENT PRuULES. 
CALLINS SEOUEmCE*.. 

call NEHS(ARSiNARu) 

ARC R»I ARRAY OF NONLINEAR PARAi'fTEKS 

ARGIl) IS the ELFCTRON DeMSITY THES, 1,E-16 
NAR6 It NUMBER OF NONLINEAR PARAMETERS PRESENT 


method. 


the HOLT2.4ARK FIELD STRENGTH IS CALC'lLATEUt ANC UsEO TO 
find the AMOUNT OF BROADENING NEEDED IN (S(ALFHA) .ALPHA) 
space, this is ACv.URAT'^ if THE final FlECTROM DENSITY IS 
close to The original estimate used here, the CAUsSIAN 
HERHITE OUADRATURc formula used herf is exact for the 
INTLGRAL of a GAUSSIAN InSTRUMFNTAL PROFILE AND A FIFTH ORDER 
curve pop t«e theoretical dROFTLE. 


2o 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 
33 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 


INCLUDE ObANK.LIST 
OIMDHSION ARG(4) 

C WE USE the M0LT7 'ARK mOR"AL FIELD STRENGTH 

C F0=2.603*Z«4.8033E~10«NE**I2/3) 

FO=i-2S03E-9*(l.H6*APS(U )♦•. 66686666867 
c WE USE the CROINATE FOR THE 3 PniNT GAUSSlAN 

c hermite ouaciatuke formula 

C DELTA=(HH«/FO)*50RT(3/(P*ALOG12) ) ) 

OELTA=1.471068S»HHI»/F0 
N=NS1 

00 10 kS=2.N 

ALPH2 ( kSki ) =AlPH1 ( KS ) 

10 SALPHt (K5+l)=. 166866606* 

(VALUE(ALPHl.SALPHiMl51t'LPMl(KS)-D£LTA.f,) 

+4.*SALPH1(KS) 

+VALllEtAi_PPl.SALPHl .NSl.ALPHKKSJA-DELTA.Ol ) 
JVLPH2(1)=-ALP‘I2CJ) 

ShLPH2(1)=SALPH?(3) 

SALPM2t2J=.6666668666*SALPHl (II 

+.3333333333»VALU£lALPHltSALPHI, MSI . delta. 0) 
ALPH2(2)=0. 

RETURN 

ENO 


H.EJCr 
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•••*«* theory (NEwT) *••«** 


» (OtM\.SI A.-El>l I) 


1 



Si/JHOUTlfJE rJEWTCASGiM'iRG) 

2 

C 



i 

c 


NAME • . . 

<♦ 

c 


MEwT 

b 

c 



& 

c 


purpose* •. 

7 

c 


TO calculate TiIE POOPILE in MAUELL4GTH SPACE* KflOwlNG 

a 

c 


THt ELECTRON DENSITY ANU THE PROFILE IM tSt ALPHA) *ALPHA> 

9 

c 


SPACE 

10 

c 



11 

c 


CALLING SEGUENCE... 

12 

c 



Ij 

c 


CALL NEWT (ARGiNAKG) 

14 

c 



15 

c 


ARG R*I APHAf OF NOIlLINtAR ARGlP'ENTS 

16 

c 


AHGU) IS ELECTRON DENSITY TIMES l.t-16 . 

17 

c 


NARG I*I NUM6ER OF NONLINEAR ARGUMENTS 

18 

c 



19 



INCLUDE ODANK*LIST 

20 



DIMENSIOJ ARG (5) 

21 

c 


WE USE the H0UT7'*APK NORMAL FIELD STRENGTH 

22 

c 


F0=2. 603*2*4. fl03E-10*NE** (2/31 

23 



]F(ARG(1) .LE.O.IGO TO 90 

24 



F0=1.2503E-9* ( l.E16*ARG( 1) ) **,o6oo6669 

25 



F0INV=1./F0 

26 



N=NS2 

27 



DO 20 KS=2*N 

28 



Ti)fcHRS)-F0I’JV*SALPH2(K5) 

29 


ao 

0EL(KS)=F0*ALPM2(KS) 

30 



T05L(1)=T0EL(3) 

ii 



DEL<n=-DEL[3) 

32 



RETURtt 

33 


90 

print 9U1*ARGU1 

34 


901 

FORHATC HEMT: density 0F**G9,3r' *1.E16 ’•] 

35 



STOP 

36 



END 


U.EJCT 
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I 

i 

i 

<* 

•j 

6 

7 

6 

9 

10 

11 

la 

i;s 

14 

15 
la 

17 

18 

19 

20 
21 
22 

23 

24 

25 
28 

27 

28 

29 

30 

31 

32 

33 

34 
3b 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 
4b 

47 

48 

49 

50 

51 


THEoHY (MLaiI) •••«•* 


, ;k<\ 
C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


'jPAC£»(l).*Jtyi) 

NAPE... 

PORPJbE... 

TO SOLVE the linear PAPT OP THE LEAST SQUARE FIT. 
USAGa. . * 

call NEVtUtARG^NARo) 


ARG 

NARC 


R»I 

I>1 


arraT of wonlimear arguments 

NUMOER OF nonlinear ARGUMENTS 


variables in common: 

SUM RfO 
CO VAR R>0 


SUM(1>3) IS THE LINE INTEHblTY 
5UM(2,3) IS THE BACKGROUMO INTENSITY 
COVMROt’Of C0VAPl3»4). and C0VAR(4,4> APE 
IfiStRTFO POR LATER USE RY PROGRAM 'PUNCT?' 


SUBROUTINE NEPUCAPCUM.NAKSUM) 

INCLUDE DDANK.LIST 
INCLUDE GROUP. LIST 
OIMENblOM ARGtlM(4) 

A=ARGUM(2)+BAS£ 

DO lb I=1.NAR0UM 
DO 13 J=I.iIARgUM 4'1 
15 SUMtl.JI^O. 

00 2b KU=NU1»NU2 

T=VALUE ( DEL . TOEL • ' 152 » A-Lt(JGTH ( KU ) . 0 > 

N=NSIIOTIKU) 

SUM(1 . 1»=SUM (1.11 +T*T*fJ 
SUM(1.2)sSUH( i. 2)+T»N 
DO 2b KSHor^l.N 

S'JMU ,3)=5UM( 1 . 3) +IMTEXP (KU.KSHOT) »T 
SUM(2»3)=SUN(2,3)+INTEXP(KU.K5H0T) 

25 CONTINUE 

5UM(2.2)=T0TAL 

C SAVE ELEMENTS FOR CALCULATION OF VARIANCE- 

C COVARIA- CE NATPIX 

C0VAR(3*3)=SUMU.1) 

C0VAR(3»4)=SUM(1.2J 

C0VArU4f4)=SUM(2i2) 

C have this SY.r-SETKIC SYSTEM SOLVED 

C all S Y 4SLV ( SUM >S)JM(l.NAHGUMtl)»MARGUM.4tO.S40) 

RETURN 
40 PRINT 900 

908 FORMAT (• singular MATr’IxIM 
STOP 
END 


a.EJCT 
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Ti'Ev/jr (sio'',u 


rjoi '•*^o-<xsPA- 


1 

2 

3 

4 

a 

6 

7 

a 

9 

10 

11 

12 

13 

14 

15 
Id 
17 
la 

19 

20 


C 

c 

c 

c 

c 


20 


ElU) .‘iti.vA 
KUHCTJON SI6 'A(A^(a) 


PUHPOiE... 


TO Fl'n TMF HEAfJ S'JljAR'" 
VALUES PRaPICTEJ USi-JS 


deviation of The data froh 
The co»io''nT parameters. 


Include OfjAtJK.Lisr 

Ol'IEJSIOfl AHRfitJ 

LOwIuAL opt 

SIGHASO. 

A=ARo(2) +BASE 
DO 20 KUSNUI.NU2 


iiv«i,UE.lOiiL.T4EL.f!S2.A-LtN3THtKU) 

U:SU4(l»3)*T+SUMt2,3) 

N=.-JSH0TIKU) 

Do 20 Kbl^orsl.rj 


rO) 


SieHA=SIGhA»t INTE*P( km»asH0T)-U)*»2 

SIGiOA=SIGrfA/<TOTAL~4J 

KETUKN 

Ef<D 


THE 


W.EJCT 
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THEwBY (STMSLV) 


2ab3/3^IH*«.0HKbPACCS< 1 > .Sr -<OLV 


1 

2 

i 

s 

6 

7 

8 
9 

10 

11 

12 
13 
19 

15 

16 
17 
16 

19 

20 
at 
22 
23 

29 

25 

26 

27 

28 
£.9 

30 

31 

32 
03 
39 
35 
3t> 

37 

38 

39 

90 

91 

92 

93 
99 

95 

96 

97 

98 

99 

50 

51 

52 

53 
59 

55 

56 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUiJBOL'TlrC 5YMSLV(A>0>NrMN> ITri) 

NAME... 

SY.^METR1C UINEAH tOUAHON SOLVEP 
CODE NAME... 

5YHSLV 


PUiiPUbE. 


TO SOLVE A linear SYSTEP AX=B WHEN THE MATRIX A Is 
SYMMETRIC AND POSITIVE DEFINITE. THE ROUTINE CAN BE 
called subsequently to perform The solution for 

A new PIGHT HA.iO SIDE WITHOUT OFCO’^POSING AGAIN. 


CALLING SEQUENCE... 


call STHSLV(A»BfN«NN.IT.S80) 
arguments on ElTpTJ 


B 

N 

NN 

IT 


S60 


matrix of coefficients, since it is symmetric. 

ONLY elements MFFOEO ARE A(I.J). l.LE. I .LE.J.LE.N 
ARRAY OF ELEMENTS FROM RIGHT HAwO SIDE 
DIMEf.SIOM OF MATRIX AND 8. 

maximum NIIhOER of rows in a (FIPST DIMENSION) 
SWITCH. ..n=l IF MATRIX A WAS OECOI-POSEO On A 
PREVIOUS CALL TO SYHSLV. AMD ONLY THE ARRAY B 
IS DIFFERENT THIS TIME. IT.NE.l IF A IS NtW. 
CONTROL will be passed TO THIS STATEMENT IF 
A PIVOT ele^'Emt is found of absolute value 
LESS THAN li£-m. 


arguments on return: 

A original matrix is destroyed. LOWER triangle holds 

LOWER TRIANGLE OP MATRIX L. (DIAGONAL ELEMENTS OF L 
ARE I'S.) diagonal elements HOLD MATRIX D. 
a SOUJTIOil ARRAY X. 

METHOD, .. 

symmetric factorization is USED TO FIND A LOWER 
TRIANGULAR MATRIX L AND A DIAGONAL MATRIX 0 SuCH 
THAT A=LDU, WHERE U IS L TRANSPOSED, ThE UNKNOWN 
VECTOR IS calculated 3Y BACK SOLVING THESE TRIANSULAR 
systems; U2=B . DY=Z * LX=Y . 

dimension B(5).A(25) . 

IF(N.GT.I)GO to 10 
B(1)=B(1J/A(1) 

return 

10 IFdT.EO.DGO TO 26 
DO 25 K=l.fl-1 

IF(AaS(A(K+N.J*K-NH)> ,LT. l.E-10)fi£TURN 6 
DO 25 JFK+l.N 



S>7 

39 

6» 

bl 

b2 

bi 

b* 

u5 

Ob 

dr 

ok 

69 

ro 


iMrour (<;rns;.v> 


00 2fl I=J,JI 

2« 00 30 j=2irj 
00 30 

30 t5( J) =IU J J -A I J+n»I*l-HN 1 *8 « T 1 
30 40 J21,M 

00 50 J-SN-1/1,-1 
DO SO 

50 D(31Sb(01-A<t'»rjN«j.NN)*8(|) 
FlETuftj 

EiMQ 


«»E4CT 
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THEoRt IVfiLUFJ 


<!U> 37 ^th«^ 0 KK 5 PACC} ( 1 > • VALUt, 

t FUNCTION VAMJE(X*T*NiXB»f'» 

2 C 

3 'C NAH£>*> 


4 

5 

6 

7 

8 
9 

10 

u 

1 2 

13 

14 

15 
lb 
17 
la 

19 

20 
21 
22 

23 

24 

25 

26 
27 
2a 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 


C VALUE 

C 

C PURPOSE... 

C TO IHTrSPotATE IN * TACLE TO FtNQ INTENSITIES FROM 

C THE theoretical line PROFILE. 

C 

C CALLINO SE 0 UENCE..« 

C 

c call VALUE(WAVErIf.T.NiHANT,P) 

C 

C XAVE 

c 

c int 

C N 

c xant 

c p 

c 
c 
c 

C METHOD... 

C dSTONO THE END OF THE TABLE. A WZ POWER LAw IS USED TO 

C EXTRAPOLATE. XIThIN THE TAPLE. AIT<F‘I*S PROCEuURE IS 

C applied USIiJO 4 pyl JTS. SINCE AN EVEN NUMBER UF POINTS 

C IS USED. THE IMTEKPOLATING FUNCTION IS CONTINUOUS, 

c 
c 

DIMENSION X<N).r<N).P(4) 

EOUI VALENCE tSAVEJ.J) 

XBAR=A 0 S(XU) 

C IF ON FAR Wine OF LIME. USE ASTpTOTIC FORMULA 

IFtXliAR.GE.XlMUGO TO 90 

C RETRIEVE POitJTER FROM LAST CALL 

SAVEJ=P( 1 ) 

C ENSURE 1 ,LE. J .LE. N 

J^AX 0 <M 1 N 0 (J.N).I) 

C DECIDE \HETHER TO SEARCH UP OR 002 M 

iFUUAR-xumo.an.ao 
c search down 

10 IFU.LE.DGO TO SO 
JSTARTcJ 

00 12 J=OSTART. 2 .-l 
IF (XUAR-X (d -1 m 2 . 75.40 
12 CONTIt.UE 
3=1 

GO TO 50 

C SEARCH DP 

20 IFU.GE.NIGO TO 30 
JSTAHTsJ +1 
DO 22 J=JSTART.N.I 
IF(XBAR-XU) > 40 . 80.22 
22 continue 
30 J=N -3 
GO TO 50 


ft. I ABRAT OF rfAVFLENSTHS (niSPLACEMENTS FROM 
LiNc. CENTER) 

R.l array of corresponding line INTENSITIES 
I. I lUMdER OF ENTRIES IN WAVf OR INT 
ft. I AAVbLENGTH at which INTENSITY IS OESIREO 
R.IO fcORK array of length 4 . ptl> IS USED TO STOKE 
A POINTER BETWEEN CALLS. SO EACH SEARCH OF 
WAVE REOINS 4 HERE THE PRECEDING SEARCH ENDED. 
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57 

M) 

M 

60 

61 

t>2 

65 

o6 

65 

o6 

67 

bS 

69 

70 

71 

72 

73 
79 

75 

76 

77 
7S 
79 
60 

31 

32 
63 
39 

65 

66 
67 
66 
39 

90 

91 

92 

95 

99 

95 

96 

97 
9» 
99 


•****« THEORY (VALUE) 

C SAV*^ THIS HOINTEK 

60 SsSA^EJ 

C SET J TO POl.lT TO FIRST OF THE if POINTS 

c IN X nearest XRAB 

J=AIi|0(MAX0(J-2.1),N-5) 

c apply AITkEN'S PHOCEDUHE USINS 4 POINTS 

50 P(l)=T(JI 

C THE GROUP OF STATEMENTS TO FOLLOW IS EQUIVALENT TO: 

C 00 60 I=2f4 

C P(I)=YtJ»I-l) 

C 00 60 L=2.I 

c Ptr.) = (P(L-l)*<XU+I-U“XUAR)-PlI)*(X(J»-L-2)“XBAn))/ 

c - UX(J*l«ll-*nA«l-IX(JtL-2)-XBAR)> 

C 60 CONTINUE 

c ..•BUT (WITHOUT LOOP CONTROL) WILL EXECUTE FASTEN 

P{2)=TCJ+2-l) 

P<2)=(P<2-I)*tX(J+£»i)-XbAR>-P(2)*(X(J+2-2)-XflAR))/ 

- l( X ( J+2- 1 ) -*e AB > - 1 X ( U+2-2 ) -XB AR ) ) 

P(3»=Y(J+3-l> 

P(3) = (P(2-l)*(X(JA3-l»-)UlAR)-Pt3)*(X{J+2-2)"XBARn/ 

- ( (X(J-^3-t>-XBAR)-(X(J*2>2)-XRAR) I 

P(3) = (P(3-l)*(X(J»3-l)-Jta«)-P(3)»(X<JT3“2)-XBARn/ 

- ICX(j«-3-I)-XBAR)-(8<J»3-2>-XPAR)1 
P(4)=Y( J>4-I) 

P(4) = (P(2-ll*(X(J4-4-l)-«oAR)«Pl4l«tXU+2-2)-XpAR) >/ 

( (X^J+4-l)-/^AR)-lKU♦^-•2)-X^ARn 
P(9) = (P(3-ll*(X(Jf4-l)-X6Ari)^I9J.{XlJ*3-2)-XBAR))/ 

- ^ (XlJt4-l)-ABAR)-lXlO»3-8)-XFAR) ) 
P{4»=(P(4-l)*fX(J+4-I)-XUAR)-Pl4)»(X(J+4-2)-XBAR) )/ 

- ( <X(J»4-l)-XeAR)-tX(jT4-2)-XnAR) ) 

VALUES? (4) 

c 5AVF THE POINTER IN THE WORK ARRAY 

P<t)=S 

return 

c XQAP APPEMiS IN the table. So WE USE THE 

C CORRESPCflDI^■'; table EHTPY 

75 J=J-1 
80 VALU£=Y(J) 

P(1)=SAV£J 

RETUNr. 

90 VALUt=Y(tl)»(X(N)/X8AR)**2,5 
RETURil 
ENO 


e.EOCT 
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